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| and, in fact, has been applied successfully to such 
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(problems. This special case is given below as the- 
orem 1, since its proof is simpler than that for the 
‘general theorem and exhibits the transformation 
jeans. matrices needed for the reduction formula. 
>. _This result is used in section 3 to give generaliza- 
arch, | tons of some theorems by A. Brauer ? on stochastic 
Arctic} matrices. Other applications are given in later 
sections. 
Radio’ As in the previous paper, for a given partitioning 
Gon }of a matrix A we shall call the submatrices, A,,, 
the blocks of A and we shall write A=(A The 
= elements of the blocks will be denoted by a4’, 1.e., 
A (a he 
pFre| 
fodu-} Vnless otherwise stated, the matrices will be 
aes. F. Haynsworth, Reduction formulas for partitioned matrices (to be 
i ‘A rer, Limits for the characteristic roots of a matrix, IV, Duke Math 
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v in physies and chemistry, and the reduction 
cessfully practical problems In addition, 
ces and of J. Williamson in partitioned matrices 


inh 


complex matrices. Also, since we will 
be dealing throughout with matrices A=(A,,) of 
form (1), we will assume, unless otherwise indicated, 
that the statements and formulas given are true 
for i.) > aes 


arbitrary 


2. Special Reduction Formula 


Kor the sake of completeness we include a lemma 
from the previous paper! which is needed in the 
proofs of theorems 1 and 3. In this lemma we 
consider three cases where there are zeros in con- 
venient places in the blocks of a partitioned matrix, 
A. In each of these cases A is reducible, and the 
proof consists merely in defining the permutation 
matrix which puts A into the reduced form, 
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ai 
(2) 
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where Cand 7 are square matrices and @ is a matrix 
composed entirely of zeros. 
LEMMA: Give na partitioned matrix A of orde r N 


with hiy hloel s A 


1. If n (7 } ero | 
are lower triangular with elements, h 
on the diagonal, A is similar to a matrix A, with 
blocks An,=(S"), kh 12,.... 7, on the diagonal, 
and zero blocks above the diagonal (this case is give n 
by Williamson 

2. If 


i 


and the blocks Ay 


' wee, * 


di 


C 0 
A , (3) 
D T’ 
where all matrices T,, are square, of order r, and the 
matrices C;; are (ny—r) X (nj—r), then A has the form 
(2), where C=(C,,) and T=(T,,). 
3. If Ay has the form (8) where the matrices Ti, 


are lower triangular, we will say Aj; is partially trian- 
gular. Then if A has blocks Ay which are all par- 


tially trianqular with submatrices, T,; of order r, having 


elements t\?. h - Pee r, on the diagonal, tr 
roots of A are roots of the r matrices, (t\? - 

J. Williamson, The latent roots of a matrix of special type, Bull. Am. Math, 
Soc. 37, 585 (1931 








The proofs given below would follow in a corre- 
sponding manner if all blocks were transposed 

PROOF: 

1. The rows and columns of A should be arranged 
in the order 


1, n+1, 2n t--l)n-+1 
2.n-+2, 2n-+? f l)n-+2: 

} 
n, 2n, 3n, ee th 


Then the new matrix A will have the matrices 
Am=(\}?) on the diagonal and A, 0 for i<y, 
so its roots are the roots of the n matrices, A 


K 
2. Let >\n N,. Then N n,. NN If we 


nd 


i=l 
arrange the rows and columns of A in the following 
order: 
| N 7 
\ L. N.t+2 \ ] 
\ .. 2 \ } 5 


we have a new matmx A in which the matrices 
Cy are together in the upper left corner, and the 
matrices 7';,; are together in the lower right cornet! 
So A will have the form (2) where ( ( and 
7 =(7,,). 

3. Case 3 follows immediately now by first apply- 
ing the permutation in 2 to the rows and columns of 
A and then applying the permutation in 1 to the 
rows and columns containing 7 

THEOREM & Sup pose the blocks A ot the par- 
titioned matria give n in (1) satisfy the equation 


AyX,=X,By, 6 


where B is ad square matria of orde? -. {j ys. 
with strict inequality for at least one value of 7, and 
P is an ti 7 matrir with a nonsingular matrix of 


order fr. Ai’, in the first r rows Let the last 
vi rrows of 4 he As r and let 
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vhere A Sq are } 4 I ¢ j The A N similay 4 
the matria. 3. 
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chere B is a parti oned matria of orde fy vith ble, 
B . as de fined in H). and ¢ has hloel X 


( | LY \ | 
with dimensions ( n } If either yn, 
n r, the corresponding block ¢ does not appear The 
By hypothesis not all n; are equal to r. or else we | sect 
would be left with the matrix Bh B which Would stoc 
be similar to A resu 
Thus the roots of A are the roots of the smaller lf 
matrices B and ¢ rect 
Proor: Let P; be a matrix of order 7 row: 
\ O 
| 
A / : 
then * 


where J, represents an identity matrix of order | 


Since by (7) and (10 


then by a) P df 


and, using (11) and (9 
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i y=P;'A,,P 12) |") 
n p elen 
So, if we let P be the direct sum of the matri 
."« 
P=>-P,=P,+P 
and let 
{—P-'!4P 5 


then A has the blocks A,, given in (12), and 
simultaneous permutation of rows and columns 
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of A given by case 2 of the lemma above, will put it | Th 
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3, Generalizations of Brauer Theorems on 
| Stochastic Matrices 


proved a number of interesting 


A. Brauet 


theorems about veneralized stochastic matrices, 1.e 


matrices A a of order n such that 
S*\a S LB» « n). 14 


Such matrices have as one root, 8, and the vector cor- 
7 < 
responding to this root ts 


U Oe oe (l4a 


The reduction formula of theorem 1 is applied in this 
section to matrices which partitioned into 
stochastic blocks, to give generalizations of Brauer’s 


can be 


results . 

Ifa matrix .1 of order N can be partitioned into 
rectangular (7 7 blocks, in each of which the 
row-sums are all equal, Le., 

A Su), a=1i 2. 8, &=1,2. . 7 
wit! 

>) 4; Bij, 3) ee 15 


Aisa hlock-stochastic Matrix 


THeoreM 2: /f AA is block-stochastic, (1) t roots of A 
are roots of the . and (Z the other N T 
are roots of the partitioned matrix of order N—t. 


matria S . 


If « ther ns or ny 18 . the hlock ( is omitted. SO 
| d ] do not NECESSA] ily take on all pal uve S e 2 . . f 
Proor: Since, by (14a) and (15) 


A,,f Sid an 17 


t ‘ 


statement 1 follows immediately, for by (17), 1f we 


set r=1 in (6) the matrices B,, of theorem 1 are 
simply the elements s,,; given by (15 Then, if we 
subdivide A as in (7), ia » Consists of all the 
elements of the block A excepl for those in the 


first row and column, i.e 


| a p==Z, } e=2, } 
(ii) —(q k=2,...,% 
Since by (17), the matrices X, in (6) are the 
ectors (. . we have 
XV 1, X / 
Then, since [ Aj? is a matrix having in each 


ow the elements 


statement 2 follows from (9 For example, the 
following matrix arose from a problem in molecular 


physics 








[¢ bh b b Cc ¢ cddy 
h a b b Cc ddce 
bbhiaod cded 
bh ba ddede 
A ce \iddieft fff 
cdicd fieffgq 
cd de J fe gt 
decd f faef 

Lde def g fe) 


The lines are drawn in to show how A is subdivided 
into stochastic blocks. Then by theorem 3, the roots 








of A are those of the 5&5 matrix, 
ey bh 0 c—d c—d d—c>) 
) a h ( d d—¢ ( d 
( c—d c—d e—dg 0 0) 
c—d 0) f—q Gun if _ 
. 0 ¢—¢ {—@q q—f] f fJ 


and those of the 4 4 matrix, 








& bh 2b c 29c+2d *) 
2h a+-b d 20 rd 
S . 
2 2d é 4f 
_ ¢ d rf d } t- 2f ; al 


It is easy to see that C can again be divided into 
stochastic blocks where n 2 7 Ll, t%=2. SO 


the five roots of C are those of the 2 2 matrix. 


a—hb 2c—2d 
( 
( d ‘ qd 
and the 33 matrix 
a—b c—d U 
( 2 dd é gd 0) x 
( d | q e—2} q 


Thus e—2/-+-¢ is a root, and there are two degen- 
erate roots since the matrix oF corresponds to the 
quadratic matrix which remains when the 
e—2}/-+q9 1s factored out of C4. 

Coro.uary (THEOREM 2): Jf 
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O—>>-., where P, and Q, are perm utation matrices of 
‘ haul i « 


supplying this example 
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order n;, any matrix PAQ is also bloek-stochastie and 4. Other Applications 


has t roots in common with A. that 
Proor: Since the row-sums in a given block would In this section we will mention several other types 
be unchanged by permutations within the block, this — of partitioned matrices, A A,,), such that , 
follows immediately from statement | of theorem 2 and 
Brauer (see footnote 2) proved the following B ' 
theorems for certain special stochastic matrices PAP ( ) on, whe 
1. Assume that 0 ( ' 
a c gnZ, 3, ...8 9<9). where P is the direct sum of transformation matrices 
P,, and B a ( , are partitioned Matrices 
Then A has the same roots as the diagonal matrix which will be defined in each case. We also giy, 
a case where the transformation PAP produces L 
B=diag (a22—C2, 33—C3, . . -, a Cup 8). a real matrix from a complex one 
Since the proof in each case is the same ay 
2. Assume consists merely in performing the indicated matris 
multiplications, it will not be given in detail. | 
Q;;—C, JF i, and a hb 
4.1. Block-Circulant Matrices 
7 Jl» J / 
This is the case where each block A,, is a circulant 
matrix of order 7 The result (21) is given b 
Williamson 
Then (6—e) is a root of multiplicity 7 The roots of A are roots of the n matrices , 


We generalize these results in the following order f, 
theorem for block-stochastic matrices 


THeoremM 3: Let A be a block-stochastic matrix with yy, is ) h—-1.2 ae 
submatrice 8 A (Ay, having constants ec, above the , j—_ : coe he i It 
diagonal in the lastr columns of each block, 0 / Tih 7 toric 
E+, where «, is one of the nth roots of Unity 

a c k>n, rm h<ok+n n;). IS 
. , 4.2. Block-Singular Matrices 
The ntr of the roots of A are roots of the matrices . ‘- 
If A is singular, of rank n—-1, there exists a veetor| (a 
. r,, Le, , wv), With components not all zero, | pon 
Cex a , (k=n aror ; 19 
Ckk ‘ Jai , , such that 

mn , : AX=X-@. 

(The second of Brauer’s theorems comes as a special since 
case of the above if we permute the rows and columns \ 
so that the columns ) : jo, 5 * are the last } columns Cc i without loss ot ag 1 Ns 

» ; ¢ ; «.* } , Qc hea , , 
of the matrix. Also it isn’t necessary that all the len, as in theorem 1, Uf . 1+A2, ¥ 
nondiagonal elements of a column be equal as in can let ‘ 

Brauer’s hypothesis U 

Proor: By theorem 2, N-f roots of A are roots l -- 
of C=(C;,). Here we have 7, j= 1,2, , t, since by NX, fd 
hypothesis, O< 7 min nm, and thus min » | 

; and, applying the transformation (16), we can show 
\ (16) 3 (1S), ¢ Sa parth ‘angular Nd 

But by (16) and (18), = 6 parteny Geenguer 4..: 1h. nemacre roote of A are roots of | 
matrix with elements a c on the diagonal of 
the triangular portion, i.e., for k=n,—r+1, 7 ; x So 
Then, using case 3 of the lemma above, we can - - ee eo eg ar 
permute the rows and columns of A so that A ts ; 
similar to the matrix If A has rank less than n—1, the same type 0 

transformation can be applied to the matrix C and jhoeg, 

P , this can be continued until a nonsingular matrix 8  tyjgn 

R obtained. Or, we may be able to partition A IMl0 pops, 

; 

0 Q blocks so that elem: 

AyjX,=X; 0. Will 

when (= ((,,) is a matrix of order tr having square — 
blocks of order ¢, defined by (15), on the diagonal, Then if we let P=S3P,, where each P; has th 
and (,,.—0 for kh form (22), and apply transformation (20), we show 
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If. in addition, the matrix (a,,) Is stochastic, with 


that the roots of A are the roots of 
row sum s, by theorem 1, A has as n of its roots the 








Lypes B B ) of order / number S l/n+-1)2t cos her (k l, = éom Oe The other 
al n(t—1) roots are degenerate with multiplicity 1, 
( ( since the matrices (’,, in (5) become, by (22) for any 
- where value ol Mi. 
(’ ad, ra’ & — _ 
a; ay; (i, 7 «oe 27) 
h=2 n k—2 n 
LUT ees 
trices . 4. imi re) 1M 
| ie 4.3. Block Tridiagonal Matrices 4.4. Complex Matrices Similar to Real Matrices 
a Let 
duces a 0 0 a> A A,,) where each block is a 22 matrix defined 
} —_ as follows: 
an 
natrin e _u 0 
ad \ lh () 
() a () () A i=] j 
0) a,— ,—I1b 
T a 23 
‘ulant ‘ ake , . ; . 
m b We use y—1 instead of 7 to avoid confusion with 
the subscript 
| 
eS 0 () () 0...a4 
ad h 
h a 
: It is well known’ that such matrices have charac- 
ristic roots which may be written in the form a b 
A / 1.2 h J h | t) 
\.—a+2 cos (kr/'n+1 k=1,2,....” 24 a b 
nd, if we let @#,—km/n+-1 it is easy to show that the a a 
ect haracteristiec vector corresponding to A, has com- A ( i. | ae 1} k) 
707r0 onents h h 
sin 6, sin 24@,, , sin né 25 
| and .1,, an arbitrary real matrix ((,j;—k-+-1,.. .,f). 
since the trigonometric identities Then. 
u asin 6é sin 26 a-+-2 cos 6) sin @. 
¥ 1+, | l—,—1 
m—1)0 asin mé sin (m 1) @ » ) 
») ; P ; 12 ] 
Zé a - COS 7 SL wg mn be a . 
l—, | I++ | 
sin (7 1) 6 asin né a-+-2 cos @,) sin né 2 2 
Snow . : ° 
old for each @, as defined above P ] i—k+1,...,8 
So if A A where the blocks A are all tri- a 4 
A T ad h a 
pe ay \ 
and hese blocks can be simultaneously reduced Lo PA P { ca | }-: 
mix triangular form using the matrix of characteristic —e rd =_— 
| inl vectors (25), since these are independent of the : 
ements d,;. Thus, according to a theorem by PA A ‘ =r" J ke veo x8 
Williamson, the roots of A are those of the » 
Matrices, AvP At =—k+1,...,t g=il,.. k 
su 
Qi;t-2 cos (4/n-+-1)), = 
show Since 
ad » COS WW hi ] . (26 
} Rutherford SOrTit mtinuant determinants il Tit tl pl Vsics and “" ‘ 
stry, Il, Proc. Roy. Soc. Edinburgh 63, III (1952 A=PAP-" 28) 
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has the same elements as «1 except in the 22 blocks 
on the diagonal for i=1, 2,..., k, and the 
in these blocks are real. (Matrices of this 
arose from an applied problem at the 

Bureau of Standards 


elements 
type also 


National 


5. Applications to Determinants and 
Systems of Equations 


The transformations given here for a partitioned 
matrix A=(A,,;) can obviously be applied to finding 
the determinant of such a matrix since if 


4=paP"=(7 ?), 29 
i | B.\¢ 


These transformations could also be profitably 
applied to the solution of a system of linear equations 
with coefficient matrix A; Le 


tions 


the svstem of equa- 
AX=K 30) 


is, by (29), the same as 


P-\AP)X=K 
or 
AX=—K 
where 
PK=K 
and 
PX=% 


and A=K Ky ar 


Then if Y= A \ 


to correspond With thre 


Partitione 
dimensions of B and ¢ 
can first solve the system 


~ ~ 


CX Ky 


and then the S\ stem 


BX ,=—kK 


So. solving these two systems ol linear equations 
and the comparatively simple system (31) (wh) 
can be partitioned as /? is partitioned 
solution of (29 

Such a transformation would be especially valuab) 
in solving a large system on an electronic computer 


we have ¢} 


where the length of time for solution goes up o \ 
the order of nm. Thus, if we cult the dimensions , 
half by transformation (29), systems (32) and (33 NI 
would each take one-« iwhth 2s long LO solve as the \ 
: yur 
original system, and if we assumed only that / 


could be partitioned Into tWo parts the total solutior 
time would be at most one-half that of the origing 
system. In many cases the time 
approximately one-fourth or less . 


could he Cut tf 


WasHineron, D.C Paper 63B2-9)| car 
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Capacity Requirement of 


a Mail Sorting Device: II 


A. J. Goldman 


July 14, 1959 
The combinatorial analysis of a mathematical model of a sorting device suggested by 
S. Henig is completed. The relevant parameters are the number of destinations for the 
alions mail, and &, the number of letters entering the device during each cycle of operation The 
whiel apacity of the device, if it is never to jam, should be between rk and ri r—1) inclusive 
ve t] irguments indicating that the latter value is preferable are given 
luabk 1. Introduction capacity greater than rk—(r—1) would be only 
puter transient. 
dhe We deal with a highly idealized mathematical The author acknowledges several improvements in 
ie model of a mail sorting device suggested by S. Henig the exposition suggested by J. R. Rosenblatt. Un- 
| (33 NBS Electronic Instrumentation Section . For our fortunately . somewhat more formal treatment than 
the purposes the operation of the device can be deseribed * that ol the orw seemed required by the 
* as follows. Mail (to any of r destinations) enters 
lution he svstem. After & letters have entered, the device 
— “asks itself’ to which destination it contains the 
ul ! et letters (ll letters to this predominant destina- 
are then dropped out of the device, another & 
tters enter, and the process continues The dev ice 
ssaid to jam if, after a dropout, it is still so “full” 
at entrance ol the next & letters would cause an 
verflow 
Investigation of the most appropriate capacity for 
such a device leads to mathematically interesting 
roblems of two general types. If the capacity is 
to be so chosen that jamming never occurs, then the 
1 oblems ure ¢ ssentially combinatorial in nature and 
B2-9)| can be treated mainly by “counting” methods. — If, 


owever, the object is the more modest one of keep- 
g the frequency of jamming down to some specified 
tolerable level,”’ then rather difficult’ probabilistic 
problems arise; this due to the fact * that the 
behavior of the device constitutes a Markov chain 
tha great many states, governed by the probabilis- 
tie distribution of mail by destinations 

The present paper completes an earlier analysis 
f tl which that 


Is 


ie combinator ial quest Ons, show ed 


sa possible capacity for the dey ice, and (if the dey ice 


gins Operation empty) is in fact the minimum 
aPacity This left open the possibility of gaining 
some advantage by using a greater capacity than 


that given in (1) and beginning operation with the 
vice partly or entirely filled. It will be proved 
below that no capacity in excess of rk should be 


ploved, and that advantage arising from a 


11h \ 


greater complexity of the analysis below. 


2. Statement of Results 


In order to describe our results, it is convenient to 


define 
rtf 


number of letters in the device just (2) 


before the tth dropout, 


and also to give a precise definition of “capacity.” 

A non-negative integer C will be called a capacity 
for the device if, should the device begin operation 
containing no more than ( letters, it can never sub- 
sequently contain more than C letters under the 
dropout rule described above; symbolically, 

r(ty))<C implies 2(t)<C forall t>tp. (3 
This condition is clearly necessary to insure that the 
device never jams 

THEOREM 1. Cis a capacity if and only if 


rk—(r—1) <C. 


A capacity ¢ will be called efhicie nt if there is at 
least one set of initial contents for the device, with 
no more than © letters, such that the device might 
wssibly contain C letters again at some later time; 
symbolically, for some set of initial contents, 


rt (’ for f 


r(f <) and max poss to. (4) 


If this condition is not satisfied, then the full capaec- 

itv (' of the device will not actually be needed after 

the start of operation, a wasteful situation. 
THEOREM 2 (is an efhicre nt capacity if and only 


rk—(r—1) <C<rk. 

A CAapac4ily ( will be called ESSE ntially efhicre net * if 
there is at least one set of initial contents with no 
more than © letters such that the device with non- 
zero probability will contain C letters infinitely 


the subscript ‘‘es n wa 
+) rX 
theory 


TT cs the tern ent her and of 
i! Tre irt 


sug- 
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often 


2(t)<C and 2(t) mex es‘ for t >t. 5 


This definition is strengthened 


essentially efficient capacity by requiring that thi 


stated condition hold for a// initial contents of no 


more than ( letters: t.e., that 


r(to <C for f lo. 6 


implies 2(t)max ess = 
" THEOREM 3. 
ity is rk (7 ] C: 


efficie nt. 


The proofs of these theorems are given in section 


3. Section 4 
estimate related 


contains a very crude probabilistic 


to theorem 3. 
3. Proofs 
3.0. Preliminaries 


It is helpful to define 


Mit max s,(¢ 
where 
r,(t)=number of letters to the 7th destina- 
tion in the device just before the fth 
dropout. 
LemMa 1. M(t) >k im plies r(t+] rt 
LeMMA 2. x(t rk implies r(t+] rt 
LemMMaA 3. M(t I implies r(t+- 1 r(ft 
Lemma 4. M(t)>k implies x(t+1)<2x(t 
Lemma 5. x(t) >rk—(r—1) implies x(t+1) <2x(t 
LEMMA 6. as(to) <rk—( | implies 
rf rk } | for all , ae, 
LemMaA 7. s(t.) >rk—(r—1) implies 
rit)<z for all t>f 


Lemma | is proved by observing that if \/(¢ ke. 
then more than & letters leave the device in the 
tth dropout, whereas only & letters enter before the 
(t-+-1)-st dropout. Since there are only r destina- 
tions, z(t) >rk implies V(t) >k, and lemma 2 
follows from lemma 1. Lemma 3 is an obvious con- 
sequence of the dropout rule, and lemma 4 is obtained 
by combining lemmas 1 and 3 

To prove lemma 5, note that since there are only 
r destinations, 


st) 


x(t) >rk—(r—1 rik—1) implies M(t k—1, 


so that lemma 4 can be applied. Lemma 6 merely 


states the fact that rk (7 | IS a Capacity. In 
7 If this condition is not f hen w 
from some point on, fail to use D 


7; symbolically, for some set of initial contents, 


proving lemma 7, two possibilities must pe 


to that of a uniformly 


The only OSNE ntially é fh rent ca pac- 
itis also uniformly CSSE ntrially 


. » COL 
sidered. If, on the one hand, 
r(t)>rk r—] forall t>t. 
» then by lemma 5, 2(f)<.2(f,) for all t>f If, on ¢ 
other hand, there is a / f, such that 
rif > rh | for ee } 
but 
r(f rh / | 
then by lemma 5, r(t)<r@ for ty <t<t,, andh 
lemma 6 
r(f rk } l)<r(t forall ¢>t#,. 
so that again 
rift rtf ior all f>T 
3.1. Proof of Theorem 1 
To prove the “only if’ statement, SUppos 
( rk } 1) and v(t) <C. If, on the one han 
r(t, 0, then (see op. cit. footnote 5, p. 79) there 
a sequence of possible events leading to a situation jy 
which 
r(f rh ; | ( t<f ; 1) (J | 


This violates condition and so (’is not a capacity 
If, on the other hand, x(f 0, then a 
possible events can constructed 
situation in which 


sequence 


be leading to 


To do this. choose the / letters entering the devi 


between the tth and ft 1 )-st dropouts, 

ty <t<t,, to consist entirely of letters to some or 
destination for which 2,(¢ 0. Then M(t) >k f 
fy <t<t,;, and so by lemma 1 the number of letters i 


the device steadily decreases until condition (7 
satisfied. From this point we may argue as in tl 


ease rt 0) (now using 7,, instead of f again ol 
taining the conclusion that Cis not a capacity 
To prove the “if” statement, suppose 
f rh | and r(f ( 
If «(ft ri } l then by lemma 6 


forall ¢>f 


while if x(t.) >rh r—1), then by lemma 7 


forall t>f 


in either case condition 


capacity. 
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s satisfied, and so Csi 


Th 
des! 


DV 
lett 


at 


3.2. Proof of Theorem 2 


eC COr 
To prove the “if”? statement, suppose 
rk—(r—1) <C<rk 
On t! : 
The initial contents of the device can clearly be 
chosen so that x(t, Cand M(t))=k. Among the 
destinations with z,(¢ k, select the one “preferred” 
by the tie-breaking rule at the f-th dropout; let & 
letters to this destination enter the device between 
the f-th and (f)+1)-st dropouts. Then M(tp+1) 
and by —k: among the destinations with 2,(f)+ 1)—&, select 
the one preferred by the tie-breaking rule at the 
] )-st dropout, and let k& letters to this desti- 
nation enter the device between the t,+-1)-st and 
4+2)-nd dropouts, ete. This construction results 
in M(t)=k for all t>t, so that, by lemma 3, 
r(t rf, ( forall t>fp. 
Hence Cis an efficient capacity. 
‘he The “only if” statement will be proved next. By 
. _ theorem 1, no C<(rk—(r—1) can be a capacity; thus 
rm no such ( can be an efficient capacity. Suppose 
lon Mn} rk and x(t.) <C. There are three possibilities to 
be considered. If 2(to) <rk—(r—1), then by lemma 6 
. x(t) <rk ; 1) <rk<( forall t>t, 
acity . 
a so that Cis not an efficient capacity. If 
ce 
v0 rk—(s 1)<2x(t,) <rk, 
.| then by lemma 7 
r(t) <a(t,) <rk ( for all t> Io, 
qevir 


30 that C' is not an efficient capacity. Finally, if 


Lf ) >rk, then by lemma 2 there is a ft, >t) (with 
* t<f ( rk) such that 
ters 
rh r(ft r(f ( for f t<f 
in Uf , 
ind 
In OD- 
r(f,) rk: 
} 
the arguments for the first two possibilities can be 
applied to tf, to vield 
r(t) <ax(t,) <rkh ( forall ¢>f? 
sO that z(t (' for all t>t). and again (C' is not an 
elicient capacity 
3.3. Proof of Theorem 3 
) First, rk—(r—1) 7s a uniformly essentially efficient 
- apacity. To see this, note that the construction 


ised above to prove the “only if” part of theorem 1 
shows that if 


then with probability one * there is a 4, >t) for which 


a(t,)=rk—(r—1): 


i 


by the very same argument, with probability one 
there is a ¢, >t, for which 


a(t.) =rk—(r—1 

and so on. In fact, for rk—(r—1) the phrase “‘with 
nonzero probability” in the definition of ‘uniformly 
essentially efficient capacity” could actually have 
been strengthened to ‘‘with probability one.” 

Second, no C>rk—(r—1) 1s an essentially efficient 
capacity. ‘To prove this, observe that by theorem 2 
attention can be confined to the case 
C'<rk. 
Suppose s(t) <C. If 

r(t,) <rk—(r—1) forsome 1t,>to, 

then by lemma 6 


a(t) <rk—(r—1)<cC for all t>h, 


so that the device cannot contain C letters infinitely 


often. Thus we may assume that 
a(t) >rk—(r—1 forall t>t, (8) 
so that, by lemma 5, 
r(t+1) <2(t) forall t>tp. (9 
If x(t (’, then by (9 
r(t) <2r(by C’ forall t>bt, 
so that the device can never contain (letters. 
Therefore assume that 2(f),>=C. If 
M(t,) >k forsome ft, >to, 


then by lemma 1 and (9 


r(t) <<ar(t, +1)<xr(t,) <2(ty) =C for all t>ht, 


so that the device cannot contain C’ letters infinitely 


often. ‘Therefore assume 
W(t) <k forall t>t. 
From (8) and the fact that there are only 7 destina 


tions it follows that .\/(t)>k, and so M(t)=k. Thus 
the only possibility for the device to contain C 
letters infinitely often is given by 


r(t OA Mit I: forall t>t: 


this clearly has probability zero. 


‘TT onstructio hows that if fo)<ri r 1), then fo inv block olf 
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4. A Probabilistic Estimate 


Suppose C’ is an efficient capacity which is not 


essentially efficient; i.e. (according to theorems 2 
and 3), 
rk—() L< € rk 
We know that with probability one, there is a 
t, >t) such that 
arlt)< 0 forall t >t 


It is of interest to estimate how early the first such ¢ 
(the last moment at which the full capacity of the 
device used) might occur. In this primarily 
combinatorial paper, we will be content with an 
exceedingly crude probabilistic estimate 

Denote by FE the event 


Is 


r(to Mit, k 


and 
By the proof of theorem 3, ¢;=t) unless / occurs 
Therefore the function of ¢; to be estimated is 

P(t,)= Prob: x(t QO for some ft >t, # 10 
where Prob: A B denotes the conditional probability 


of A given B. 


Let the relative frequencies (or probabilities of 
letters to the different destinations be denoted by 
ox §, <I / 
let 
( rik ] \ | () \ } 
and let 
Pe=fits f 
It will first be shown that, for any t >t, 
Probiz(t)=C, Mi(t)=kaelt—1I=C, Mit—1 i 
| I 1] 


Suppose 
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since 
r(f | ( r(f / } | sk. 

it follows that 2,(t—1)=f for at least s+1 desting. 
tions The letters Lo some one of these destinations A 
leave the device in the (f—1)-th dropout, and 
W(t) >k if any of the & letters which next enter th, 
device vO to any of the remaming destinations 
with & letters each Since the sum of any s fre. 
quencies f, is at least I the inequality ll holds 

By the proof of theorem 3, for any one ¢ 


both 2(t—1 (’ and MW(t—1 k hold if and only 
r(t’ ( ana V(t’ } or t££.4F<43 
Thus the probability on the left side of (11) can by 
rewritten 
Prob) «(ft CO. Mit kor (t’ C. Mit’ ke foy 
<2 
so that the inequality (11) vields 79 
n 
IT Prob rift $2 VJ Port ( \V/ t’ ht, / ols 
ol 
oO? r ; ; | fk sor 
| 
This last product is equal to _ 
hol 
Probiz(t{)=C, Mith=k for th<t<tiE = 
WSS 
which, according to the proof of theorem 3, is a _ 
uppel bound for Pat and so we have proved that , 
Pit |—f 1? r 
If for example k=2 and / 0.001, then after 500) #PI 
dropouts the probability that the device will make | ™" 
any further use of its full capacity is less than (0.50, mT 
ai 
10! 
WaAsHINGTON, D.C Paper 63B2-10 " 
of 
ui 
we 
Ol 
\ 
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ations 
and 
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Analytic methods are 


niques are applied to a specific numerical problem; sé 


are quickly eliminated, and a “hybrid 
comes Within four percent of optimum (if 


1. Introduction 


A number of plausible configurations for auto- 
matic mail sorting equipment have been suggested. 
In this report methods for the comparative study 
of such configurations are developed, with the object 
the configuration which 
sorting of mail at a given (required) rate at least cost. 

The solution to this ‘‘optimization” 
course depends upon the distribution (by destina- 
tions) and volume of mail at the sorting installation, 
and also on the operating parameters and costs 
associated with the various components of the sorting 
system. Therefore, in order to present our methods 
in the most comprehensible way (i.e., in the context 
of a specific numerical application , It Was necessary 
to work with a definite set of realistic “input” data. 
These inputs are described in sections 3 to 6 and 
appendix A We wish to emphasize that we work 
with these particular data only to illustrate the 
methods empl ved; the reader who follows out the 
ealculations for this specific case should have little 
trouble carrving out the corresponding calculations 
for other data 

The results of the analyses are summarized in 
section 2. In brief, we are able (7) to rule out several 


achieves 


of the proposed configurations, (77) to show that 
one of the proposed configurations will (if a certain 
parameter is choosen correctly) come within 8 per 


ent of Optimum, and (777) to exhibit a “hybrid” of 
two of the proposed configurations which (if a 
certain pair of parameters is chosen correctly) will 
come within 4 percent of optimum 

We also determine how accurately the parameters 
m (7) and (777) must In view i 
analysis ol more complicated configurations did not 
seem worthwhile 

For a comparative study, it is necessary only to 
onsider which vary appreciably from one 
configuration to another. Our analysis is set up so 
that these variable costs are all absorbed into costs 
the devices which inject mail 
we will generally use the 


bye chosen. of ‘di 


Costs 


associated with a 


nto the sorting system 
Office Department 


Re 


“au cerlall 


problem of 


5. | 959 


ce veloped to aid in determining the equipment configuration 
which achieves sorting of outgoing mail at a given 


The tech- 
suggested configurations 


required) rate at least cost 


veral of the 


of two of the proposed configurations is found which 


pair of parameters is chosen correctly 


term /oading complex for such a device, instead of the 
longer “distributor loader complex’’), (b) the 
ceptacles (6¢ns) in which the mail ends after passing 
through each stage of the system, and (c) the opera- 
tion (sweeping) of removing the mail from these bins. 

A natural method of determining the cost of any 
proposed configuration, therefore, is to determine 
the numbers of loading complexes, bins, and sweepers 
involved, and then to multiply each of these numbers 
by the appropriate cost coefficient and add up the 
resulting partial costs. We might call this the 
add up method; it is described in section 7 and at the 
end of section 8. We have also developed a follow 
through method based on the idea of following an 
individual letter through the sorting system, charg- 
ing it an appropriate fraction of the cost of each 
loading complex through which it passes, each load- 
ing complex by which it passes (and thus prevents 
from operating at maximum rate), each bin it 
cupies, and each sweep it receives. The follow 
through method is not strictly applicable to some 
configurations (though it seem to provide 
initial approximations which may shorten some of 
the work of the ‘add up” method); its main ad- 
vantage is that it permits a much closer estimation? 
of the theoretical minimum cost described in the next 


OC- 


does 


paragraph. 

An important step in the anglysis is the derivation 
of a theoretical minimum (more precisely, a 
“lower bound’’) which is independent of the configu- 
ration. This cost is minimum in the sense that “‘any 
actual configuration must cost at least this much” 
and is theoretical in the sense that “no actual con- 
figuration can cost quite this little’; it therefore pro- 
vides a yardstick against which the costs of specific 


cost 


configurations can be measured, the deviation from 
optimality of a configuration can be assessed, and 
the permissibility of plausible approximations can 
be checked. A rough estimate of such a minimum 
cost is given in section 8 (without using the specific 
mail distribution): in section 9 the “follow through” 
method is emploved to derive a significantly more 
result 


accurate 
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The detailed analyses of the configurations studied 
are given in sections 10 to 14 and appendixes B and 
C. We hope that our rather full exposition of this 
material will prove helpful to the reader faced with 
the problem of carrying out similar analyses. 


2. Summary of Results 


From our specific numerical results (which are valid 
only for the particular data and assumptions used 
in our analysis) we can infer certain qualitative results, 
believed to be valid for post-offices (7) for which 
the sorting system must be able to handle about 
1,000,000 letters/hr or more and (7/7) in which local 
mail is removed before the outgoing sort. 

Of these two types of results, the qualitative ones 
are more informative, and so we list them first: 

(a) The simplex scheme and Christmas Tree 
scheme * can be eliminated from further considera- 
tion. 

(6) The residue scheme (if properly 
quite good; some multiple input schemes may pos- 
sibly be competitive. 

(c) A promising approach is to combine the basic 
ideas of the residue and multiple input schemes; 
the resulting “hybrid”’ scheme has lower than 
either of the two original schemes if the relevant 
parameters are correctly chosen. 

We turn now to the specific numerical results. 
The “inputs” leading to these results are given in 
detail in sections 3-6 and appendix A; at this point 
we only add to (7) and (77) above the facts that 
(i117) only costs which vary between different equip- 
ment configurations were considered, (iv) the use of 
automatic transfer equipment was not considered, 
and (v) the possibility of memory-sharing between 
oading complexes was not considered. 


chosen Is 


cost 


Condition (777) is a natural one, since we are 
primarily interested in comparing different configu- 
rations. Auxiliary calculations (not given in this 


that the qualitative results reported 


paper) show 
and 


above do not depend on (iv 

As noted in section 1, a “theoretical minimum 
used as a vardstick which the de- 
sirability of any proposed configuration can be meas- 
ured. Kor this purpose we use the ratio 


cost’’ is against 


The Or 


Theoretical minimum 


Cost of proposed config tical minimu 


“Near-optimal’’ configurations are those with sma// 
values of R. In the notation of the body of the re- 
port, (2.1) becomes 

R C'— Cn) /Coin- 2.2 


as follows 
(\ early 


Our specific numerical results are 
(a) The theoretical minimum 
cost is about $2,050,000. 
(6) For the simplex scheme, 


variable 


R is about 160 percent; 


These schemes are described they ‘ ter the report 


Le. the scheme costs about 2.6 times the theoretieg 
minimum cost. 
(c) For the Christmas Tree s heme, R is at leas 


69 perce nt: 1.e., the scheme costs at least 
the theoretical minimum cost 

(7) For the residue 
peree nt 


There 


1.69 times 


best scheme, R IS about Q 
are a many ways in which th, 
ideas of the residue and multiple input schemes 
can be combined. An analvsis of the full range of 
possibilities would be bevond the scope of our study 
We have, however, examined a relatively simpl 
svstems of this tvpe, in which each sub. 
involves only two loading complexes jj 
optimal configuration within this elass 
has R about 4 percent, which (4) compares favorab) 
with all other configurations studied, and (jj) jp. 
dicates that investigation of more complicated 
systems would not be worthwhile 

(f) All multiple input with more thay 
four subsystems (i.e., more than four series of loading 
complexes) have R at least 8 percent (but probably 
substantially more); multiple input schemes 
four subsystems have R at least 4 
probably substantially more). We 
complete proof that multiple input schemes wit} 
fewer than four subsystems have excessively hig! 
values of R (in comparison with (d) and (e)), but 
strongly believe that this is the Each sub- 
svstem of such a configuration would involve nine o; 
more loading complexes in series, 


creat 
basic 


class of 
system 
series: the 


schemes 


wit! 
percent but 
do not have a 


case 


3. Equipment Cost Data 


Only those equipment costs which appeared likelh 


to vary appreciably from one configuration 

another were considered The cost of coding the 
letters, for example, was eglected on the frounds 
that in a code sort system, all letters to be sorte 
must be coded once and only once, regardless of 


undergo 
whose 


how many readings or sorts they 

The pieces of equipment 
considered variable are the loading complexes al 
the modules of bins Hereafter, the 
the sorted mail will be 


costs of the corresponding 


two costs wel 
cost of bins 10! 
include 

cart pockets, et 


considered to 
COnVeEVOr, 


The following tentative cost estimates were Givi 
bv a representative of the manufacturer 
$2? 000 module of 30 bins. 3.] 
$1?5.060 first load ig complex of a sorter, 32 
$50.000/each of the next 4 loading 
complexes of a sorter, 3.3 
SLO.000 pall of end preces 3.4 


Jecause ad 


total equipment costs are somewhat 
lower than total personne! costs (see the three sum- 
mands of (8.6), for example) and (b) the estimates 

(3.1) to (3.4) are rough it seemed per- 


missible to be rather loose in our treatment of equip- 


only Ones, 


[sin 


Forn 
naly 


a | ical 


least 
Imes 


ut & 
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emes 
ze Oo} 
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sub- 
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in- 
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thar 
ding 
ably 
Wil! 
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hig! 
but 
sub- 
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sum- 
lates 
per- 
| UIp- 


ment costs. First, the end-piece costs are small 
lative to other costs, and are difficult to handle 
because of the way in which the number of end- 
pieces depends on the particular placement of 
equipment For these reasons, the end-piece costs 
have been neglected altogether. Second, there was 
apparently some doubt as to whether the extra 
loading complexes of a multi-input sorter can actually 
bp used in the manner assumed in the estimate 
»2)- for this reason we have disregarded (3.3) and 
taken the cost of every loading complex to be 
625.000. Our numerical work for specific sorting 
schemes is fairly sensitive to these decisions, but 
the results of comparing different schemes are not. 
In summary, the equipment costs actually used in 
what follows are 


») 
6 


$? OOO module of 30 bins (; 
$125,000 loading complex 


These are initial costs; yearly costs are based on a 


ij-y7 am whi ation period 


4. Space Cost Data 
The only space requirements taken into account # 
were those for the loading complexes themselves, for 
pins and for the sweepers work area around them: 


8.75 ft=length of a loading complex (4.1 
875 ft=—leneth of a module of 30 bins (4.2 
25 ft=width of a module, including 
sweepers’ work space (4.3 
[sing the tentative cost firure 
$2.20/ft ?/vr (4.4 
suggested by L. Allison (NBS Electronic Instru- 
mentation Section) we obtain (8.75 11.25 2.2 
22, or 
$220 loading complex/yvr 1.5 
$220 module of 30 bins/yr. (4.6 


It is convenient to combine space and equipment 
sts. From 3.6) and $.5 remembering the 
l)-vr amortization applying to (3.6), we have 


$12,500 + $220 ~$12,700/loading complex/yr. (4.7 
Similarly, from (3.5) and (4.6) we have $200+- $220 
$420 module of 30 bins/yr, so that we have 
$14/bin/vr. 1.8 


formulas (4.7) and (4.8) are the ones used in our 


halysis. 
5. Personnel Cost Data 


The personnel which appears to vary appreciably 
om one configuration to another is composed of 


ee 


pi ed likely t vary appreciabl 


Iration t nother 


sweepers and pouchers. The tentative data used 
below were suggested by L. Allison. 

Sweep operators remove sorted letters from the 
bins and either tie them out or place them in trays 
preceding a further sort. We take 

60 letters=number of letters a sweeper can be 
expected to handle in one sweep, 


(5.1) 
180 sweeps/hr=working rate expected of a 
sweeper. (5.2) 


It is assumed that the sweeping is carried on in such 
away that: 


Each bin is swept at least once per hour, (5.3) 


if more than 60 letters/hr are expected in a 
particular bin, then (subject to over-rule by (5.3)) 
the bin is not swept when it contains fewer than 
60 letters, and (5.4) 


no more than 20 sweeps (i.e., 1,200 letters) /hr 


is expected for any one bin. (5.5) 


The rules (5.1) to (5.5) permit computation of the 
number of sweepers required, if it is known how 
many letters/hr are expected in each bin. This 
quantity can in turn be computed (as will be done 
in later sections) for any particular sorting scheme 
and mail distribution. 

Pouching operators toss tied bundles of mail into 
pouches; their working rate is 


540 sweeps pouched hr. 5.6) 


Comparing (5.2) and (5.6), we see that, for sorts 


which are immediately followed by pouching, 


number of pouchers=1/3% (number of sweepers). 


(5.7) 


As for sorts which are not followed by pouching (i.e., 
primary sorts to be followed by secondary sorts), 
we Observe that the additional work involved in 
handling the trays and feeding the loading com- 
plexes for a secondary sort appears to be of the same 
order of magnitude as that of pouching a similar 
quantity of mail. Therefore, since we suppose the 
transfer between nonautomatic, we 
will take the cost of the personnel needed for this 
transfer to be equal to that of a number of fictitious 
pouchers riven by (5.7 
The average salary figure used is 


sorts to be 


$11,000/man-position/yr, (5.8) 
where a ‘“‘man-position”’ requires more than one per- 
son because of the several shifts worked and the 7-day 
week involved 


6. Mail Distribution 


The particular mail distribution assumed in our 
numerical work is a hypothetical one obtained by 


required for initial loading on the primary sort is 


Ir ntrast. the personnel 
] t f ition, and so is not considered 


equipment configur 
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modifying data of outgoing mail from Los Angeles.' 
It was decided that for a code sort machine, at least, 
all local, postage due, uncancelled mail, ete 
removed at the coding station. 
go backs and residue (about 4%) was excluded from 
sorter input because of the difficulty of distributing it 
properly by destination. The mail distribution to 
the sorting system, therefore, consists approximately 
of what is left after these deletions, with the percent- 
age of mail to each destination upgraded to bring the 
total to 100 percent, and the amount of mail to each 
destination upgraded to bring the total to 


1,000,000 letters/hr. 6.1 
Our distribution involves 
1.600 destinations. (6.2 


The details of the distribution are given in appendix 
A; they involve some inconsequential grouping in the 
“tail” of the distribution. 


7. Cost Formula 
The three variables in the cost formula are 
L=number of loading complexes, 


B=number of bins, and 
S=number of sweeper man-positions. 


The values of these variables can be found for any 
particular mail distribution and specified arrange- 
ment of equipment. 


The total yearly cost C is given by 


C=C,4+C,+¢ 

where 

(’.=yvearly equipment cost 

C. vearly space cost, and 

C, vearly personne! cost 
Using (4.7) and (4.8), we have 

Cc, ( 12.700 L 4B 
using (5.7) and (5.8) we have 

C’,=11,000 (S+3S 14.670 S 


The cost formula we will use is therefor 


C= 12,700 L+-14 B+-14,670 S 7.4 


where L, B, and S are defined above and 


O= variable dollar cost/vr. 7.5 
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8. General Minimum Cost Estimate 


We will first derive a theoretical MINIMUM cosy 
estimate which IS dele ral in the sense that it does 
depend on the particular distribution by desting. 
tions of mail, but only on the volume of mail: i.e. 
the fact that the system must 
handle 


» OD 


sorting be able to 


1,000,000 letters/hr. 8 ] 


i 


Using this and the fact that the maXimMum Input rate 
Is 


36,000 letters/loading complex/hr, 


we that least 1.000.000 
complexes are required; since L 


36,000 
isan integer, 


see at a7 


loading 


L>2s8 83 


| 200 letters 
the number of bins re- 


By 5.5) each bin accounts for at most 
hr, so that in view of (8.1 
quired is at least 


1.000 00071. 200 —S833 


since bins come in 30-bin modules, B must be 4! 
multiple of 30, and so 

B> S840. R 4 
By (5.1) each stack of letters contains at most 6 
letters. so there must. be at /east 

1.000 000 60 stacks ‘hy 
and thus in view of (5.2) there must be at /east 
L. OOO OOO 60), 1S0 Qt sweepers; 

thus S > 95 85 


We now apply en S.4 
taining 


> (12 760 IS 14 N4() 14.670 93 
or, rounding off, 
('>356.000 12.000 564.000 1.732.000. (86 


Thus we have proved that $1,732,000 is an (approx 
able yearly total cost to) equip 


mate) minimum ra 


ment, space, and personne! 


The three summands 


to costs associated with 


in (8.6) refer, respectively, 
loading complexes, bins, and personne! We enm- 
phasize that this is a theoretical minimum cost; DP? 


actual system can cost this little 

Next we will be a little more realistic (and a littl 
less ceneral and use the following one fact about 
the particular mail distribution with which we deal 
Kach of the 700 least frequent destinations receiv’ 


60 or fewer letters/hr. On the one hand, they re 


hot : 
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aes 
tot 


lett 
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l COst 


PS not : 


stina- 
C., On 
le to 


8.1] 
trate 
(8.9 


ading 


tters 


is Te- 


we a 


8.6) 


rox 
(julp- 
ands 
with 
em- 
t: no 


littl 
hout 
deal 
celyt 
y re 


quire at least 700 sweeps hr and therefore vield aut 
least 700 stacks/hr. On the other hand, these 700 
destinations account for about 1.8 percent of the 
total mail, and thus for 18,000 letters hr: if these 
letters were arranged into stacks of size 60 (as was 
assumed in deriving (8.5)), then 


18.000/60—300 stacks/hr. 


eather than 700, would result. Thus the argument 


loading to (8.5) underestimates the minimum pos- 
jeu Ing , & 
sible number of stacks/hr by 700—300=400, and 
thus (using (5.2 underestimates S by 400/180 ~2 
sweepers. 8.5) should be replaced by 
S>95, (8.7 
and the corresponding modification of (8.6), after 
rounding, 1s 
('> 1,761,000, 
«yo that we have a minimum cost estimate of 
C_...=$1 461 000. S.S 


The estimate (8.8) 1s still too “general” to use as 
a vardstick, and we shall use instead a “detailed 
minimum estimate” on the detailed 
properties Of the specific mail distribution) which 
will be derived in section 9. This care in choosing 
a vardstick might seen unnecessary, since a “‘yard- 
stick” is only a unit of measurement whose choice 
cannot affect which of two proposed sorting systems 
appears less costly. The choice of vardstick does, 
however, affect our decisions as to what constitutes 
a siqnifeant difference in cost (either between two 
systems being compared or between a single system 
and a hypothetical “minimum-cost’” system), and 
what constitutes an allowable error in 
making simplifying approximations 

We shall use (8.8) primarily as an aid in caleulat- 
ne the the various systems studied If, 
sing (8.3 _and (8.7), we define 


cost (based 


to 


] 
wisO as 


of 


S.4 


costs 


2s 
S40) Ss 


x 9 
LO 
X11] 


AL=L 
ABB 
AS=& 
enit turns out that a convenient wav to calculate 
he cost ( of a system to find these numbers 
AL, AB, and AS of ‘extra’ loading complexes, bins, 


Is / 


and sweepers, (77) to calculate an “extra cost,’ 
AC=C-C 
by .A¢ 12.700) AL+14 AB 14.670) AS (8.12 
this formula follows from (7.4 and (ii) to find 
{ by 
( ( A ¢ 8.13 


This constitutes the “add up”? method mentioned in 
section | 
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9. Detailed Minimum Cost Estimate 


Our main goal in this section is to derive a mini- 
mum cost estimate which makes detailed use of the 
particular mail distribution we are studying. The 
methods used in this derivation will also turn out 
to be helpful in the analysis of some of the systems 
considered later in the report. Reading the first 
paragraph of section 12 may be helpful here. 

All the systems studied later have the property 
that: 


For each destination, either all mail to that 
destination is sorted by destination on the pri- 
mary, or all mail to that destination goes into 
residue on the primary and is then given a sec- 
ondary sort. (9,1) 


We shall therefore assume this property in (mentally) 
constructing a hypothetical “‘minimum-cost”’ 
tem. Since this is to be a minimum-cost system, 
we can also suppose (to minimize the cost of residue 
bins) that: 


sys- 


All residue bins are operating at their maximum 


capacity (see (5.5)) of 1,200 letters/hr. (9.2) 


Once this assumption is made, it follows (in order to 
minimize loading complex costs) that: 
A letter which to into residue will be 
dropped into a residue bin before passing another 
loading complex. (9.3) 


. r 
Is go 


The method of estimating costs used in section 8 
depended on counting the numbers of loading com- 
plexes, bins, and sweepers, and adding up the result- 
ant The approach used below is fundamen- 
tally different, in that it involves following every 
letter through the system and assigning an appro- 
priate cost at each stage of its progress. We first 
consider a letter which enters the system through 
some loading complex %; and then (without being 
dropped) passes a succeeding loading complex ». 
This letter prevents 7, from working at its maximum 


costs. 


capacity of 36,000 letters/hr; some other letter, 
which could otherwise have entered the system 
through 4 will now have to enter the system 


through some ‘“‘extra’’ loading complex (in addition 
to the minimal 28 complexes known to be required 
by (8.3 If we suppose (to minimize the number 
ol extra complexes and thus the total cost of all 
extra complexes) that all extra complexes are oper 
ating at capacity, then our original letter, by passing 
/, has created a requirement of “1/(36,000) of an 
extra loading complex.’’ Since (by 7.4) each extra 
loading complex adds $12,700 to the system, we are 
led to the following rule for use in estimating the 
cost of our hypothetical system: 

Assign a cost of 12,700/36,000 = .35 


3 whenever 
a letter passes a loading complex. (9.4) 


Suppose now that it is possible to divide each 
extra loading complex into 36,000 parts, each able 
to handle 1 letter/hr, and that it is possible to add 
enough of these parts to each of the minimal 28 
‘“nonextra’’ loading complexes required for the pri- 
mary sort to ensure that these complexes operate at 








their full 36,000 letters/hr capacity. Such a policy 
would tend to decrease bin indirectly, 
sweeper costs (since bins which previously received 
fewer than 60 letters/hr might now receive 60 or 
more ; (5.3)). We therefore assume, for our 
hypothetical minimum-cost system, that this policy 
has been adopted, so that the primary sort involves 
28 loading compleres (each possibly augmented by 
parts of extra loading complexes) all working at their 
maximum capacity. If we also make the ‘“‘minimum- 
assumption that all loading complexes used 
are Operating at maximum 
have a situation in which: 


costs and, 


see 


cost’ 
for the secondary 
capacity, then we 


sort 


All loading complexes operate at their maximum 
rate of 36,000 letters/hr. (9.5 


Thus any letter entering the primary or secondary 


sort uses 1/(36.000) of the services of the loading 
complex through which it enters, so that we have 
the following rule analogous to (9.4 
Assign a cost of .353 whenever a letter enters 
the primary or secondary sort 9.6 
The essential assumptions made so far are (9.1 


(9.2), (9.3 and (9.5 Before proceeding further, 
we point out that two of the (physically realizable 

types of sorting configurations to be analyzed late 
actually do satisfy these assumptions. ‘The simplex 
scheme (treated in 10) satifies them exactly, 
obeying vacuously the conditions referring to residue 

the optimal residue scheme found in 12 
satisfies them very nearly. Thus these schemes can 
be treated by the method developed below.’ There 
is no need (see sec. 10) for such a detailed treatment 
of the simplex scheme. We shall treat the residue 
scheme, however; in fact, since the methods of this 
section are less obviously correct than is the “count 
the sweepers, bins, and loading complexes” approach, 

we will later analyze the residue scheme using each 
approach separately and verify that the same answer 
is obtained in both cases. 


sec. 


section 


Returning to the analysis, we define 
f,—fraction of the mail which goes to the 7th 
destination 9.7 
and consider whether or not letters to the 7th des- 
tination should go into residue on the primary sort. 
For this purpose it is convenient to define 
\7 1 000 ,000F if 1.000.0007f, < 60. Q.8 
60 otherwise, 
and to note that (by (6.1 
1.000.000 f,—number of letters/hr to the 7th 
destination 9.9 


In addition, from (7.4) and (5.2) we are led to the 
rule: 
het tual sches mes involving multiple input fail to obe 
sche s by the methods developed below t tak 
ane tes of the loading complexe 


Assign a (personnel) cost of 14,670/180 8150 
for each stack to be swept. (9.10 
If letters to the 7th destination vO into residue 
then the total cost to be assigned to each such lette; 
as it travels through the sorting system can }, 
calculated in the following way. The letter enters 
the primary and secondary sorts, so (9.6) vields » | 
loading complex cost of 2 (.353). By (9.3), ther 
is no contribution from (9.4 By (9.2) and (74 
the letter should “charged” 14/1,200 for its ys 
of a residue bin. By (5.5) and (7.4), assuming 
10°f,< 1,200, it should be charged 14/10°f, for its 
use of a bin in the secondary sort. By (5.1) and 
(9.10) it should be charged 81.50/60 for the sweep 
of its residue bin, while by (5.1), (5.3) and (9.10) ; 


50/V, for the sweep of its bis 
Adding together these partial 


should be charged 81 
in the secondary sort 
costs and defining 


cost to be assigned to a letter to the ith 
destination, if such letters go into residue. 
9.1] 


we have 


S1.50/V 2? 076 
1.200). g 19 


14/10" 


if 10°/ 
Next want to calculate the analogous cost 
letters ith destination do not VO Into residus 
for a secondary The primary of ow 
hypothetical System now of possibly 
augmented) loading complexes in series, each followe 
by a row of bins (some of which may be residue bins 
this series arrangement involves no loss in generality 
since any other arrangement can be 
special ease. We first that 


we 
to the 
sort sort 


consists 28 


obtained as 


note 


bins for the 7th destination are spaced uniformly 
through the primary (except possibly at the end 
9.1 


‘To see why this is so, suppose for example that it is 
found that the first bins for the 7th destination shouk 
be placed after the third loading complex Then 

view of (9.5), the situation regarding this destina- 
tion is the same beginning with the fourth loading 
complex as it was beginning with the first complex 


and for the same reasons as before we would fin 
that the next bins for the destination should be 
placed far after the fourth complex as the first 


bins were placed after the first complex Thus Wt 


can define 


n number of loading 
cessive appearances in the primary of bins 


complexes between su 


for the 7th destination 9.14 
In view of (9.5), a letter to the 7th destination i 
equally likely Lo pass oS. 2. a i l complexes 


each with a oe l/n,, so that the average c0s 


contribution due to (9.4 


Is 
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eff 
det 


In 


COl 


$1.50 353) X(O+1+2+...+(n 1) )/n,=0.1765 (n 1), ive., for n,<0.1231/yf,, 
(9.10 ; 
; where the formula for the sum of an arithmetic and is increasing for higher values of n;. Since 
Sidue. progression has been used. Chere Is a contribution n ~ se. the minimizing value of Ny is 
letter of 353 due to (9.6). The system contains 28/n 
aN be! cots of bins for the ith destination, each set (dis- 0.1231/y7 if this is <28. 
enters pegarding fract ional effects) receiving ; 
ldani 28 otherwise, 
a 

ther 1.000.000 7,/(28/n,) letters/hr: 
(7.4 1.€., 8 
is Use if we define : ' 
iming eae \/ if *£,>4.439 10 
~ , 1,000,000 f,/(28/n if this is <1,200, L 28 otherwise. (9.21) 
and 1.200 otherwise, (9.15 
pr By (9.20) and (9.19 
od the letter should be charged (again ignoring fractional 
aaa effects 14 \, for its use of a bin Similarly, if we ('(f,. 28 5.1185+- (0.955 & 1074) /f Cf 
artial wr ; ~ ; ; . . 

denne 

. : S P if > ape P Irmative P (O91) 3 . 

| 1.000.000 f,/(28/n f this is <60. oO that if the econd alt rnative of (‘ 21 hold the n 
th J 60) otherwise. (9.16) mail to the 7th destination should go into residue. 
lue, If the first alternative holds in (9.21) then 
9.1] 

hen the letter should be charged see 9 11 ; on 7 

81 50 1’ for the sweep of its bin Adding together C'(14 4.21) nts Cit... O.tzol \/ 0.1765+ (.04345 \/ 

_ these partial costs and defining ; ' 
9 19 and (from (9.21) and the fact 1,000,000 /,<60 
CT... N cost to be assigned to a letter to the 

OSE II ith destination, if such letters are t.459 & 107° Sy, 57.746 10 
sidu sorted directly on the primary, 
Ou 9.17 129.1 <2=(1/yf,) $225.3. 9.22 
sibly | we have . 
_— iii 1” =a /\ From (9.19) and our expression for C(f,, %)mm, We 
ins C(fi, u O.lébor \4 SI oy then have 
ality 0.1765 Q 18 
as % Nin Cf 0.04345 r—0.955 & 10742" 


We have now developed the formulas needed to 1.8955: 
nvestigate whether or not mail to a given destina- , 


rmiy | tian a — Pig mie , , 
m should be put into residue. It is convenient to pic tynetion of x is increasing (i.e., its derivative is 


“Y divide the destinations into three classes: positive in the range (9.99) and is Oat ez 129.1; 

| Class 1. 1,000,000 *,<60, thus it is positive throughout the range (9.22), which 
it is Class 2 60— 1.000.000 f, < 1.200. shows that even if the first alternative of (9.21 holds, 
oul Case 3 1 200< 1.000.000 7 mail to the ith destination should go into residue. 
en il hus in our minimum-cost system, a// mail to class 1 
tina-| In each class, the “into-residue” cost C(f,) should be = stinations should go into residue. 
ding} compared with the “‘not-into-residue” cost C(f,, For our particular mail distribution (see app. A) 
plex, evaluated at that value of the system design param- the class 1 destinations are destinations 901 to 1,600, 
find) eter n, which minimizes it , and these 700 destinations receive 18,080 letters/hr. 
1 b Analysis of class 1. Since the first alternative of These letters lead via (9.6) to a cost of 

first vs holds, 9.12) vields 
Ss wi 2 (.353 1S,080 

('(f 2 O76 0.955 10 f 9.19 


: associated with loading complexes, of 
Since the first alternatives of (9.15) and (9.16 apply, 


SUC- 

9.18) vields oa 
bins yield 14 18.080/1.200) 4-14 700 
114 eee = nee 
’ ( (fon O.1/65 nH 2 674 10 nt 0.1765, 

9 20 associated with residue bins and secondary bins, 
yn is . ; 
; and of 
XS sothat OC (f..) on 0.1765— 2.674 * 10 n> 
COSI 
81.50 1S.080/60) + 81.5 & 700 
irom which if follows that €' ~R is a decreasing : 


lunction of ». for ; 
associated with sweeps of residue bins and secondary 


n?<1.515X 107*/f;, bins. Adding these, we find that in our minimum- 
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cost svstem class 1 destinations involve a total cost 
of approximately 


$104,400. 9.23 


Analysis of class 2. 


of (9.8) holds, (9.12) vieids 
C'(f 14,10°7 3.434 Y 24 
The first alternative of (9.15) holds, but either alter- 


native of (9.16) may apply, so that (9.18) yields 


CU ten 0.1765n (2.674 LO nl 0.1765 
(ifn 1680/10°7,). (9.25 

Ct. 0.1765n 0.392 10 (nf 1.5348 
if n |. 680/ 10° Q 26 


As in the analysis of class 1, the function given in 


(9.25) is a decreasing function of n, for 
n, <O.1231/yf 9.27 
and is increasing for higher values of 7 On the 


other hand, the function given by (9.26) is a decreas- 
ing function of n, for 
n, <0.04713/ yf 9.28 
and is increasing for higher values of » 
For further analysis it is convenient to divide 


class 2 into subclasses, depending on the relative 


positions of 7 0.1231/\f; (where the function 
given by (9.25) has its minimum), ;=0.04713/,/ 
where the function given by (9.26) has its mini- 
mum), and n,;=1,680/10°/; (where the formula for 
C'(f,, n,) changes from (9.25) to (9.26 This sub- 
division leads to 

Class 2a, 60< 1,000,000 7 186.3, 

Class 2b 186.3< 1,000,000 1.200, 
corresponding respectively subject to 60 10° 7 
1.200) to 

O.04713/\f O.1251 7 |. 6S0/10°7 Y POn 
O.04713/\/ 1 .680/10°/ O.1231/\f Q 2Oh 

Analysis of class 2a Here C(f,. n must arise 

either by setting n,=0.1231/\f, in (9.25), vielding 
0.1765 0.04545, ,/ 

or by setting n,;=1,68010°/, in (9.26), vielding 
1.768] 296.5/10° 

It turns out ‘that the first form is the smaller, so 


that, combining the first form with (9.24), we have 


Since the second alternative 
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Clr. C'(f 0.043457 14.10" 3.258 
9.30 

This function of f; is positive for 
60 1 O00.000 7 168.2. (9.3) 


which shows that mail to the 7th destination shoy// 
vo into residue if (9.31) holds. 

For our particular mail distribution, destinations 
526 to 900 are the ones obeying (9.31 372 
destinations 10,250 letters/hr, leading to g 


cost of 


these 
receive 


10.250 


associated with loading complexes, ol 


12 10.250°1.200 14 375 
associated with residue bins and secondary bins 
and of 

81.50 $0. P5060 S150 10. 250/60 


associated with sweeps of residue bins and secondary 
bins. Adding these, we find that in our hypothetical 
minimum-cost system the mail to destinations obey- 
931 


Ing has a total associated cost of approximately 


$143 .500 Q 29 


sefore proceeding further, we will write out ex- 
plicitly a principle which will be helpful in much of 
the following work. The correctness of the prineiph 
could be proved analytically in situation in 
which we invoke it, but such proofs would be lengthy 
and and 
with stating the principle (which is intuitively evi- 
dent) anyhow Informally, the principle simp 
recognizes the fact that bins for destinations receiving 
a good deal of mail should appear rather frequently 


each 


repetltlous, SO We content ourselves her 


in the primary of our minimum-cost system (1.e., 1 
will be small), bins for destinations receiving less mail 


should appear less frequently, (i.e., n,; will be larger 


and there is a “‘threshhold ¢” such that the destina- 
tions whose mail is put into residue are precisely 
those which receive fewer than ¢ letters/hr. Formally 

Suppose it is best in our minimum-cost system 


that mail to the 7th destination not go into residu 
and suppose that is the best value of 
Then for any other destination with /,;>f;, it 1s 
best that mail to this destination also not go into 
residue, and furthermore () a, 7 Q 33 
Returning to the analysis, we now consider th 
class 2a destinations not obeving (9.51); Le., those 
for which 
168.2 < 1.000.000 7 186.3 Q 34 
For these destinations, the function of f,; given in 


9.30) is negative, but we cannot automatically con 


5.258 
(9.30 


9.3] 
hould 


it ions 


LO 4g 


bins 


L eX- 
‘h of 
cipk 
nD in 
gthy 
her 
eVi- 
nply 
Ving 
ntly 
mail 
ver 
Lina- 
isely 
ally 


stem 


au 


n Ww 


con 


lude that mail to these destinations should not LO 
{ 


» residue. The difficulty is that we treated n; as 


jnt . . 
P eontinuous variable in finding C'(f;, n min» whereas 
‘. fact n; must be a positive integer. We therefore 


define (for fixed 7 


for all positive 


(9.35 


minimum of C(f,;, 7 
integer values of » 

fa, a will be one of the 
gers between which the n vielding 7 oo 


the value of »,; yielding C' 


inte mir 
lies. ke tee 
We begin at the “top” of (9.54), with 
1.000.000 f,— 186.3 9.36 
The value of n,; vielding C(f,, % min IS 


so that (opt is either 9 or 10. From (9.24) we have 
C'(f C1186.3 X 10" 3.51 
Substituting n,;=9 and n;—10 into the applicable 


choices of (9.25 and (9.26) (since 1,680 10°F q 


9 goes into (9.25) and n;—10 into (9.26 we 
yytain costs of approximately 
3.34 for n Q, 3.51 for n LO, 


so that (isjopp=9 and mail to destinations obeying 
9.36) should not go into residue in our hypothetical 
minimum-cost system. We note for future reference 
that, by (9.33), all destinations in classes 2b and 3 
should not have their mail sent into residue, and that 
2. 
through (9.34 
Throughout (9.34 


tl roughout these classes (vi , 
We now continue ‘‘down” 
sidering lower values of / 


con- 


/ () 1231 \ / +] 
so that (7 opt Is either 9 or 10; also, throughout 9 34 
ve have 1.680/10° 4 so that 7 Q must be sub 


stituted into (9.25), 7 10 must be substituted into 


9.26), and the results compared with the result of 
24 Since this vields 
C(t, 9 1.7650 297 .1/10°f 
ex F.. 16 , 2QOS 392/107 
ve find (using (9.24) also) that mail should go into 
esidue for 
168.2 < 1.000.000 7 169.6 9.37 
end should not go into residue for 
169.6 L.OOO O00 7 186.3. Q 3S 
and also that 
n 8) in (9.38 9.39 


For our particular mail distribution, there are no 
destinations in the narrow range (9.37 The cost 
associated with destinations obeying (9.38) will not 
be found here, since it can more conveniently be 
combined with the cost found for class 2b. 

Analysis of class 2b. We know that mail to these 


destinations does into residue. Also, 


not go we 
know that at the “bottom” of class 2b (i.e., 10°F, 
186.3) we have (n,) 9, and from (9.33) we can 


opt 
expect that as we search “upward” through class 
2b (i.e., examine successively larger values of f,) 
we will reach a poimt where (7 opt changes from 9 to 


8. Our immediate aim is to find this turnover 
point 
We find that for 186.3<10°f,<186.7, C(f;, 8) 
and ((f;, 9) are both given by (9.25), so that — 
Cf, 9) = 1.7650 + (297.1/10°f 
C'(f;, 8) = 1.5885 + (334.3/10°Ff, (9.40) 


Throughout the range C(f;, 9) is smaller, so that 


Ns)ent=9 for 186.3 < 10°, < 186.7. 
For 186.7<10°7,<210, C(f;, 8) is still given by 
9.25) but C(f,, 9) is now given by (9.26), so that 
C:, 9 3.1253 + (43.56/10°F (9.41 


Comparing (9.40) and (9.41), we find that 


for 186.7<10°f,< 189.4 


but 


10°f; immediately above 189.4 

We can continue to search upwards through class 
2b, looking next for the turnover from (n,)op.=8 to 
(Wi) opt 7. ete. The results of this search are given 
in table 1 which includes the class 2a destinations 
obeying (9.38). The quantity 10°/,; has been rounded. 

The first three columns of table 1 are independ- 
ent of the particular mail distribution. The first, 
third, fourth, and fifth columns can used to 
derive a figure in the following way: For 
each row, the fourth entry shows how many destina- 
tions are involved and the first entry enables us to 
find how many bins to each destination appear in 
our minimum-cost system. If (9.25) holds, then 
these bins are each swept once an hour, whereas if 
9.26) holds then the number of sweeps involved is 
the same as that used in deriving (8.7); we can 
the fifth entries of the rows involving (9.25) to find 


for values of 


be 


cost 


use 


the numbers of “extra sweepers’ involved (see 
8.11)) and then multiply this contribution to AS 
by 14,670 in accordance with (8.12 Finally, we 


can use the first and fifth columns, together with 
9.4) and (9.6), to assign a loading complex cost to 
each row. The result® is a total cost of about 


gral numbers of bins witl 
9.42 


Different methods of associating int 


which are not divisors of 28 were tested and found to change 


values of n; 
only negligibly. 








$401,200 (9.42 ing of a single loading complex followed by 1609) of « 


bins, one for each of the 1,600 destinations. of a 
associated with class 2b destinations and the class 2a The argument used to derive (8.3) shows that for | loa 
destinations obeying (9.38). this system L= 28, so that tion 

TABLE 1, Optimiz ng m n class 20 ] eacl 
aL=0 10.1) | and 
Niopt 10° / Formula Destinatio Letter sort 
thus there are 28 subsystems, so that ! ‘I 
9 170 to 187 2 401 to 525 eh) of n 
u ISS to TSU 1 26 > IR B00 RAI 
8 190 to 209 9.25) 451 to 490 7,700 AB - 1,6 4() £4,000. 10.2 and 
8 | 210 to 217 1.26) 441 to 450 2, 100 sort 
7 » 0 234 1, 25 21 to 440 1 Ow 
ats to a e21 to am Assuming the 1,000,000 letters/hr divide equalh 
7 | 240 to 254 9.26) 401 to 42 1, 900 4 rage) yr the 2S s svstems » Sn, 
ime Oa ae oe oe 0 (on the average) among the 28 sub ystems, we find 
i «280 to 305 9, 2¢ 351 to 380 &, SOK that i bin corresponding to one ol the 114 most so t 
5 | 306 to 335 2 $41 to 350 , DM . . 
336 to 384 ) 2¢ $11 to 340 10, 800 frequent destinations will receive OU Or more let 
{ 5 to 41 , ‘ ters hr, whereas a bin corresponding to one of thi 
3S. Oo 4G 1.2 iM) l : : } 
1 | 420 to 516 1.26) 265 to 299 12¢ | ,600— 114 1486 least frequent will receive fewer | and 
$ | 517 to 559 9. 25)| 253 to 264 SUK , : There are 
$ 560 to 786 ) 2 203 2 2 TH than b0 Thre re are 
2 787 to S3Y », 25 200) A 2 4 
2 840 to 1, 200 ) 154 1” 46), 40K 


| 486 * BS $1 600 so t 
otal number 

bins in the system which correspond to these las 
1.486 destinations, and by (5.3) these bins give ris and 
to 41,600 stacks/hr. These last 1.486 destinations 
receive only 27.6 percent of the total mail 276,00 
letters/hr and so, 7} all stacks consisted of 60 letters On 


Analysis of class 3. Here the different possible 
alternatives in (9.15) and (9.16) lead to three possible 
formulas for C(f,;, 


Cf, n 0.1765n (2 674 10 n f 0.1765 the basis on which (8.5) was derived), they would Seh 
(if n~<1,680/10°Ff,), 9.43 give rise to only 
C fi, n 0.1765n,;+ (0.392 « 10 (n,f 1.5348 276.000 /60— 4.600 stacks/hr. 
Gf 1,680/10°F <n; <33,600/10°7;), 9.44 

Thus we have 41,600 — 4,600 = 37,000 extra stacks hy 
C(fi, 0.1765n¢+ 1.546 leading via (5.2) to lett 
Gif nm; >33,600/10°F,). 9.45 a re 
, ” 37 000/180 ~205 all 
We can work “upwards” through class 3 just as we rela 
did in class 2b, beginning with (nj)op.=2; the only additional sweepers, of which two were taken into | the 
new complication is that possible transitions to g¢eount in passing from (8.5) to (8.7 So , pur 
(9.45) must be allowed for. It turns out that the nur 
change from (Niopr=2 tO (Niopr=1 occurs for AS— 202 10.3 des 
10°f, = 1,614 (.e., at destination 119 in our particu- leac 
lar mail distribution) and that the cost associated and by (8.12 | 
with class 3 destinations is whi 
- ( 14 14,000 14.670 203 3.594.000 SO 
$1,397,000. 9.46 2 , Bs O94, 1.4) | are 
lef 


Finally, we add up (9.23), (9.32), (9.42), and (9.46 
to obtain (after rounding off) a total of 


C*  —=$2,046,000 9 47 c= AC 9,399,000, 
as the approximate cost of our hypothetical mini- and by 9 47). the 
mum-cost system. This is then our theoretical — 
minimum cost, to be used as a vardstick in dealing ( Hh Be. 1629 
with proposed sorting configurations. It is signifi 
cantly larger than (8.8). Thus on a cost basis the simplex scheme should 
The general approach used in this section consti- definitely be rejected The excessive cost comes 
tutes the “follow through’? method mentioned in primarily from . personne! i.e., from the second ‘ \ 
section I, summand in (10.4 = 
10. Simplex Scheme : Th 
11. Christmas Tree Scheme 9. 
The first sorting scheme we examine is also the wns 
simplest. The sorting system consists of a number This might also be called the ‘“‘square root’ 
of essentially independent subsystems, each consist- scheme In our case the square root of the numbe r 
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of destinations is 41,600=40. The system consists 
of a number of subsystems, each containing a single 
loading complex and 40 bins." The 1.600 destina- 
tions are divided into 40) groups ol 40 destinations 
each: the primary sorting of mail is done by group, 
and the mail to each group is then given a secondary 
sort by individual destinations. 

The key to the analysis is the fact that each piece 
of mail is dealt with twice. By applying (8.3), (8.4 
and (8.5) to both the primary and the secondary 
sorts, and adding the results, we have 


L>56, B> 1,680, S> 186. 
so that 
AL > 28, AB> S840, AS>9] 
and by 8.12 
AC'> 1,702,000, 
so that by (8.8) and (8.13 


('>3.463,000 


and 
C~C*,)/C* 


1 min 


O07 
O09. 


basis we can reject the Christmas Tree 
Again personnel costs are the major factor 


On this 
Scheme. 


12. Residue Scheme 


In this system the primary sort consists of sorting 
letters directly to those destinations which receive 
a relatively large fraction of the mail, while dropping 
all letters to the less “‘frequent”’ destinations into a 
relatively small number of residue bins, which are 
then given a secondary sort by destination. The 
purpose of this maneuver is to avoid having a large 
number of sweeps (corresponding infrequent 
destinations) resulting in small stacks; such sweeps 
lead to excessive personnel costs. 

The system is determined by stating definitely 
which destinations are to be considered “infrequent r 
so that mail to them goes into residue) and which 
are “frequent.”” We therefore 
define a system design parameter 


to 


Lo be considered 


frequent if 
(12.1 


the «th destination 
infrequent if 10°%;<t; 


l threshold: 
10°F >t. 


Is 


the value of ¢ is to be chosen so as to minimize the 
cost of the system As before, 
}=fraction of mail which goes to the 7th destination. 


12.2 


We shall determine the optimal value of ¢ by two 
methods (in order to check their agreement First 
we apply the “follow through” approach of section 9. 
lhe residue scheme certainly satisfies (9.1 and 
4.3); 1f well-designed, it will very nearly satisfy 
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(9.2) and (9.5) as well. We will therefore proceed 
as in section 9. The number of letters/hr to the ith 
destination in each of the 28 primary subsystems is 


1,000,000 7,/28, 
and if this quantity is >60 then by (5.4) there 
would be no reason to put mail to the ith destination 
into residue. Thus the only destinations about 
which there is any question are those for which 


(1,000,000 7,/28) < 60. 


We can split these destinations into classes, according 


as 
1,000,000 7, <60, (12.3a 
or 
60< 1,000,000 7, < 1,200, (12.3b) 
or 
1 ,.200< 1,000,000 7, < 1,680. (12.3¢) 
For cestinations obeying (12.3a), (9.12) yields 
CU 2.076 + (95.50/10° f,), (12.4a) 
while (9.18), since n,;=1 in the residue scheme, 
vields 


(2674/1 O°f 


Similarly, for destinations obeying (12.3b) we have 


C(f) =3.484+ (14/10°f, (12.4b 
C(f;, 1) =0.353 + (2,674/10°f,). (12.5b) 


For destinations obeying (12.3c), it is not clear in 
view of (5.5) what a letter going into residue should 


be “charged” for the sweep of its secondary bin; 
it should be at least 14/1,200, however, so that 
see (9.25 we have 

C'(f 14/1,200) +-3.434=3.451, 12.4¢ 
and 

C(fi, 1) =1.71138 + (392/10 (12.5¢ 
We find that C(f ('(f;, 1) is negative (i.e., it is 


less expensive to put mail to the ith destination into 
residue) for 10°;<863 and is positive for 10°F, 
864. Thus the optimal ¢ is 
lopt S64, (12.6 

which for our particular mail distribution corre- 
sponds to destination 196. 

Next we apply the ‘add up” method 
8.12 First we define two relevant 
depending on f, by 


based on 
quantities, 


N,=number of frequent destinations, 
\’,=number of letters/hr to frequent destinations, 


so that (since our situation involves in all 1.000.000 








letters/hr and 1,600 destinations 


1,600—N, 
1.000.000 


the number of infrequent destinations, 
V.=number of letters/hr to infrequent 
destmations. 


The number AL of extra loader complexes needed 
for the secondary sort is given (using (8.2)) by 


AL = (1.000.000— V,,) /36.000: 12.7 


more precisely, AL is the integer next above the 
right side of (12.7 The infrequent destinations 
are divided into AL grou ps whose total expected 
hourly mail volumes are approximately equal. 
The sorting system contains 28+ (AL 
systems, each consisting of a loading complex 
followed by a number of bins. Each of the 28 sub- 
systems needed for the primary sort contains 
bins for the frequent destinations and one residue 
bin" for each of the AL groups of infrequent desti- 
nations. Each of the AL subsystems used for the 
secondary sort, sorts the entire primary residue of 
some one of the groups; thus these AL subsystems 
together contain one bin" for each infrequent desti- 


sub 


nation, and thus contain 1,600— N, bins in. all 
Therefore 
B=28(N,+(AL 1.600—N 
28(AL) +27.N,+ 1,600, 
AB= B—840=28(AL) +27N,+-760 12.84 
Having found AZ and AB, we must find AS. The 


frequent destinations receiving 60 or more letters/hr 
on each of the 28 primary subsystems, require no 
more sweeps than was assumed in calculating the 
minimum cost. The same holds for secondary sweeps 
of mail to the infrequent destinations. The 
sweeps therefore arise (a) from sweeping the residue 
bins and (6) from the primary sweeps of frequent 


eartra 


destinations receiving fewer than 60 letters/hr on 
each of the 28 primary-sorting subsystems; the 
frequent destinations deseribed in h are those 


receiving fewer than (28 60)= 1,680 letters/hr in 
all, and thus those whose mail frequencies are less 
than 0.168 percent. In our actual mail distribution, 


114 destinations have frequencies of 0.168 percent 


or more, and these destinations account for 72.40 
percent of the total mail. Thus 
N, 114 number of destinations ol the type ly 
12.9 
V,—724,000=hourly volume of mail to these 
destinations, 12.10 
The destinations described in (+) thus lead to 


approximately 


114)/1S0 sweepers 


of which ((V 724.000) 60)/1S0 were accounted fo 


in ecaleulating (8.8 The number of extra Sweepers 


due to source (6) is therefore approximately 


1.680.N \ 532.500) /(60) (180 


The residue bins together receive 1,000,000—Y. 


letters/hr and each receives more than 60 letters hr ! 
so that the number of extra sweepers due to soure 
a) is approximately 

LSO 


1.000 000 / 60 


adding the last two expressions gives 


AS 1.532. 500—2V 1. 6SO.N,)/(60) (180 


12.1] 
substitute (12.8 


We now and (12.11 


obtaining 


into (8,1? 


AC 13, 100( AL) + 2,660.\ 2. 717V 4 2.092.500. 


) 


Theerrorinvolved in using (12.7 
13,100 in (12.12) and thus 1 
percent in (¢ bp PRY 5. 


12.12) getting 


ing to an error of at most 
to an error of less than | 


We therefore substitute (12.7) into 
AC 2 660 A 2081 V,4+2.456.400. (12.13 
To minimize AC’ quickly (we omit the rigorous 


justification of the following method 
differential of the right side of (12.13 


equate 
to zero: 


2 660(dN ,OS1(dV 0) 


Since NV, increases in steps of size 1 (each step involy- 
ing shifting the status of one destination from 
“infrequent” to “‘frequent’’), we set dN,=1 an 
obtain 

dV 2 660) 5.081 863 letters hr 12.14 


Since the increment in V, due to one extra “frequent 
destination is simply the expected hourly volume of 
mail to that one destination, (12.11) tells us that for 


the approximate minimization of AC, the last of th 


“frequent” destinations should obey 10°, = 863, 80 
that (see (12.1 

f St) *® 1? 15 
in near-pertect agreement with the result 12.6 
obtained by the ‘ follow through” method For ou 
particular mail distribution this ‘cutoff’ occurs 


after obtam- 
196 from 
find 


between destinations 1°6 and 197, and 
ing the value of V 
appendix A and 
rounding) that 


corresponding to .V 


substituting into (12.14) we 


Af 140.000 
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ane 


Using (12.11) we find that the cost of personnel is 
again the main cost factor 

The optimum given by (12.15) (i.e., given by a 
eytoft’” between destinations 196 and 197) is not a 
very sharp one; the cutoff can occur as low as about 
destination 165 or as high as destination 225 without 
raising (C—CRin)/Chin to more than 0.5 percent 
above its minimum. 


13. Multiple-Input Schemes 


In these schemes, the sorting system consists of a 
number of identical subsystems, each receiving its 
input from a number of loading complexes arranged 
in series. More precisely, each subsystem consists 
of a first loading complex followed by a first row of 
hins for some (but not all) destinations, then a 
second loading complex’ followed by a second row 


t 


, of bins for some destinations (not necessarily the 


same ones as in the first row and finally a 
last loading complex followed by a last row of bins 
for all 1.600 destinations. 

Consider some particular sorting system of this 
type; eall it system 1, and let 


\f/=number of subsystems of system 1. (13.1 


Let system 2 be obtained from system 1 by replacing, 
with residue bins, the last 700 bins of the last row of 
each subsystem these bins correspond to destina- 
tions 900 to 1600, which for our mail distribution are 


he ones receiving fewer than 60 letters/hr The 
esidue then requires a secondary sort. Let 
R=number of residue bins/subsystem for 
svstem 2 13.2 
b) 


We will compare the costs of system 1 and system 2; 


et L,, B,, S,, C), L., Bo, S,, €, denote the respective 
values of L, BLS, C, for the two systems. 


Destinations 900 to 1,600 receive about 1.8 percent 
18,000 letters/hr: this is less than the 
6,000 letters hr capacity of a loading complex; and 
SY) milly one complex is needed for the secondary sort: 


The 700 7 bins used for destinations 900 to 1,600 
ithe last rows of the subsystems of system 1 are 
replaced in system 2 by RM residue bins plus 700 
bins for the secondary sort, and so 


a i 


B,— B,=700M— (RM -+-700 (13.4 

The above-mentioned 700 .\/ bins for system 1 
each received fewer than 60 letters/hr, and so to- 
gether required 700 M sweeps/hr. For system 2, the 
700 secondary bins together require 700 sweeps/hr. 
As for the RM residue bins, ? is chosen as small as 
possible, so that each residue bin (except possibly for 
the last one on each subsy stem) receives at least 60 
letters/hr: the residue bins receive 18,000 letters/hr 
and thus require approximately 


18,000/60 ~300 sweeps/hr. 


Thus we have, using (5.2 


S, —S. = (700.M— (300 + 700))/180. 


The analog of (8.12) which applies to our situation 


Is 
(',—C,= 12,700(L,—L,) + 14(8,—B,) + 14,670(S,—S, 
(13.6 
and so, by 13.3 (13.4). and (13.5) we have 
( ( 67 000.\17— 14RAT— 104,000. (13.7) 


Since FP is chosen as small as (5.5) permits, we have 


RM = 19,400/1,200 = 16, 


and so (13.7) can be rewritten as 


67 OOO.AS— 104.000. 13.8 

One rather sweeping conclusion which can be 
drawn from (13.8) is that no multiple input scheme 
with more than Ole subsyste mn (1.e., with M 1) can he 


optimal. ‘To prove this, take the system in question 


as the “system 1” of the above discussion. Since 
\J~>1, (13.8) shows that C;—(C,>0 and thus that 
( (’,; since system 2 costs less than system 1, the 


latter cannot be optimal. We could‘use this argu- 
ment to eliminate multiple input schemes with 1/7 >1 
from further discussion, ‘f we were going later to 
examine the “systems 2.”’ Unfortunately, the 
analysis of such systems, which combine elements of 
the multiple input and residue schemes, appears too 
complicated to be attempted here (the difficulties of 
even relatively simple systems of this type will be- 
come apparent in sec. 14 and app. C We can, how- 


ever, use (13.8) to eliminate all multiple input sys- 
tems (1.e., all “systems 1”) with M>?7. To do this, 
we note first that (,>C,,,, so that 

AC’ &s a c,— Cs; 
together with (13.8) and \/>7, this implies that 


AC ( (’, = 67 000 AZ— 104,000 > 365,000, 


so that system | is at best negligibly less costly than 
footnote 13, p. 96) the 


actually more costly; see 








system to be found in section 14, for which (see 


(14.24)) 
('=366,000 (13.9 
Now we want to examine multiple-input schemes 
with <6 (i.e., with six or fewer subsystems). It 
is convenient to define 


Ny=number of destinations receiving fewer 
than 60M. letters/hr, (13.10 
Wy,=number of sweeps associated with these 


destinations in calculating (%,,. (13.11 

Kach such destination receives fewer than 60 letters 
hr (and thus requires at least one sweep/hr) for each 
subsystem, so that together these destinations require 


at least MNy sweeps/hr for the entire system. Thus 
AS > (MNy— Wy) /180. (13.12 
For WM=6 we find MNy ~7,620, Wy, ~2.140, so 
that, by (8.12), 
AC'> 14,670 * 5,480/180 ~ 447,000 
(for 1=6 (13.13 
By comparison with (13.9), systems with M=6 are 
eliminated.“ For M=5 we find AINy ~6,150, 
Wye ~ 1,925, as well as 
AB= B—840>5 1,600—840=7.160; (13.14 
thus by (8.12) 
AC’> 14,670 * 4,225/180+ 14 7,160 ~ 445,000 
for M/ 5 is.ko 
By comparison with (13.9), systems with M=5 are 


eliminated. 
These relatively simple arguments are not ade 


quate to deal with the schemes with 1/<4. Con- 
sider first the situation \J=4. Here Ny ~4,720, 
Wy ~1,710 so by (13.12 

AS > 3,010) 180, 13.16 


If each subsystem contains 
eight or more loading complexes, then 


Two cases are possible. 


AL=L—28>4 x 8—28=4, 13.17 
AB=B—840>4 1,600—840—5,.560, (13.18 
and by (13.16), (13.17), (13.18), and (8.12), we obtain 


AC’ > 374,000, 


so that in this case the system can be eliminated by 
comparison with (13.9 The other possible alterna- 
tive is that each subsystem contains no more than 
seven loading complexes. Seven loading complexes 


rhe same argument, if applied to M=7, vields AC 2,000 
result than the one found earlier (the d splay first alx 13.9 
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can handle at most 


252.000 letters/hr. 


7X 36.000 
and each of the four subsystems must handle 


1,000 000 /4—= 250,000 letters/hr. 

so seven complexes subsystem are needed, all work- 
ing close to capacity; i.e., the input of each loading 
complex in a subsystem must be nearly all dropped 
in the following row of bins in order for the next 
loading complex to have a nearly empty conyeyo 
belt. This clearly requires at least 200 bins after 
each of the first 6 complexes in each subsystem 
the seventh loader is followed by 1,600 bins (on, 
for each destination) so that 


AB=B—s840>4 6 200 1.600) —840= 10.360 


13.19 
8.12) we obtain 


from (13.16). (13.19) and 


(’>390.000. 


so that the system is eliminated by comparison with 
(13.9 Thus syste ms with M t are eliminated. 
We have not been able, within the limits of time 
and effort reasonably assignable to this particular 
point,'* to devise a mathematical proof that multipk 
input systems with 1/—1, 2,or 3 can be eliminated 
because of excessive Cosi The difficulties en- 
countered, and the reason for their occurrence for 
small values of AM, ure discussed in the next par- 
agraph. Nevertheless, we strongly believe that 
such svstems ure excessively costly This belief is 
based on auxiliary calculations which will not be 
reproduced here, and also in general on our ex- 
perience with the othe parts of this study. Roug! ly 
the situation is this: In each subsystem, most of the 
mail must be dropped out fairly soon after it enters 
the subsystem, for otherwise the input from many 
of the loading complexes would be cut substantially 
below 36,000 letters’ hr, so that number of 
loading complexes cost eXCESSIV! 
would be required to pass the required 1,000,000 
letters hour into the system Such an early dropout 
of most mail, however, would require a large numb 
of bins, many of which would receive fewer tha 
60 letters hr and thus (see (5.3) and (5.4)) involv 
“inefficient sweeping’; the cost of bins and sweepers 
especially the latter) would then be 
The difficulties for small 
Vf (i.e., for systems with a small number of sub- 
systems fact that 
subsystem must contain a relatively /arge numbet 
of loading complexes; for 1/=3, for example, eac! 
subsystem contains at least 9 such complexes, whil 
for \J=1 we have a single long subsystem with at 
least 28 loading comple xyes in series The first caus 
of difficulty is that the formulas involved in the 
analysis of a subsystem with NW loading complexes 
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become more and more complicated as 1” increases. 
See Appendix B for more details. Kor n>9, for 
example (i.e., for M<3 . it is still possible to use 
these formulas in making calculations for any one 
iven subsystem, but is extremely difficult to use 
them in comparing a large number of possible sub- 
evstems; this latter problem is of course, the one 
which actually arises in our work. The 
cause of difficulty is the enormous number of systems 
to be compared, when the value of » is high (i.e., 
when M is low); this number is in fact 


second 


n—1)° (/=number of destinations 

and since d= 1,600 it is clear that even if 99 percent 
of the possible systems could be eliminated on some 
common-sense grounds, the number remaining for 
(if m is moderately large) would still 


analysis be 


astronomical 


14. Multiple Input and Residue: Double 
Loading 


The “multiple input and residue’ schemes, com- 
binations of two of the proposed schemes which 
we have analyzed earlier, are like multiple input 
schemes except that some of the bins may be assigned 
to residue; this residue then requires a secondary 
sort. We therefore speak of “primary subsystems’”’ 
and “Secondary subsystems.” 

There is a great variety of subclasses of this type 
of scheme, and within any one subclass the analysis 
required to determine the optimal choice of the 
relevant parameters appears to be quite difficult. 
We will make a detailed analysis only of a relatively 
simple subclass; the result turns out to be a system 
so nearly optimal (14.25)) that investigation 
of more complicated systems isclearly not worthwhile. 

We consider “‘multiple input and residue”? schemes 
which are determined by two parameters, 7 and k 
with ;<k), in the following way: Each primary sub- 
system contains two loading complexes The first 
omplex of each primary subsystem is followed by 
bins for each of 


See 


the “type 1” destinations 1,2, .. . , J, 14.1 
and also by re siduce bins for 


1.600. 
(14.2 


‘type 3”” destinations b+-1. 


Thus that part of the first complex’s input consisting 


of mail to 


the “type 2” destinations } ba ee hs (14.5 


| 


does not get sorted until after the conveyor passes 
the second complex. The second complex of each 


Primary subsystem is followed by bins for type 1 


97 


and type 2 destinations, and residue bins for the 
type 3 destinations. We wish to choose j and k so 
as to minimize the cost of the system. 

We will use the notation 


J:=Tfraction of mail to destination 7; (14.4) 


the destinations are so ordered that 


> fs > 


>} . 
< 1, 600+ 


We also set 
fraction of mail to the first 7 
destinations. (14.5) 


The input to the first loading complex of each 
primary subsystem is the usual 36,000 letters/hr, but 
only a fraction F';+-(1—F;) of these (corresponding 
to types | and 3) get dropped before the second com- 
plex, so that the input to the second complex is 


36.000 (1+ F FF.) letters/hr 


and the input to each primary system is 


36,000 + 36,000 (1+ F)—F;, 

36,000 (2+ F,—F;,) letters/hr. 
Therefore the required number of primary S ubsyste mis 
is about 


1 .000,000/36.000 (2+-F,—F,) 
(14.6) 


so that the number of loading complexes for the 
primary subsystems is about 


56/(2+-F,—F, 
while the residue of 1,000,000}(1—F;) letters/hr re- 
quires 
1,000,000 (1—F,)/36,000 =28 (1—F, 14.7 
loading complexes for secondary sort. Thus 
L = (56/(2+-F,—F +-28(1—F’, 
or AL= L—28 = (56/(2+- F;—F,))—28F, (14.8) 


The type 3 destinations are divided into groups as 
for the residue scheme 12); the number of 
these groups is given by (14.7). The secondary sub- 
systems together contain one bin for each type 3 
destination, or 1,600—k bins in all. Each primary 
subsystem contains 7+ bins for separate destina- 
tions (j bins before the second loading complex, & 
bins after it), and also two residue bins for each 
group (one before the second complex,"one after it). 
Hence, using (14.6) and (14.7 


(sec, 








B= (1.600—k) + (28/(2+- F J 


st) 


AB= B—840 =~ (760—k 


28 (7 tft 56 | l 14.9 

Having found AL and AB, we still have the more 
complicated task of finding AS. The sweeps of bins 
of the secondary sort require no 
The extra sweepers are required (a) in sweeping the 
residue bins,'® (6) in sweeping those bins (if any 
before the second loading complex which receive 
fewer than 60 letters/hr, and (c) in sweeping those 
bins after the second complex which receive fewer 
than 60 letters/hr. There are 


ertra Sw eepers 


1.000.000 (1—F letters/hr 


the contribution 


and a2 


coing to the residue bins, and so ! 
of (a) to AS is approximately (using (5.1 
1,000,000 (1—F (60 180 14.10 
The analysis of the contribution of (6) to AS 
depends upon the fact that, for our mail distribution, 
only destinations 1 to 114 receive 60 or more letters 
hr out of 36,000 letters/hr (the input of the first 
loading complex of a primary The 
contribution of (b) to AS is 


subsystem 


() if ysl, 14.14 
If 7 >114, then the j7— 114 destinations 115, 116, 

J have their bins before the second loading complex 
receiving fewer than 60 letters/hr, leading (via 
(14.6 to approximately 

28 (j—114)/(2+-F fF.) sweeps hr 

Not all of these are extra sweeps, however; the volume 
of mail involved is 


1.000.000 (F Fs) /(2+F hI 


letters hr, 
which in calculating the minimum cost would receive 
60(2 + F fk 


1.000.000 (F I sweeps hr, 


Thus we have a contribution of (+ 
mately 


to AS of approxi 


[28(7—114 1 000 000 (F Fy) /60)/180(2—+F h 
(if 7 >114). 14.12 
The contribution ot ( to ANS arises trom fio 
SOUrCES, I rst, there are the type | destinations (if 
any) which receive fewer than 60 letters/hr after the 
It i imed t I 
For 1 2+ J / 
ug! rat 


second loading complex of each primary subsystem | By 
To handle these, we define nnew variable mm, depend. 
ent on } and k, by \ / 


last destination for which 
36.000 (1 I FP) 60). 14.12); whi 


destination m 


The contribution to AS from the first source of («| Ay 
is then 

0 if m >); 14.14 
if m<j; then the (j—m) destinations m+1, m+9 


: j are the ones under consideration. leading Wilt 
Via (14.6 


to approximately 
28 (j—n 2+Hh I 
sweeps, of which (see the derivation of (14.12 Ne 
1.000.000 (F I | +f F.) /60(2+ F I hi 
fro 


Thus we have a contribution to AS from 
, of approximately mul 


are extra 
the first source of (« 


14 
28()—m 1.000 .000(F km)(1+-F I.) /60 whl 
IS0(2+-F I if m<) 14.15 N 
Second. there are the type 2 destinations whiel LO 
receive fewer than 60 letters/hr after the second com- 1 
plex 1.8.. the input from both complexes adds up to alN 
fewer than 60 letters/hi To handle these. wi 
introduce another new variable n. also dependent on pal 
and k. by 
destination 7 last destination for which 
36.000 (2+ Ff F.)f 60 14.16 viel 
The contribution to AS from the second source of SU 
is then ee 
() if n>k: 14.17) Th 
reasoning as in the derivations of (14.12) and (14.15 “a 
we find the contribution to AS from the second ny 
source of (c) to be approximately 
| 
I (I 7 1.000.000 (/ / Ye IF.) /60 el 
180(2-+-F I if j<n<h 14.18) #PI 
and to be approximately e 
ap] 
IS( bk —} 1.000.000(F KF) (2+F Fk.) 60 ma 
1S0(2+ F / if 7 14.19) “"s 
This completes the derivation of the approximat 
formula for AS : 
At this point it is convenient to split AC into two 
parts, y 
AC A, A. 14.20 in 
\ 


where A,¢ represents the effects of AL AB, and the nn 
contribution of (a) to AS, whereas A. represents the _ 


effects of the contributions of hb) and ( to As 


98 


‘VStem 


epend. 


14.14 


Mi } 


eading 


) from 


14.15 
which 
com- 
up to 
. Wwe 
it on 


4.16 
of 
$.17 


£.15 


Cond 


4.19 
mat 
. two 
$21) 
| the 


s the 
AS 


By (8.12) and (14.10 
4,0 =12,700(AL) + 14(AB 
14.670(1,000,000(1— F 60)(180 14.2] 
which, together with (14.8) and (14.9), yields 
4,C=[37.48 & 10° 1.234 « 10°)7 +- (3.931 & 105k 
14.78 10%)F 52.04 ™ 10°) F 
18.51 * 10°) (F 2—F,F (60)(180 
2+F;—Ff 14.22 
where the lower-order-of-magnitude quantity 


1.512 10)k(F I 

has been dropped from the numerator of (14.22 
\o single formula can be given for A.C. This 
pnecause Apt 


Is 


which is obtained (by multiplication by 14,670 
from either (14.11) or (14.12) (according as 7<114 
ol i14 the second of which is obtained (by 
multiplication by 14,670) from either (14.14) or 
14.15) (according as m } or m< } . from the third of 
which is obtained (by multiplication by 14,670) from 
either (14.17) or (14.18) or (14.19) (according as n >k 

Ss k orn | These diverse possibilities lead 
to 2X23 12 cases: the case-by “Case analysis Is 


quite complicated, and we relegate it to appendix C, 
riving here only the result: 

The « AC’ is approximately minimized for our 
particular mail distribution by choosing 


‘Ost 


p== $30, k= 264, 14.23 
vielding AC’ = 366,000 14.24 
so that ( Cate ae 1% 14.25 
The minimum is a rather insensitive one; if 7 is 


hosen any where between roughly 90 and 150 then 
assuming = properly Can) /Ce 
will be less than 0.5 percent above its minimum. 

The analysis given and continued in 
C) has employed the ‘add up” method only, 
deliberately avoiding any use of the ‘follow through” 
approach of section 9 We conclude this section by 
showing how the “follow through’’ method can be 
ised (1) to reduce substantially the calculations of 
appendix i. and (i) to proy ide rather vood upproxi- 


chosen ( 


Is 


abov © 


upp 


mations to the optimal (j, &)-pair, (14.23 Our 
irgument will show that J and j-: should be chosen 
to obey 

/> 105, k<281 14.26 
| 7 . 
tse of this information would have permitted 


significantly less work!’ in treating cases 8 to 12 
the most difficult cases) in appendix C.  Further- 


more, if we regard (14.26) as suggesting j—105, k 
“Sl as an approximation to an optimal choice, we 
nnd . Bs. Is less 


iu 


that for these values (¢ ( 


is the sum of three terms, the first of 


99 


percent above the value determined by 


We recall the definitions (14.1) to (14.3) of the three 
types of destinations. Our first assertion is that: 


863, then the 7th destination should not 
(14.27 


if 10°/ 
be a type 1 destination 


For, we found in the “simple”’ residue scheme of sec- 
tion 12 that if 10°f,<863 then mail to the 7th destina- 
tion should go into residue rather than receive its 
final sort in the primary. Since this applied after one 
loading complex operating at full capacity (and thus 
applies after the first complex of each of our primary 
subsystems here), there is even more reason for it to 
apply after a loading complex operating below capac- 
itv (such as the second complex in each of our primary 
subsystems). Thus it would be better to have the ith 
destination as type 3 than as type 1, and so (14.27) is 
proved. Of course, it might be still better to have the 
ith destination in type 2, and this is the next question 
to be considered. 


If the 7th destination is taken to be type 2. then 
(according to (9.14) and (9.17)) the cost associated 
with a letter to it would be denoted (fi, 2 The 


formula (9.18) for C(f;, n,) was derived assuming all 
primary loading complexes operating at capacity; 
this formula therefore provides a lower bound for 
C(f;,n,) in our actual system, so that the special case 
(9.25) of (9.18) vields 


S40) 
14.28) 


(0.52954 (if 10° 


1O7°/7 


If, on the other hand, the 7th destination is taken to 
be type 3, then (according to (9.11)) the appropriate 
cost is denoted C(f Except for a negligible correc- 
tion due to possible violation of (9.2), the formula 
9.12) still applies, vielding as in (9.24) 


134+ (14/10°f Gf 10°f,>60 14.29) 


Krom (14.28 and (14.29) we find that 


>) for 60<10°7,<455; 


i.e., if 10° 155, then the 7th destination should be 
taken as type 3 rather than type 2. 14.30 
From (14.30) we have 10°, >455, which for our 


particular mail distribution yields k<281. 


recall that in the discussion in section 9 
the change from (7%) opp=2 (Ni)oor= 1 
occurred for 10°f;=1,614; that for 10°/,>1,614 
it was better to sort a letter to the 7th destination 
directly after it enters the system, rather than either 
to put it into residue or to send it on past another 
loader complex. In our current situation this con- 
dition must be altered to take account of the fact 
that the second loading complex of each subsystem 
operates at only a fraction 1+ /,—F, of capacity. 


Next 


we 


16 > to 


see 


Is, 





ith destination should be of type 1 rather thar 
type 2 if 
> 864 1+ f f 


We find thus that for 


(1+-F',—F;,) 10°f,>1,614 (10°F 10°7,> 1,614, 


ih oA . ' articul: ea ave 
on the basis of input to the second loading complexes = *° that in particular we should hay 
only, the 7th destination should be of type 1 rather 1+F—s 
than type 2. 


1O°7 1.614 


which for our particular mail distribution Vields 
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follows that 


the 
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> 105 


Thus 


14.26 


is obtained 


Appendix A: Mail Distribution by Destinations 
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Appendix B: Analysis of Multiple Input 


We consider a multiple input system in which 


Systems 


\f=number of subsysten 


n—=number of loading complexes in each 


subsystem 


is 


(Bl 


(B2 


In order to specify the system completely, we must 


also specify, for 7 


D 


A 


i 


set of destinations whose letters are dropped 
in the row of bins after the 7th loading com- 


plex In each subsystem. 


(B3) 


i t fractior mail to the first i destinations so that 
F; i fi 

om 052 8H, 623 2s] O45 87. Oo 

62 O52 gH. O75 282 045 87. 654 

iS 052 86. 727 283 O45 87. HUY 

24 2 86. 779 2s4 045 87.744 

A) SH. SU 285 4 &7 su 
t “) a6. &7TO et 043 87. S32 

t 4) aH. 920 IRT 043 87. 87 

Jt a 86. O7U JAS 043 87. G18 

( “) 87. O20 JSU 043 87. UO] 

ul 87.079 200) 043 RA O04 

0) 87. 129 201 043 RA. O47 
2 OM §7. 179 ap 4 043 RA OUD 
; OM) 87. 220 203 0438 BA. 145 

{ Om) 87. 279 204 043 RA. 17¢ 
4) &7. 320 208 043 RA VIY 

1; 87. 376 2H 043 RA OAD 

47 &7. 423 297 043 RA BD) 

» 47 &7. 470 JUS 043 RA. 348 

47 Ss) | 204 043 mA BO] 

s 4 S87. 504 iM) (41 SS. 452 

A\ (s i] 
626 O13 0. 32 05. 042 
bol Ol2 5M W), 242 
Hib 7M oll 27 06, KIT 
olt ft 27 6. 792 
72t 4 Olt 2) 17,042 
T51t Ta 22 17. 2¢ 

vie iM (4s 200) W7. 467 
MOLTO’ OOS +A) U7. HHT 
F » SH OO7 07. S42 
R51 5 OT VS. O17 
x to GO Oy ] oR. 1902 
901 to 92 OF 1m) WS. $42 
126 to OO) Ow 125 US, 407 
U5i to ¥ Oo 12 Us. Suv 
176 to 1.001 OO > OS 717 
LOwlt ‘ OO4 200 OS. 17 
1 O51 to 1, 1 (ni4 JK) "4117 
1101 to 1.1m 003 1M) 09, 267 
115i t 1M 003 1™) WY. 417 
1,2 l,2 OOS 1s”) 090. 567 
Lg 1. 300 on2 100 09. 667 
1 2o1T 1 3 oo 100 Ww. 767 
13511 1K mr ”) V9. SO7 
1. 401 to 1,45 OO On 9. 917 
1, 451 to 1, 64M ool ORS 100. 000 


(D,, will necessarily consist of all the destinations.) 
The sets Dy, . . D, may overlap (for example, 
some destinations may be included both in J), and in 
D),), so that it is also necessary to deal with combina- 
tions of the D's such as D D, (the set of destina- 
tions included in J), but not in J),). If Gis any set 
of destinations (for example, Y might be J), or 
D,—D.), we use the notation 


fraction of the total mail which 
goes to destinations in &. (B4) 


We wish now to analyze the inputs to the various 
loading complexes, and also the input to each sub- 
system as a whole. Let 


101 











T,=ratio of the actual input to the 7th loading 
complex of each subsystem to its maximum 


possible input of 36.000 letters/hr s5 
The quantities /;, J, . . . , 7, (and thus the inputs 
36,000 J,, 36,000 Jo, . 36,000 7,) can be com 
puted one by one from the formulas 
I,=1 7 
I,=F(D 
I,=F(D.—D,)+I,F(D,), 
I,=F(D,—D.—D L.F(D,—D L,F (Ds), , 
Bb 
T,=F(D, D), .—D 
I,F(D Danas LD), 
+ J,-,F(D } 





To obtain a more compact notation for these equa- 
tions, we adopt the convention /,=0 and define the 
sets of destinations 


a | j<i-1); B7 


then (B6) can be rewritten as 


I,=>) F(Dt,)1 BS 
The input to each subsystem is given by 
(input/subsystem) =36,000 (/,+ J, as ae BY 


Appendix C: Proofs of Results Asserted in 
Section 14 


The following material familiarity 
with section 14, to which frequent reference is made 
We recall that AC had been split into two parts, 


presupposes 


AC=A,0+4.C, 
that a formula (14.22) had been derived for A,C, but 
that A.C could apparently be given by any one of 12 
possible formulas, leading to 12 possible cases re- 
quiring analysis. 

First we write down the conditions defining the 5 


cases (out of these 12) which can be treated most 
easily: 

Case I: y< 114, m >), no? 

Case 2: ) 114 m >) } ] 

Clase 3: y<114 mJ / j 

Case 4: z>114 m >) n>k 

Case 3: j>114, m >), gens 


proved logically im possible, 
14.16) it ean be deduced that Ta / 
ruling out cases | and 2 


) would \ ield 


These Cases will be 
14.13) and 
that m Nn, 
third conditions of cass 


From 
The first and 
using (14.16 


1/600 =fi4<f;<f 1 /600(2 + / / 
which is impossible sinc 
2+/ h,>2+/ 1>1+/ 1.724; 


thus case 3 is ruled out. The first and second eon. 


ditions common to cases 4 and a vield, using (14.13 
1/600 =f), > f 1 60006] / ] 
which is impossible since (using 7<k to deduce 
F< F,) we have 
l+F I | +f I 

thus cases 4 and 5 are ruled out 

Next we deseribe the technique to be used in 
handling some of the remaining seven cases. The 
parameters 7 and *& will be treated as continuous 
rather than integer-valued variables, so. that 


caleulus methods can be used in searching for the 
minimum of AC. As the two inde pe ndent variables. 


it is convenient to choose not 7 and k, but rather 
J and ul | i I C] 
Then & becomes a dependent variable, and from 
(C1) we have 
OF’../0) PP C2 
so that 


For our purposes, it is sufficiently accurate to replace 
14.13) and (14.16) by 


f 1 600 f 1 /600(1+ wu). C4 


and 7 hence m,n, F,, and Ff depend only on 


vu and not on J; vielding 


sO] 


Om/Oj= On 0) = OF, 0) = OF ,/ 07 =0 C5 
A necessary condition for a minimum of AC is 


O( Af O/ =, 


and if we introduce the svmbol J) solely as an 
abbreviation for the frequently-occurring quantity 

D=f,.(1+u)0(AC)/d), Cb 
then it follows that a necessary condition for a 


minimum of AC ts 
D—0 C7 


We will be able to eliminate a number of the remain- 


ing seven cases by sl! that they are imncom- 


patible with (C7 


OWilnY 


102 


nol 


ui 
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rom 
» SO 
and 
16), 


COl)- 


uce 


In 

} 
rhe 
WOUS 


that 


ntity 
(‘6 
for : 
C7 


main- 
conm- 


4 useful preliminary step 1s obtained by rewriting 


14.8) as 
AL s ] I u La) ] u 
rewriting 14.9) as 
Ab 760 / 28(7+k 5O(u I | -tw 


14.10) can be written 


as 


noting that 


LO°(u F,)/60 180. 


and then using (14.22) and (C1) to (C3) to obtain 


| +u)O(A,©)/ O07 17.139 < 10°f,f,(1 + 
147,(1--u 2.195 « 10°F, f, +392 } -f,). (C8 
The next 3 cases (out of the remaining 7) to be 


treated are 


(‘ase 6: ) <11 } MW ] / } 
( asé 72 | | } Mi / } } 
CLasé o-: yee 114 ms] / h 


These cases are logically possible, but we will prove 
that thev do not give rise to the minimum of AC, 


In case 6, AC =—A,C, so that by (C6) and (C8), 
D 17.139 & 10°F f.(1 +4 14/,(1-+-u 

2.195 x 1047; f,+-392(f fi ("9 
Since (C4) and the last condition of case 6 imply 

f,.(] u 1/600, 
e can deduce from (C9), using the fact 7 > f., that 
D 17.139 X 10°7,f,(1 +4 7TS4/, 

17.139 LO°/ 600 )7 7S84/ 0). 


so that in case 6, 0, contradicting (C7 
In case 7, AC’ is obtained by adding to A,C' the sum 


multiplied by 14,670) of (14.12 (14.15), and 
14.19 Using (C1) to (C6). we obtain 
D=—30.722 « 10°f f,(1 + u) — 14fju— .22 « 10°F, 

+ 2, 6607 + 2,674f,. 


since the third condition of case 7, together with 


C4), implies 


f.(] u 1/600, 


while the first condition (together with 7 <k) imphes 


i <9,<7 1/110, 
\)) ave 
D 30.722 * 10° 600) f, — 22 0007 ft, + 2,646 
2,074), = 2,646, —2,446f, — 22 ,000f f, >200/ 
22 0007} 2007.01 1107, 0. 


so that in case 7, D 
In case 8, Af 
product of (14.15 


we obtain 


0, contradicting (C7 
is Obtained by adding to A,C the 
by 14,670. Using (C1) to (C6), 


DD 


30.722 


LO°f f.(1 +a) + 13.363 > 


LO°F Af, 
+-378f,+2,674f,—14f,u. 


The third condition of case 8 (together with (C4)) 
implies 


f.(l-+u 1/600, 
so that (using the fact that /, </ 
D 30.722 10 600) 13.363 LO°7 53,0527 
/ (13.363 & 10°f,— 2,068). 
Thus DOO if 13.363 10°f,—2.068<0: ie, if 
/ 00155, or equivalently, if k>120. Hence, in the 


remainder of the discussion of case 8, we can ssume 
that &< 120, so that Fy < Fy, The second condition 
of case 8, together with (C4), yields. 


1/600 >f,1 + F,—F;, 
so that, using the fact F< F 


. we obtain 


1/600 >f,(.2662+ F.). C10 


By the first condition of case 8, f;>f\4, so that by 


(C10), 
1/600 >fii4(.2262 + F 

or FF’, <.7142, implying 7 S108, 
so by (C10) 1/600 >fios(.2262 + F 
O1 F< .6347 unplying I5738, 
so by (C10 | /600 >f73(.2262 + 
or I’, <.2949 implying 7213, 
so by (C10 1/600 >> f\3(.2262 + F 
or 

F< 0.1552, 


which is impossible since F, >0. Thus 
eliminated. 
There are now only four cases left. 


umiune 


case SS is 


We now ex- 


Case 9 yS114, m<j, gnsk 


It turns out that D=0 can occur in this case, but 
only in the subcase defined by 

99<7<114, I87<k< 281, C11) 
To prove this, we note first that in case 9, AC’ is 
obtained by adding to A,C the sum (multiplied by 
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14,670) of (14.15) and (14.18); using (C1) to (C6), so by (C18 
this vields 1/600 >7,(.12391 +4 
ol 
D 14.305 LO? fy f,u 14fu—30.942 LO°7, 7 f O020S8 
2,660 Jj 2,004 J implying 7>99. as asserted in (C1] 
To prove the remaining part of (C11), we hote t} 


. » gee » : . saa Q rao _ / 
By the second condition of cass together with the third condition of case 9, together with (C4 


(C4 (C19 vields 
which together with «<1 vields and since 99 and k <281. we have 
DD $4.305 ™ 10°/600)7 50.942 0 10°F f 2 646) f (2+ } Fe.) < 1/600. 
+2. 674f,—=2,646/7,—f,.(30.942 * 10°F 1710), Wielding 
OOO9UTG 
so that J >0 if 
s> that &>1S7 as asserted in (C11 
{<2 646/ 1710+30.942 ly Phe 3 cases not treated so far are 
bh (Clase 10 114 
or equivalently, if . 
( asN¢ pie 114 / / / 
fe <1/(1,169.4+ 1.780 f 
(Case 12 114 h 


Thus in the remainder of the discussion of case 9 
in all these cases, JJ=—0 can occur, and no reduct 


like that in case 9 upper ared possible We therefo 
resorted to numerical exploration for these thr 


we can assume that 


/ 1/(1,169.4 1.780 , 
cases and also for the range (C11 of case 9 
From the first condition of case 9 we have f,> fi. op ol hag ut 7 were used = — h the val 
, ‘Tl : () Ss PT pure $ ne co Ss why 
and combining this with the last inequality yields SS eee te panel: tand t] t-bahas. ye 
- ing value of AC’ was computed The results of t 
f 00045 calculations are given in the 1 lowing table. Ww! 
, ’ supplies the conclusions given in (14.23) and (14.24 
so that &<2S1 as asserted in (C11 Thus ; 
F< 
which together with the version i 
1/600 >f,.1+F I ; a sein 
of (C12) yields 1/600 >f,(1 +4 F'o,), O1 ” 
1/600 >f,(.12391+ F,). C13 — ' - 
Since 7>1, we have F,>F, so that by (C13 st 
1/600 >f,(.12391 + F 
or 
f 0102 implying J2 20, 
sO by bie 
| 600 >f,(.12391 + F 
or 
f << 00332 implying J 264 
SO by hk 
1/600 >f,(.12391 + F 
or 
f,<.00229 implying j291 Wasuineton, D.C, Paper 63B2-1! 
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New Method of Solution for Unretarded 
a Satellite Orbits 


John P. Vinti 
June 16, 1959 
(n axially symmetric solution of Laplace’s equation in oblate spheroidal coordinates is 
found, which may be used as the gravitational potential about an oblate planet. This 
potential, which makes the Hamilton-Jacobi equation for a satellite orbit separable, has an 
expansion in zonal harmonics in which the amplitudes of the zeroth and second harmonics 
can be adjusted to agree exactly with the values for any axially symmetric planet and a fourth 
irmonic which then agrees approximately with the latest value for that of the earth. 
Phe net result is therefore a reduction of the problem of satellite motion to quadratures, with 
ise of a potential field that is much closer to the empirically accepted one for the earth than 
any heretofore used as the starting point of a calculation. It may thus be possible to do the 
ravitational theory of a satellite orbit very accurately without use of perturbation theory 
Che method can take into account a first harmonic in the potential, in case observations are 
reduced to a center which does not coincide with the center of mass of the planet 
ction 1. Introduction 
eto! 
thre The purpose of this paper is to develop a method of solution of the central problem of 
de satellite astronomy, VIZ, the theory of the effect of the oblateness of a planet on the orbit of any 
non satellite, bound or unbound, but near enough to the planet so that the forces of other bodies 
if t! may be neglected A preliminary note has appeared in 1] 
whi To do so, let us consider the gravitational potential V outside an axially symmetric oblate 
a4 planet of mass 17. We call the axis of symmetry the polar axis; this will ordinarily be indistin- 
ruishable from the axis of rotation, except for entirelv negligible periodic effects. We also 
choose some point 0 on the polar axis as the “center” or origin of a coordinate system and let r 
be the corresponding position vector of a field point With 7 r|, the declination @ as the com- 
plement of the angle between r and the polar axis, @ the gravitational constant, ./ the mass of 
the earth, and w=GM, the potential V(r, @) can be expressed by means of an expansion in 
spherical harmonics 
V ur 1-4, Rir)"P <in 6) | (1) 
Here 2 is the equatorial radius, i.e., the radius of a section through O perpendicular to the polar 
‘3 axis. The notation J, for the coefficients was introduced by Merson and King-Hele [2]. 


For points on the polar axis (4 r/2) at the distance 7? from the center, the ratio of anv two 
| | . 


terms in (1) 1s equal to the ratio of the corresponding J’s 

If the planet has symmetry with respect to the equatorial plane, all the odd harmonics 
drop out; Le., Jj;=J3;=J;=. () In any case, even if no such symmetry exists, the first 
harmonic drops out, t.e., J, vanishes, if O is chosen at the center of mass The coefficients 
J, and J, are given by 


J,=25/3=(U,—1,)/MR, 2.1) 
mF 87/35. 2. a 


where the quantities J and 1) were introduced by Jeflreys [3], and where J, and J; are the 


vas supported by the U.S. Air 1 through the Office of Scientifie R the Air Research and Development Command 
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moments of inertia of the planet about polar and transverse axes through O. For the earth. 


with O taken at the center of mass, Jeffrevs [3] gives the values 


P 1637.5)10 Dd 10.7)10-°. 2 


j 
- 
jf 


while King-Hele and \lerson } derive from recent satellite observat 
J 1624.6+0.3)10~°, iz 6+1)10-°. oi 0] 510 39 


It should be noted that until very recently it has alwavs been assumed that the earth has 


symmetry with respect to its equatorial plane strictly speaking, that the geoid, or smoothed 





earth, has such svmmetry. That is, all the odd J’s have been assumed to vanish Recently, 
however, there has appeared some evidence [5] for an observable asvmmetry with respect to 
the equatorial plane. Such evidence, if confirmed, would have the following results If () 
is taken at the center of mass (, then J, would still vanish, but J., J would have small 
nonvanishing values. If, on the other hand, O is taken, either by error or by design, as non- 
coincident with (, but still on the polar axis, then we should have J, #0, with J, J. still 
small but nonvanishing. The only conceivable difficulty that might arise with this latter 
choice for the center has to do with the inertial forces acting on the satellite This difficulty 
is discussed and removed in the next section \ reason for deliberately choosing O as non- 


coincident with @ might be to represent any equatorial asymmetry largely by means of a first 
harmonic, and such a procedure might work if the third and higher odd harmonies then turn 
out to be negligible in comparison with the first 

In the present paper we first derive a solution V7, of Laplace’s equation, in oblate spheroidal 
coordinates, having axial symmetry and leading to exact separability of the Hamilton-Jacobi 
equation for the motion of a satellite It contains three adjustabl constants, bo, 6, and e, 


which may be so chosen that the expansion of V, in spherical harmonics is 
Vi=—yur 'fl—J.( R/r)*P,(sin 6) + J2( R/r)4{Py(sin 0) — J3( R/r)®P (sin 0 

J (Rr) Py (sin 6) + Jy J. R/r)8 Ps (sin 0) — J, Jz Rr)? Ps (sin 6 } 
where J, is given by (2.1) and where ./J; is the same as the J; in (1), vu 


Jim h 4] 


0 being the northerly displacement of the center of mass (' from 0 
(‘omparison of (4 with | shows that the fit is exact for the zeroth harmont 
and for the first and second harmonics It would be exact through all harmonies if the 


planetary distribution of mass happened to be such that 


J Ji, JJ J J 5 
J Side, Jg=d di, — Ly" Sad 5.2 
To obtain some idea of the adequacy of the potential } let us considel only the even 
harmonics. From (3) and (4) we may construct the following table of their amplitudes, with 
signs 
‘| ABLI 
H 
j 
| ” 
frey J (yn 
King-Hele and Merson [4 * 
j mi 
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The fit of the zeroth and second harmonics is necessarily exact If we fit to Jeffreys’ ” o 
the resulting fourth harmonic has an amplitude amounting to 49 percent of that given by 
Jeffrevs; the resulting error of (4) then amounts to 1.26 ppm, corresponding to an error of 1.26 
parts in LLOO for the oblateness potential If we fit to King-Hele and \Merson’s -/, the resulting 
fourth harmonic has an amplitude amounting to (86 16) percent of their value for the fourth 
harmonic, so that the agreement is within the observational error of them determination 
Even the sixth harmonic as fitted also lies within their estimates, which are, however, very 
hazy In anv event the resulting error of (4) then amounts to only 0.20 ppm, corresponding 
to an error ol only one part in 5500 for the oblateness correction 

The potential function that we find in this paper, leading to Hamilton-Jacobi separability, 
is thus much closer to the actual potential, at every point of space, than any heretofore used 
as the starting point of a calculation. If the results of King-Hele and Merson [4] are con- 
firmed, it may be possible to do orbit theory for unretarded satellites very accurately without 
perturbations of the solution of the Hamilton-Jacobi equation If not, there is still a rood 
chance of minimizing the use of perturbation theory, by distributing over the zeroth and 
second harmonics the secular effect of the fraction of the fourth harmonic that is not accounted 
fol 

2. Reference System 


As a reference system we first choose a right handed system of rectangular coordinates, 


with origin at the center O, with Oz pointing towards the planet’s north celestial pole and with 
Or pointing towards the planet’s vernal equinox. The directions of the axes are then fixed in 
an inertial system. For the computation of satellite orbits a slow precession of the axis of 
svmmetry about the celestial polar axis may be disregarded, if the precession is sufficiently 
slow and if the angle of the cone of precession amounts only to about 0.1’, as in the case of 
the earth [8, p. 202 Then Oz may be considered coincident with the axis of symmetry. 
The forees per unlit Mass acting on a satellite are then, relative to such a system and in 


the absence of drag 


F vravitational field of the planet VV. 6.1) 


F,— vector sum of the gravitational fields of the sun and of any moons that the planet 
mia \ have, plus the inertial forces on the satellite produced by acceleration of the center of 


mass of the planet by the sun and the moons, (6.2) 


F inertial force on the satellite produced by acceleration of the point O in the planet, 
relative to its center of mass C, 6. where 6 1s the vector from C to 0. 6.5) 


If in (1) we disregard all higher harmonics except the second and calculate —VV, we find 
that the ratio of the amplitude of the oblateness force Fy to the total gravitational force is of 


the order 4.5./ 1 200 Thus, for the earth 


F 5 cm/sec 71) 
Furthermore [6], 
F,| ~ 10-4 em/sec (7.2) 
and 
6 Ia /T’.)25 sin a 1.4)10~-°°6 see (7.3) 
for a period T, of precession ~433 days |3, p. 202] and for an angle @ of the cone of precession 
=0.1'". Thus 
F../} (] 4)10-"5 6 1n kilometers (7.4) 


so that / 2 Is utterly negligible ec ympared to F%, for any conceivable displacement 6 of O from C. 
In turn, £ amounts to less than 1/50,000 of the maximum oblateness force. 
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We are thus left with F, as the only appreciable force in the absence of drag, and we may 
choose our coordinate origin either at the center of mass C' or at any near point 0 on the polar 
axis that leaves mF essentially unchanged but simplifies the odd harmonic part of (1 


3. Oblate Spheroidal Coordinates 


We next introduce the oblate spheroidal coordinates S. @, and ?, le fined by 


r=c|(é+1)(1—7n*)|? cos @ Q< ¢ ©), S.1) 

y c| &-+-] | n~) |? sin @ l<nsl 8.2) 

CEN, S.5 

where ¢ is a distance to be fixed later For many of the following calculations, the book by 
Morse and Feshbach [7} is a very useful reference The surfaces é constant are oblate 


spheroids; sections perpendicular to the z-axis are circles and meridian sections are confocal 


ellipses of semiaxes ¢ ] and Ce, the focal separation being 2 The locus of the foe] is 
thus a circle of radius Cc in the plane : () The surfaces n constant are one-sheet hy pel boloids 
of revolution, asymptotic to cones with vertices at the™6rigin, and the surfaces with right 
ASCENSION o—constant are planes contaming the polar axis Note that 7 l along O-2 or ()- 
respectively and that for sufficiently large values of 7 ry? the coordinates & and y 
behave like =r/c and y=sin 6, where @ is the declination |7, p. 1292 The plane 0 consists 


of the regions 7—O0 and =O, outside and inside the foeal cirele, respective ly 


The element of are ds is given by 


ds* = hidt&’ + hidn- + hid¢’, Q 
where [7, p. 1292 
h c(é n é | 10.1 
h,=c(&+7°)*(1—1 10.2 
h ele | ] n 10.3 


4. Hamilton-Jacobi Equation 


Henceforth all quantities proportional to the mass of the satellite are taken per unit mass 
| pro] | 


of the latter. Bv (9) the kinetic energy of the satellite is then given by 


so that 


where the veneralized momenta are 


p h?=? —OS/Ok 13.1 
Pn =hin=OS/On 13.2 
p h 100) Os Od.” h cbead 


S(&, n, @) being the action function, 
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For the case of axial symmetry V=V‘(é, 7) and the Hamiltonian /7 becomes 


"£2+1 .,1 Is , 
I] ae" - D+ - ” pet f + V(E.n). (14) 
eT? - £°+- 7° e+riUd—-s 


The Hamilton-Jacobi equation is then 


. . . 9 ' 
t 1) (OS/0E l—*)(OS/On)°4 na OS /Od 207 (E°-+-*) V(E,n 2c*(&-+-n*)ay, 
. 7 ” 
(15) 
where a, IS the constant energy, 
To try to separate variables in (15), we express S as 
S=5,(E)+So(9) +83(9), (16) 
whereupon LdS becomes 
l l 
& LS l—7-)S +( )S 20*(£*+-*) V(E,y 20° (*+- n° )ay (17) 
l—n- ¢ l 
Here 
Si=ds,/dé, S3=dS/.dy S;=—dS,/do@ (18) 
From (17) it follows that S% is expressible as a function of — and » alone, so that 
S;=a,;—constant. 19) 
Further inspection of (17) shows that the variables € and » can be separated if and only 
if the potential 1(é, n) has the form 
VE, 9 E-+-n°)~"[ f(E) +-9(n)] (20) 


The question then arises: Is it possible to find a solution of Laplace’s equation, of the form (20), 
that will adequately represent the known facts about the earth’s potential? As indicated in 
the introduction. the answer is essentially ves. We accordingly devote the next section to 


finding all the possible solutions of Laplace’s equation that have this form. 


5. Solutions of Laplace’s Equation That Lead to Hamilton-Jacobi Separability 


For axial symmetry, Laplace’s equation in oblate spheroidal coordinates becomes [7, p. 


O _ ol O oO} ; 
bs t l—yn Q). (21) 
Of Of On On 


Insertion of (20) into (21) then shows that the functions /(£) and g() must satisfv the equation 


}?PQ0 


EE WF" (E) —2E(E +2) fh" (E) +21 +2) fh (& +07°(1—n") 9’’ (9) +20(9?—2)Q’ (9) —2(9?—2) g(n) 
+ M(£,n) 0), (22) 
where 
VE, 9 nk E“(7(n 23) 
Fk g Lyf?’ (é 2s (és 4,6: 24.1) 
Gm n> —1)q’’(n)+2nq'(n 2q(n). (24.2) 


By (22) a solution of (21) of the form (20) can exist only if 7 ¢,7) 18 separable. Now the neces- 
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sarv and sufficient conditions that \/(é, 7) be separable are that /(&) and G(») have the forms 


where A. B. and Care constants and where the signs and the factor 4 are chosen with a view to 


later convenience The conditions are obviously sufhicient: to show that they are also necessary. 
note that if \/(£, 7) is separable, there must exist functions @(£) and Y(y) such that 
nh é)— 2G(n)=V(n)—- (8 26 


On first placing n»=0 in (26) and then &=0, we find respectively 


b(t W(0 GO 7 7 
Win a0 F(O)\n 24.2 
Then. on mserting the eXpresslons 27 into (26 we find 
n° | t FO fl (r(n GO d(0 Wi) IR 
Placing either &=0 or n»—0 in (28) now shows that @(0)=—y(0 
so that 
e[F (i ko e--1G(n) —G(0 1/ 29 
where Cis a constant. bv a familiar argument Thus 
ke KO 1": ). | 
(i(n G0 1/7 0.2 
which have the general forms indicated in (25). so that the latter are also necessar\ 
(‘omparison of 24.1) and (25.1) and of (24.2) and (25.2) shows that the functions /(£&) ane 


g(m) must satisfy 


1)7’’ +-2&f’ —2f/—4C% te—0 LA 


By (23) and (25) it now follows that 


Wie. 9 LAig t/3y 2 


Insertion of (32) into (22) then shows that the functions 7(£) and g(nm) must also satisfy 
&*(¢ 1yf''(é 2E(; 2)f'(é 2 2) f(% tA 

7 Liq Y 2 2 a 2 2)\q t/} >? 
The usual argument about separability then shows that each side of +> s equal toa constant 


tA Thus f(g and q(n) must also satisfy 


We may avoid solving second-order differential equations, as follows: To eliminate /’ 


multiply (31.1) by & and subtract the resulting equation from (34.1). To eliminate g’’, multiply 
31.2) by 7 and subtract the resulting equation from (34.2 The resulting first-order equations 
are the linear ones 
Ef’—f £+-1)—c'¢ A+B)P?—K (35.1 
ny q n l i( 7" A B n K}. 35.2 


f(g Ki +-bt+Ce& A+ B—¢ AYE tans 36.1 
l+y 
q(n K— bn + Cn? +5 (A+ B—C+-K)n In 7 (36.2 
Z n 
where 4) and 4 are constants of integration. Insertion of (36.1) into (31.1) now shows that 
B= A. (37 
Then, with B= A, insertion of (36.2) into (31.2) simply yields a check 
Thus 
f(é K+ bEs- CE 2A—C'+A)E tans 38.1] 
, l . l+-y 
q(n k—b n--Cn--4 : 2A—( KC )n In 38.2 
y - ” 


It is easy to verify that these solutions for /(£) and g(m) satisfv both sets of second-order equa- 
tions, viz, (31) and (34), as they should Moreover, since f/(£) and g(m) must satisfy (35) and 
since the expressions (36) are the most general solutions of (35), it follows that there are no 
other solutions for /(£) and g(y 

From (38) it now follows that all axially symmetric solutions of Laplace’s equation that 


have the form (20) in oblate spheroidal coordinates are expressible as 


, l+-n 
Vien ( é ) LY hoe —b yn +b. ( QE tan~'&+n In ; 39 
; , bal In +2) 


where we have replaced 3(24—C+-A) by &. The Hamilton-Jacobi equation is separable, 
when the potential is a solution of Laplace’s equation, if and only if the potential has this 
form in oblate spheroidal] coordinates 

If as usual we take V to be zero at infinity, the constant C vanishes. We are thus left 
with three explicit adjustable constants, bo, 6;, and 6 It is at once clear, however, that 5 
must vanish if the solution is to represent the potential outside a planet, for the logarithmic 
term has a singularity evervwhere on the polar axis, for which l. 

We thus finally arrive at a solution of Laplace’s equation which makes the Hamilton- 
Jacobi equation separable and which contains two explicit adjustable constants, bp and 4y, 
and one implicit adjustable constant c, the latter appearing in the equations of transformation 


from oblate spheroidal to rectangular coordinates. This solution is 


V S n I bh, b ” by Re (& A} b, Im bS qn : +) 
The reader may readily verify that (£+-jn is a solution of the Laplace equation (21) and thus 
that (40) is a solution We here use j instead of ¢ for the imaginary unit, saving 7 for use In 


later papers as orbital inclination 
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6. Expansion of the Potential in Spherical Harmonics 


From (8) we obtain 


where 7 is the distance from the origin to the satellite and where 
sin @, $1.3 


6 being the declination of the satellite Then 


g yn is n 2jén—re 2)sre 12] 
where 
h - 12.3 
Then 
sTJ/T cl | hs h- 133 
Since we shall choose e< <r, it follows that {/ land (1—2hs+/ is then the generating 
function for the Legendre polynomials P (s Thus 
£ Lin l__ pp! eR h"P(s ep7! >> j Menyp—-MP (g) 14 


(£°+ n°) 'n Im(é+jn er fer" P, (sin 6 er P.(sin 6 er °P.(sin 6 15.2 


From (40) and (45) it then follows that 


V=boer '[1—e?r-?P. (sin 6) +e'r'*P,y(sin 6 . .|—byer™' [er P (sin 0) —e'r Ps (sin 6 
$4) 
7. Evaluation of the Adjustable Constants 
Comparison of (46) with (1) shows that to fit the zeroth harmonic we must choose 
ies mite 17.1 
that to fit the second harmonic we must then choose 
ct Ry, 17. 2) 
and that to fit the first harmonic we must choose 
be ulJ, =p. 17.3 


That is all the flexibility that (46) permits; the higher harmonics are then all fixed, as indicated 
in the relations (5), which follow from (46 


These choices lead to the form 1) for the potential, which has been amply discussed in 


the introduction. For R—6378. 388 km we find from (47. 2) and table 1 
e=210.7 km (Jeffreys), 18. | 
e=209.9 km (King-Hele and Merson 18. 2) 
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Although eqs (8) appear to place this constant ¢ in so central a position in the theory as to make 
numerical predictions depend strongly on its accuracy, it is not really so. Actually we shall 
see that a knowledge of the value of ¢ is needed only to find the deviations of the orbit from a 
conic section. 

From (40) and (47) we obtain as our working expression for V, in oblate spheroidal cO- 
ordinates 


V po~*(£-+- x? &+-nd/c). 19 
8. Formal Solution of the Hamilton-Jacobi Equation 


Insertion of (19) and (49) into (17) now results in the fully separated equation 
é Lys s 1) ~*ag— 2uck—2ay,c"é l— 17°)» | —*) ~'a2 +-2ydn + 2a,c*n*? =k. 50) 


By the usual argument, each side of (50) is equal to a constant &. 
At this point it Is necessary to discuss complete and incomplete orbits We define a 


complete orbit as the orbit that would be traversed by the satellite when the time ¢ is allowed 


torun from to ,on the assumption that the equations defining the orbit hold for all values 
of t. In the cases of a planetary satellite or of a meteor coming in from infinity with its orbit 
being bent around the planet without hitting it, the actual phvsical orbit Is complete. Any 
projectile fired from a gun and returning to earth, however, exemplifies an incomplete bound 
orbit. Any missile which escapes from the earth exemplifies an incomplete unbound orbit, as 
does any meteor with energy a, >0 which hits the earth. Note that the complete orbit may or 
may not intersect the planet If it does, we define it as that complete orbit which would be 


described by the satellite, in the absence of planetary matter, but in a gravitational field cor- 
responding to (49 For rc, this field can be represented by the multipole field (46). 
sy (50) it follows that 


k+-a l—7n*)S, —7*(1—n ag + 2y6n + 2a,c-n-. 51) 


Thus uf the plane n—O0 is reached bv the complete orbit 


In the case where 6 is chosen to be zero, the plane »=0 is the ordinary equatorial plane through 
the center of mass of the planet Thus, if this plane is reached by the complete orbit, (52) 
is satisfied. 


It is clear that (52) 1s only a necessary condition that the complete orbit reach the plane 


n—0. It is difficult to find a corresponding sufficient condition. By analogy with the Kepler 


problem ( potential), we may suspect that when 6=0 the equatorial plane will always be 
reached by the complete orbit of any bound satellite. If 640, however, the plane »=0 may 
not be ren hed if the orbit Is very close to the equatol This latter circumstance is illustrated 


in figure J 


| 

| 

j@— POLAR AXIS 
PROJECTION OF ORBIT 
ON PLANE THROUGH } ® 
POLAR Axi<s | 


iN 





PLANE 7=0O 





* 
i 
> c 
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The difficulties in finding a condition sufficient for »=0 to be reached are illustrated 
by the existence of cases where the complete orbit may not reach either the region 7»=0 
outside the focal disk or the region €=0 inside it. For example, if a particle is dropped from 
a point on the polar axis into a matter-free region satisfying (49), so that a,—0, 1, %=0, 
it follows (with 6 0) that the complete orbit Is a piece of the polar aXis, penetrating neither 
E—E=0 nor 7»=—0 

In the rest of the paper we exclude the possibility of orbits that lie only on the polar axis, 
thus ruling out, e.g., the case of a missile fired vertically from either pole. We then consider 
only those orbits whose completions either reach the plane n 0 or reach at least one venuine 
extremum n, for which 7=0, 740 

From (50) it now follows that 


S; t & ] ae t kh: 2uct 2ayt ya At a ] a3 1) 
So( l—yn |— ag he +-2udn +-2a,¢7n*)(1—7 a3.2 


where by (13), (16), and (18) the upper sign is to be used whenever the corresponding co- 


ordinate is Increasing and the lower sign whenever it 1s decreasing Also 


S: oO a 19) 


From (53) and (19) it then follows that 


AY Q2O 


hey 


| n ; ag i: 2u6n + 2a,c"n ’ dn+const 54) 


Before defining the lower limits &, and », it is well to note that eq (54) is somewhat deceptive 
| | 


In appearance, because of the presence of the “small” distance ¢ in terms that are clearly very 


important Let us now remove this difficulty, aus follows Replace the variable é by the 
variable p, where 
p Ce 5D 
By (41) it then follows that 
0 } ] n~)¢ yf). | 
pn rsin @ D6.2) 
be | r ( ( 9 
so that p=r on the physically realizable part of any complete orbit. On substituting (55 
into (54) and dropping the constant term, we obtain 
S QO p ( Ca I 2 p 2a po p ( do 
| ] cr / 2 ud7 rar | 7 dy 71 
where 
D Ce Zi 
If the complete orbit intersects the plane n—0, we choose the lower limit n, to be zero, 
corresponding to @ 0, in accord with the usual canonical treatment [8] of a bound orbit for the 
cause i" ul this uppears to he i vood chon efor an unbound orbit also lf tiie complete 
orbit does not Intersect the equatorial plane, we choose 7 nsan extremum ol al the complet 
orbit It will be a zero of the n-integrand satisfving 7 
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To choose a value for p,, we reason as follows. The p-integrand of (57) has the factor 


FF p C7 ae" (p--+-c*)f(p). (58. ] 


where 
I(p 2a,p* + 2up+k (58.2) 
will be almost equal to those of f(p), which always has two 


The appropriate real zeros of F' p) 
>, one zero of f(p) will be positive and the other 


in any physically realizable case. If a 
negative (approaching— © as a approaches zero). If a,< 0, both zeros of f(p) will be positive. 
Any other real zeros of F(p) will occur for inappropriate values of p, of the order ec. Thus for 
unbound orbits, we choose p, to be the largest positive zero of F(p), this being the only one 
accessible to complete unbound orbits. For bound orbits, we choose p, to be the next to the 


largest positive zero of F(p), equal approximately to the smaller zero of f(p). In either case 
p, is then the minimum value of p dy namically accessible to the complete orbit; the correspond- 
ing spheroid may or may not lie within the earth. (Any value of p less than p,; as just defined 


is dynamically inaccessible to the complete orbit in the field V(é,»), even if the planetary matter 
could be removed without change in V(é,9): to reach such values the orbit would have to 
pass through regions in which F(p) is negative and p; complex.) Of course p,>F and thus 
outside the earth for any ordinary satellite in a bound orbit. 

Whenever the plane n=O 1s intersected by the complete orbit, we have seen that k: Qs, 


so that we can then antroduce a real a. such that 


h as (59 1) 
where we can assume 
A () (59 2) 
without loss of generality. Then 
a5 a (59 >) 


If we put 6—0 and c=0 in (57) and then use p=r, n=sin 6, which would then follow from (56 


and also use (59.1), we find that (57) reduces IS} to the correct expression for the action func- 
tion for the case V jl Thus the terms in (57) that contain c or 6 give rise to deviations 
from a conic-sectional orbit 

From (13.3) and (19) it follows that the constant a, in (57 
The constant & has no name, but is related to the magnitude of the total angular 


is the z-component of angular 


momentum 
momentum, which is, however, not conserved in the non-central field produced by oblateness. 


To see the connection, note that whenever /&<s a”, we ean write (59.1) and (59.3). Then 


when 6 and ¢? are both allowed to approach zero in (57), a approaches the magnitude of the 


total angular momentum, which is then conserved This result follows most easily from the 


fact already mentioned, viz, that (57) then reduces to the correct action function for the case 


| ul 
9. Kinetic Equations of Motion 


Implicit equations for p, n, and @ as functions of f are now given by 


OS/Oa t+-B8 (60.1 
OS/O0k— B, 60.2 
OS/Oa,— 8 60.3 


where ¢ is the time and the 6’s are constants. Since the lower limit p,; in (57) is a function of 


a, k, and ay, the derivative of the p-integral with respect to any of them is the sum of an in- 
tegral and a term like | .(p,) Op, /Oa Here I,(p) is the p-integrand, which vanishes for py, 
so that this second term vanishes. If the complete orbit intersects the plane n=0, so that 7 
can be taken to be zero, on second term arises In a derivative of the n-integral; if not, m is to be 
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taken as an extremum value of 7 on the complete orbit and the second term again vanishes. 
since 7, is then a zero of the n-integrand. 
he kinetic equations are thus 
p l¢ a h 2up 209° p Cc il p ( ” a hk: 2udn rf | n adn f b;. 
6] | 
+3 ca I: =“p 2a p p ( dp 1 a } 2 udbn C7) ”) dn 3, 
‘) » 
Or Ca p ( era (k 2 up + 2a)p p ( i p 
a | ” a } 2udn Pa dn Db 1 ° 
If we restrict considerations to those cases where the plane n 0 is interseeted by the 


complete orbit, then kSa We can then substitute / a 57) and replace 
(60.2) by 
OS / Oa Db 2 
From (57 60.1] 60.5), am 62) there then follow the kinetie equations 
|’ | : : ) l” 
é ) ) ) 
i p°|caa a -Up-> <a)p p ( dpe n a Q5-> on 
a Pl e 
Jat ” f ) 6:3. ] 
a | l¢ a es <up a p p f lo cx a a 2 Za | ly 3 iff 
63.2 , 
OC a p ( Ca a 2up 2a p-)\p ( dp 
; { 
a | ” a a 2 On Jaye ” 17 ) 0.3.5 pr 
Future papers will be devoted to evaluating the integrals in these equations, soly ing them mT 
simultaneously for p, 7, and @ as functions of time, and deriving the astronomical results that su 
they imply, both for bound orbits and unbound orbits 
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Effect of Sudden Water Release on the Reservoir Free 
Outflow Hydrograph 
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The free outflow hydrograph is studied in the case of a sudden water release from a 
reservoir. The outflow hydrograph is called free when it is not affected by the tailwater 
vels Both the effect of the steep negative wave, created by sudden water release, and the 
ffect of the flow resistance are analyzed through the use of a fictitious inflow hydrograph. 
The water accelerated by the steep negative wave movement, whether or not the flow resist- 
ince is taken into account, gives this fictitious hydrograph. The procedures are given for 


celerities and t 
computation of wavefront heights and 


t he computation of the 


wave 


determination of the fictitious and total inflow hydrograph 


pV! imidal reservoirs of rectangular and par 
determination of the 


d 


thre reservoir free 


rive 


1. Introduction 
1.1. Definition of the Free Outflow Hydrograph 


If the outflow hydrograph from a reservoir is un- 
affected by the tailwater levels or if such influence 
Is neghgible, it will be called the free outflow hydro- 
graph. In that case both the outflow hvdrograph 
and the downstream propagation of the wave cre- 
ated by the sudden release of water can be analyzed 
separately, first the hydrograph and then the wave 
propagation 

The reservou free outflow hvdrograph in the case 
ol sudden water release depends On: (da) The area 
and shape of the outlet; (6) the type of flow (free 
surface flow or flow under pressure the shape 
characteristics of water body (storage function, Cross 
d¢) the inflow hydrograph; (e) the time 
interval in which a full opening of the outlet is 
f) the movement of negative waves, 
cally of steep negative waves; (q 


C) 
sections, etc 
created: espe- 
the resistance to 
flow along the reservoir 

The free outflow hvdrograph can be approximated 
by taking into consideration only the most influen- 
a) to (d \ny classical reser- 
vor routing procedure will suffice for the computa- 


lon 


tial of these factors. 


he analytical solutions ure possible Ik Cases 
when the outflow rating curve, the storage function 
the inflow hydrograph are given and can be 
approximated by simple expressions mathematically 
tractable [3] In that the epening time 1s 
as zero and the effects of steep negative 


nia 


case 
assumed 

ive movement and of flow resistance are neglected, 
assuming that at the very moment of fast opening 


he water in the reservoir was instantaneously ac- 
| Bi tw 
Ay I hir ( Ky , Ari 


outflow hydrograph in 


water velocities along the reservoir, for 
maximum outflow discharges, and for the 
The examples are given for the 
‘abolic cross sections \ general procedure for 
an approximate form is de- 


he new 
the 


celerated to give maximum discharge. This assump- 
tion is the source of the errors in the shape of the 
free hydrograph, if the steep negative wave is suffi- 
ciently high, as when the opening is large in compari- 
son with the reservoir cross sections. 


1.2. Purpose of the Study 


This study seeks the development of a simple 
procedure for the computation of free outflow hydro- 
graph, when the movement of the steep negative 
wave and the flow resistance are taken into considera- 
tion. The procedure must be suitable for use with 
the available data, to the accuracy demanded of the 
resulting hydrographs, and it must require a reason- 
able computation time. 


1.3. Scope of the Study 


This study will investigate the influence of the 
steep negative wave and of flow resistance on the 
shape of the free outflow hydrograph. The knowl- 
edge of the shape of reservoir (cross sections) and of 
its roughness coefficients is needed. Since the open- 
ing time is very short in comparison with the outflow 
time of the major bulk of water, it can be taken as 
nearly zero. 


1.4. Basic Approach to the Problem 


The shape of the free outflow hydrograph, com- 
puted by neglecting the effect of the steep negative 
wave and of flow resistance, is to be considered as the 
first approximation of the true outflow hydrograph. 
This first shape is sufficiently accurate for small nega- 
tive wave heights in many practical cases, because 
the accuracy of background data, the routing pro- 
cedure and the accuracy demanded do not justify 
more precise computation. The first approximate 
shape of the free outflow hydrograph is thus con- 
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ceived as a gradual emptying of the body of water, 
starting from already accelerated water [3]. If the 
neglected factors were taken into consideration and 
if more data about the reservoir were available, the 
subsequent approximations would be the 
true hydrograph 

Two approaches can be used 
of the subsequent approximations of the true outflow 
hvydrograph: (1) To develop a special procedure for 
computing the hvdrograph directly from. all v1ven 
data and (2) to the previous approximate 
hydrograph by introducing the effects of new factors 
In the second approach the effects are represented 
by a new fictitious inflow hydrograph. The varia- 
bles are time, the starting reservoir level, the charac- 
teristics of the steep negative wave, the length of 


closer to 


n the computation 


correct 


reservoir, the flow resistance, ete 

The new approximation of the true hydrograph 
can be computed from the first approximation by 
simple transformation using the fictitious funetion 
of inflow Next, the sums of the real inflow hvdro 
graph and fictitious inflow hydrograph can be made 
The new inflow hydrographs can be applied directly 
for the computation of the other approximations 
of the true hvdrograph This procedure will be used 
in this study 


1.5. Effect of the Steep Negative Wave and of Flow 
Resistance on the Free Outflow Hydrograph 


The outflow from a reservoir in the case of sudden 
water release does create a steep negative wave whith 
travels upstream along the reservoir. 

Every increase of outflow discharge also mcreases 
the head losses along the body of water 

When a steep negative wave is created, its 
ment is influenced by 

a) The change of the wave area AA 
shape (A, is the vertical projection of the wavefront 


Ao 


move- 
and of its 


and the change of the total cross section area 


of the area under the wave (.4-—A,). figure 1. which 
change the propagation celeritv: 4 the change ol 
the direction and the widt! ot reservou of the 
‘ 
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|.Steep negative wave & Y 
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2.Very gradually varied negative wove “> / 
h= front height of steep wove $ 
C 1C, ,C, =celerities of water body at rest VA 
H,, Ho,H,= mean depths “@ 
FIGURE 2 Two type eg ! : 
shape and the rough ess of its shores which CUuse 


the reflection and the refraction of th 
to the acceleration of Water 


and affect 
Wave; (¢ 
due to the movement of the wave, and the resistance 


the resistance 


of the accelerated water moving downstream betweer 


the wavefront and the outlet. which decrease thy 
volume and the kinetic energy of the accelerate 
wate The flow resistance decreases the outfloy 


discharges in the first part of the hvdrograph and 
increases them in its tail part 

As the propagation celerity of the negative wav 
gradually in the direction 
because ol the eradual decrease of mean depth Ina 
reservorr with sloping bottom, 
per unit of time thus correspondingly de- 

As the level at the outlet would be highe 
t than the level corresponding 
the next 
made lo 


decreases upstream 
the accelerated (ils- 
charge 
creases 
atany previous mome! 
to the arriving 
moment, a new level adaptation must be 
satisfy the continuity equation by a lowering of th 


accelerate discharge | 


level ni the outlet which Causes thre additio a 
acceleration of water the reservol The volum 
of water involved in this lowering ts. taken fro 


surface lavers of the bodv of water between tl 
outlet and the Steep hegative wave Therefore, a 
continuous propagatllol ol small nega waves 


changing very craduatlly take s place from the outlet 


upstream For each time unit a smal] drop in leve 


at the outlet is created which propagates upstream 
with the celeritv ¢ ,gll, where F/, is the meat 
depth (fig. 2 As H, near the steep negative wav 
is lower than the meat dept! of the lowest part 0! 


steep negative wave the small waves can not reac! 


\s the ait 


CONCEILN ed aus a Group of vers 


process is continuous, 


the steep wave 
wave movement can be 
small waves making a new long wave, very gradually 


fie. 2). with // 


changing decrease 


1.6. Approximations of the Free Outflow Hydrograph 


The combination of the movement of two types 


of negative Waves ofa very steep wave depression 
] 
and ofa very long waive corresponding to gradually 


unsteady flow with quite opposite char- 


varied ’ 
the complex reflections and refractions 


acteristics 





CUUSe 
of tl 


Water 


Stance 
‘tweet 
se th 
erated 
utflow 
h and * 


Wave 
CLIO! 
hina 
d dis. 
A cle- 
higher 
nding 
next 
ide lo 
ot tl 
tional 
olume 
fron 
nh wu 
pre, a 
vaves 
outlet 
level 
tream 
meal 
wa 
art ol 
reac! 
s. the 
very 
lually 


yraph 


types 
ss1on 

lually 
char- 
“tions 


of waves In a natural body of water, and the re- 
the flow make the exact computation 
of the free outflow hvydrograph very complicated. 
Therefore, the solutions must be obtained in ap- 
—_ ations of a different degree of accuracy. The 
first approximation neglected the movement of 
wative waves and the flow resistance (as in [3]); 
second approximation will include the ad lition 
effect; the third approximation 
of the additional effect of flow 
resistance These approximations are treated 
in the subsequent part of this paper Further 
approximations ¢an be considered when practically 


sistance to 


tne 
of the 
will take account 


steep wave 


LWOo 


fe as sible. 


9. Effect of the Steep Negative Wave on the 
Free Outflow Hydrograph 


4.1. 


As the breach of a dam or an accidental opening of 
some types of large gates are relatively fast processes 
. the negative wave created 
by sudden water release is a steep wave traveling up- 


Quite different is the case be the opening 


General Assumptions 


a matter of few seconds 


stream 
of outlets is very slow, the outflow rating curve 
changing with time as the opening becomes larger. 


(isa function of opening dimensions and of time 
This type of opening will not be studied here 

The steep negative wave undergoes three types of 
changes as it travels upstream: (a) Reflections and 
se of propagation celerity due 
and Cc 


refractions; h decrean 


to the change of cross section; loss of energy 
resistance 
and refraction of the steep negative 
wave are complex because of the factors which Cs- 
sentially continuously along the reservoir. 
The process 1s, somewhat less complex in 
channels of nearly I effects of the 
eflection and refraction of the wave 
will be assumed in the second 
as negheible for the 
The effect of the negative wave ends when 
ration celerity ( in water at rest 
flow velocity Vo already existent (C= V 
arrives at the shallow if water 
during the movement 
The effect of resistance will be nes elected In the see 
approximation, which is justified especially for 
re bodies of water and moderate heights of steep 
It will be considered in the third approxima- 

3 of this paper 

The main factor in this second 
the additional influence of the existence 
ment of ste ep negative Wave 


due to the 
The reflection 


change 
however, 
egular shape The 
steep negative 
and third approxima- 
tions valley reservoirs 

its prop- 
mean 
or when the 


equals the 


wave regions is in rest 


Waive 
al d 


Wives 
1 (sec 
! approximation. is 


and move- 


The flow through reservoirs or lakes, where the 
mean velocities are relatively low and the cross sec- 
tions change gradually, is a gradually varied flow, 


etl 


ler steady ol unsteady 

Before the analysis of the effect of the 
live wave on the outflow hydrograph is begun, 
roblems will be studied: (1) The propagation of the 
wave along the water body the steep wave 
height and the outflow shortly after the 
sudden water release. 


steep hega- 
two 


and (2 


discharge 





Sch 


oO} the 
due to diffe rent 


transformation on 
celerities 


FIGURE 3 ematic 


the front s 


re prese ntation 


ha pe of sleep negative wave 


cave elements 


The shape of a steep negative wave moving along 
the body of water is not stable (fig. 3). The reason 
is that the parts of wavefront travel with different 
celerities. The steepness of the wave decreases with 
the distance traveled, and the final shape will be ob- 
tained when the slope of the wavefront is essentially 
equal to the slope of the bottom. Because of these 
differences in celerities of different points of the wave, 
it is necessary always to define C (as that of the high- 
‘st point, of the lowest point, of the center of gravity, 
of the mean height of wave, ete.). The approxima- 
tions in this study deal with the schematic negative 
wave, considered as a discontinuous steep wave. 


2.2. Mean Celerity of the Steep Negative Wave in 
a Cross Section 


The steep negative wave arrives at a cross section 
of the body of water, X distance from the outlet 
at time ¢ from the opening time, with wavefront 
height A. The celerity of the wave in water at rest 
is ( at XV. The mean velocity of water because of 
inflow is V. The eross section characteristics are 
fig. 1 Area .1 before the wave arrival; area A, 
of the wavefront height; distance f of the center 
of gravity CG of area A, from the new level (or from 
the starting level / h f area 1 A,, after the 
wave has passed. The effect of the bottom slope 
is indirectly imtroduced by the change of A, A), 


fo f;, and uA A 

Two equations, continuity and momentum, give 
the celerity of the steep negative wave and the 
increase of the mean water velocity in a cross 
section of reservoir, due to the acceleration caused 
by the wave movement, as functions of section 
characteristics [1, 2]. For a small time interval dt 


the celerity can be considered constant in the vi- 
cinityv of the section at .Y, regardless of the 
bottom slope which has little influence on the cross 
section characteristics for short path length (fig. 4 
The steep negative wave travels upstream a 
distance (C’—V)dt during the time dt (fig. 4), and 
the water travels downstream the distance (V+ L’)dt 
An observer, traveling at the section YY with 
the velocity Vo downstream, after time dt, will be 
at the section Y— Vdt, or a distance Vdt downstream. 
To him the distance traveled by the wave in the 


CTOSS 
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Volume |-2-3-4-| =A At 


i x 4 


= Volume 456-7-4=(A-A,)U,at X 
A, Cy= (A-A,) U, 


Fiat RI 1. Scheme for the der ation ot cor nuily eq satior n 


case of movement of steep negat 


upstream direction is Cdt and the distance traveled 
by the water downstream is (df 

The continuity equation gives for the observer 
(or also in case V=0 


AC dt A—A,)l dt, 


so that generally for any 


celerity. of the 
\"is the prop 


where (' is the 
negative wave in water at rest and ¢ 


propagation steep 


agation celeritvy in moving water with velocity 
V (V: (' Equation | can be obtained also by 
equating the total water discharge entering” the 
small volume cidr, to the discharge leaving that 


volume 


For the momentum equation, it is assumed that 


the acceleration of water from the velocity V to the 
velocity V+U’ of the volume of water (1—.1,)Cd?. 
in time dt and at the section .Y figs ft and 5) is 


made without loss of energy. 

The volume element (A—.1,)Cd?, which ts accele- 
rated from V to the velocity +l times 
t and ¢+dt, has a momentum (4li—.A1,)Cl dt and a 
rate of change of momentum (.1—.1,)C7 The 
force F which creates this change is the difference 
of pressure on the areas «1 and w1—<1,, (fig. 5 
F=gh(A—A gfA, with f already defined 

This force is equal to the rate of change of momen 
tum, so that 


between 


(7 Oh at . 
gh+9I a—F 


The celerity (at rest 


after using eq (1), Is 


any section im water at 


C=ygh(a—| gf, 2 


120 





— 





> »'F 
> A-A 

ae 
nw - na 
pg és 

- AF B. Nm 

x A= B, Hm 
A= area |-23451 
A,sarea 123-54 
A-A, #area 3453 
| Ic;U RI 5 Scheme fo / mt OF Mmome? ” equatltor 1 
where a=.1/.1 The ratio a>1, and theoretical] 
can be between the limits 1 and infinity if 
mcreases, a ] decreases faster than / and 


so that the mean celerity of the wavefront decreases 
with the increase of / With the mean depths 0 
areas ol and J1,, // AUB, and h A,/B, (fig. 5 
with #4, the width at the surface level, and putting 
/ h j, eq 2 


hee OMCs 
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a, 


and 7f, can be 
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/; him Is heariy 


When 


i comparison with // for 


hecrectler 
a4 


small values of A) the 
( \gH : } 


which 
of small surface disturban 
the mean depth I] 

The imerease of the mean 


| and 2 


Is the know! eXpresston iol the 
ces in shallow 


propagauio 
water, wit! 
water velocity is. there- 


lore, [rom eqs 


or trom eqs 


which approaches ZeTO W hi 1 h approaches ZCTO 
The discharge P AA A IS the accelerate 

movement of the 

isidered as the fietitious inflow 


discharge by the steep necatve 
wave, which will ay Col 
in the reservou 
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characteristics and 


ety 
functions of seco 


wavefront A, bot 


CTOSS 


and j are 


ot height of 


( and are functions of those parameters The 
cross section characteristics are functions of the water 
level For a given section, level, and height / th 
values ¢ and l” are Ung ely determined 


Wally 
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hot! 
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93. Computation of the Initial Front Height of Steep 
‘Negative Wave and of the Maximum Outflow 


Discharge 


If an opening in a reservoir is created with constant 
area and shape, its hydraulic characteristics are 
defined by the outflow rating curve, generally of the 
type @ bk’, where 6 and r depend on the area, shape 
and position of the opening, as well as on the type 
of flow, with & the total energy head above the level 
of the outflow zero (fig. 1). The opening height, as 
the difference between the level at the moment of 
opening and the level of outflow zero, is ho. It is 
the sum of two heights: A, the front height of the 
steep negative wave and ky, the difference of the 
level just upstream from the outlet (where the 
depression due to the outflow is negligible) and the 
level of outflow zero, so that hyo=h-+-k The energy 


| ead 1S 


where Vis the velocity near the opening before it is 
created and ais the kinetic energy correction factor 
For the small time interval, Af, after the opening, 
water flowing through the opening is the sum of 
e inflow discharge P=AV and the 
discharge due to the moving steep negative wave 
The interval Af must be sufficiently long for the 
negative Wave to develop but not so long as to 
nerease the outflow by the further significant drop 
water level of the stretch of reservoir between the 
outlet and the steep negative wavefront 
The continuity equation between supply and out- 


flow gives 


accelerated 


Q,=—=AV+A,(¢ \’ Ak’, S 


‘ 


hich ean be written as [4, p. 213, without last term] 


Q,.—=AV-+A, (¢ V 


t 


For small value of velocity head in comparison 
th ko, eq (9) can be expressed in the approximate 


T ( second term of 
} 


t, p. 213) in form of (V3/2q 
1 takes here the form 


h k | l | Ob 
| 


d should be added on the right side of eq (9 
Ot neglected for outflows in 
with the first term The celerity By 
“ , and velocity 

For given h : h. and r (values or functions depend- 
nt only on the outlet characteristics), and given 
the magnitudes A,, C, U, and ky=ho—h 


, Can 
reservou comparison 
is viven by eq 
or (3 


by eq (5) or (6). 


V and a, 
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are functions of A and of the cross section character- 
istics near the outflow, at distance AY =—(C— V)At. 
When the section at distance AY from the 
outlet is selected, eq (9) enables the computation 
of the value A. As the discharges around the maxi- 
mum discharge change slowly (fig. 1), the selection 
of A.V is easier. It is, therefore, sufficient to select 
an average section some hundred meters 
upstream from the outlet. For given h, the values 
ky and UL’, or A; and C can be obtained, so that eq 
(9) also gives (Yo, the discharge outflow at the time 
At (which is the maximum outflow if AX is well 
selected), 


Cross 


CTOSS 


The computation of h is to be made by successive 
approximations until eq (9) is satisfied for the com- 
puted fh. Numerical or graphical procedures are 
feasible for that computation. 

A simple graphical procedure for the determination 
of h, Y@ and downstream tailwater level 47, is shown 
in figure 6. Three curves as the functions of h and 
independent of the height ho of the opening are first 
plotted 


Q,=A,(C—V)+4AYV, 


and 


Q,=f(,), 


where Q, is the rating curve of the tailwater channel 
Kor the given opening height h,. the fourth curve 
().=bk’ is plotted. For some types of outlets, the 
coefficient 6 is a function of kp/(77f>—hy) and of son.e 
other parameters, so that approximate values must 


be first introduced. For convenience the velocity 
head Ak, is plotted with the same scale as ho, h, k, 
and 4 For an assumed value é or ky (point 1) the 


discharge () A, (C—U 
point 3 


AV is obtained (curve 1, 
, as well as Aky (curve 3, point 2). By adding 





Ak, to k, (point 1) the value & is obtained, and 
(o=bk" (curve 2, point 5). The adjustments have 
to be made for A until Q,=Q.. If Q;=—Qo, the value 
of Af is obtained. From that value the discharge 
‘ yy = 
, * ke 
nVrvuU 
ik 7; 
eda — ‘ 
ai eon depth upstreor 
f the outlet 
4 
’ 4 
| 
t ZLa_o ft _ 
FIGURE 6 & aph ( procedure for the determination of front 


he ght of sleep negali e wane of marimum d scharge, and oft 


j j 


maxrimun a waler leve 


lischarge plus inflow P 








(), is determined (point 6). For a given QJ» (curve 4, 
point 7) the maximum downstream level //, (point 8) 
is obtained. When //,<8(//,—hy), the outflow ts 
uninfluenced by tailwater level. The coefficient 
8>1 is a correction factor which takes into account 
that, due to the effect of kinetic energy of outflow 
water, the tailwater levels somewhat higher than 
(F1,—ho) do not influence the outflow. 


2.4. Relation of the Height of Steep Negative Waves 
to the Opening Height of a Rectangular Outlet 


For the given opening, the values 4 and Q») depend 
on the shape of three curves, 1, 2, and 3 of figure 6 
The accuracy of curve 2 (V=6k") depends on the 
accuracy of the dimensionless outflow coefficient m 
which enters into the expression for the coefhicient h 
with water at rest, of 

BH 
bottom 
assumed 


The case of a reservoir, 
approximately rectangular cross section A 
near the outlet will be treated here The 
of small length of reservoir A.Y can be 
horizontal, especially under natural conditions with 
deposited sediments. The simple case of a rectangu 
lar sharp-edged opening will be analyzed here as an 
example, with area Bho, where Ao is much smaller 
than /7/,, (fig. 6), and B< AB). Equation (9) will be 


used, and according to [4] and [5], with different 
symbols used here, 
Q) i) } bh B 


.¢ t 


Mi 5 > > “s, 
Bhy gh B voh ? b Vd b 


B _ 
m PelG7, j R,W) 10) 


with kh h h and R and 1 the Reynolds and 
Weber numbers, respectively 

In this expression Bb is taken in the denominator 
instead of the usual 2, so that the right side of eq 
10) is multiplied by 4/4, apart from the contraction 
effect dependent on B/B,. m order to facilitate the 
changing width 


solution of eq y regardless ol the 


B of the opening 

Putting h Il sand h }] P, 
of the wave height A and the opening height / 
be related to s For the rectangulatl sharp-edged 


Opening 


the ratio Ach p 
will 


B, (Bs 
B | By | 
where mis B,/ B times higher than the usual coefficete 
m Other functions of m 
have to be used, depending on the outlet 

The coefeient m can be thus expressed asa lune 


for other shapes of outlet 


tion of dimensionless numbers, so that for each outlet 
and its characteristics the special relation of eq (11 
or of other tvpe can be obtained. For given b, B 
HT,, and ho, the coefficient m can be expressed as a 
function of A, or of p AH 

Equation (9) under the above assumptions and in 
terms of the dimensionless numbers introduced here 


can be written as 


If the function of eq (11) is introduced in eq (12 
the relationship p=f/(s) ts obtained For the practi 
cal solution it is suitable to take m as au paramete 
and to develop a family of curves for eq (12). Figyy, 
7 gives this family of p=f(s) for m=0.25 0.50, 0.73 
1.00, 1.25, and 1.50, and figure 8 gives p/s=fls) fo 

, 0.50, 0.75, and 1.00. The families of curyes 
/ (“) for given values s are more 


soth are given in figures 7 and 8 


1.05 was take 


ceonventent to use 


as dashed lines. For these curves a 
as the most applicable constant valu 
When eq 11) is determined for a particular easy 
of rectangular opening with s=hy // 
intersection of the curve of eq (11) and the ecury 
m fh ip) orn f(hih f(p/s) trom figure 7 
or S tor that Ss, ZlVeS thre solution either ior p=h H 


or pP N h h 


known, th 





It the velocity V 1th the body Ol water near th 
outlet, before the sudden release of additional water 
can not be neglected in comparison with the new 
velocity | {’), it must be introduced in dimen- 
sionless form in eq 12 In that case, for 
ia ygll l,/¢ . with ¢ Voll and | P BAH, 

Pi 
as ; J 
* r 
f 
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For each given value of w a family of curves 
fh i or ” fihoh can be developed 
Equations (12) and (13) in parametne form can 
developed ior any regular-shaped outlets rec 
gular, triangular, semicircular, circular, trape 
oidal, parabolic, semielliptic, et With s, or with 
and w, as parameters, the families of curves 
fh h or WwW / h I] can be developed once 
all cases With these families available, the 
etermined real funetion m fihsh lor a particular 
itlet of given regular shape will be sufficient fon 
determination of the height 4 of the steep nega 
re wave, and of the maximum discharge of outflow 
drograph QQ) by using eq (8 The nomographs 
also be developed for the relations in eqs 12 
1 
If mis selected to be only the function ol h, then 
l2) would give a stable relationship between p 
or eq (13) would give a family of curves with 
is parametet! 
2.5. Procedure for the Computation of a Free 
oe Outflow Hydrograph 
6 In the case of no inflow imto the ress rvoi, the 
WH.) celerity C depends on the wave height 4 and on the 


ea depth H, of reservoir at the outlet before the sudden 


ater release The celerity is, therefore, the veneral 


function 
h H 


"WT Hi, (14) 


N/L is the relative distance of a cross 
section the upstream end of reservoir, Y is 
the from the outlet, ZL is the length of 
reservoir, /7 is any possible reservoir depth at the 
outlet, //, is the maximum operational depth, A 
the height of the steep negative wave. 

For a given relative parameter //7///), a family of 
curves ¢ f(r, h/ Hy), can be plotted. 

For a large reservoir, Hy and / will be assumes 


where r=] 
from 
distance 


i ibs 
the initial or boundary conditions. For a given 
cross section « the values A, A,, a, f, C, and UU’ 
entering into eqs (2), (3), (5), and (6 are also 


viven 
The reneral procedure for the second approxima- 


tion of taking into consideration the effect of the 
movement ol steep negative wave in the case of 
valley reservoirs ts: 


1) To compute the wave height / for given water 
depth /7 at the outlet, and to determine the maximum 
outflow discharge Q, as is shown in subsections 2.3 
and 2.4, to (13), figures 6 to 8 

2) To verify that the outflow hydrograph is free, 
or that the tailwater level is lower than the maximum 
tailwater level which not affect the outflow 

fig. 6 


eqs (7 


does 
3) To compute the celerity C along the body of 
water out of given cross section characteristics (A, 
A, and f and (3 

1) To compute the velocity V along the body of 
water for the inflow P? during the wave movement, 
if there is an inflow; 

5) To determine the difference (C—V 
body of water; 

6) To compute the accelerated discharges by 
steep negative wave movement along the body of 
water 2? (A, with inflow zero, or P,.=(C—V)A 
with inflow P: 


4 To determine 


. eqs 


along the 


the time, when the wave arrives 


at the section _Y and when the acceleration takes 
| by f¢ [ ds for 1-2 

ace there, = yr or — 

pia ere \ 7 / | ol I j i 
L length of reservoir by 

“lL dia 
f 5 
L| ny 15 


with /, the total time of wave travel from the outlet 
the end reservoir (from s=1 
I () 

\ 


to upstream of to 
To find the relation between the accelerated 
and the time acceleration, ~P? f(t 
For given « the pair of values ?, and ¢t gives P f(t) 
The integration of P, over time from t=0 
t—t the volume S, of the reservoir laver h, 
determined simply from the storage function S=f(/7), 
figure 10. which ean check the function 


r f(t 


discharge of 


to 


CiIVeSs 


serve to 
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If the inflow during time f) cannot be taken as 
constant (i.e., the natural inflow hvdrograph or the 
release of water from upstream reservoirs are very 
unsteady , the velocity V at a cross section « is also 
dependent upon the time as are (’ V, Ay and r 
The inflow can be either a steep or a gradual wave and 
the determination of ?,,—/(t) becomes more complex. 

If the inflow is in form of a steep positive or nega- 
tive wave, the procedure is similar to the above given 
movement of steep negative wave until two waves 
traveling in opposite directions meet The resulting 
waves will then move in both directions. 

In the case where the inflow hvdrograph P=/(t) is 
gradually varied, and where the steep negative wave 
moves faster upstream than the other negative waves 
created downstream of the steep wave, a resulting 
new inflow hydrograph P?,,=/(t), composed of P=f(t 
and the accelerated discharges by steep negative 
wave movement, can be determined either by exact 
or approximate methods. The position of the wave 
P=f(t) along the reservoir has to be determined 
Therefore, the velocity V for given places and the 
times at which the steep negative wave arrives at 
those places must be computed in order to obtain the 


(x ' h if ), 
"iE 
with a new variable, ¢ 

The time ¢, of wave travel from one end of the 
reservoir to the other is the factor which will deter- 
mine whether or not the inflow hvdrograph can be 
replaced by a constant inflow Py. Since ¢) for the 
majority of valley reservoirs is measured in hours, 


celerity C' V The celerity is now ¢ 


constant inflow can often be assumed during the 
movement of steep negative wave 
Due to the neglected factors resistance to the 


flow, reflection and refraction of the wavefront and 
transformation the wavefront due to different 
celerities of its parts the computed function ? At 
is an approximation of the true relationship 

The procedure for the second approximation of the 
free outflow hydrograph would involve the following 
values and functions: 


(1) The head k—/ “(2 


24 

outflow rating curve Q=6 k’; (3) the characteristics 
of cross sections along the reservoir, and (4) the total 
inflow function ?,,—f(t), which gives for each time ¢ 
the discharge ?,—?+/P,, and for the corresponding 
r, also, the place of inflow P? The use of the two 
basic de Saint Venant equations for gradually varied 
unsteady flow will enable the computation ol the free 
outflow hydrograph in cases where the continuity 
equation alone can not give reliable results 

In vallev reservoirs, where the height of the nega- 


of 


a( V+ 


outflow the 


tive wave is not too large, the use of the continuity 
equation alone, instead of both continuity and mo- 
mentum, gives an approximate free outflow hydro- 
graph, sufficiently accurate in the light of background 
data, assumptions and the needed accuracy of results 

The higher the steep negative wave and the inflow 
P, the greater ts the velocity (V+ 0"), so that the use 
of the continuity equation alone will give less accurate 


results. Since the effect of neglected factors in steep 


negative wave movement is the highest in shalloy 
waters (upstream part of reservoirs) and lowest i 
the deep waters lower part ol reservoirs and Sines 
P, is low in the first and high in the second regio) 
the use of only the 
justified In that case the integration of the 
differential equation COntIMmUEY equation 
tvpe P—Q=dS/dt, with P replaced i 
hvdrograph P?,=— P+P 
graphical, or any other appropriate procedure to ob- 
tain the free outflow hvdrograph A i(f 

Two simple cases will be taken as examples: (a) 4 
pvramidal reservoir with rectangular sections B 
Hy), with the bottom slope //,/L, where H, is jts 
depth at outlet and Z reservoir length. and with B 
the width of reservoir at the outlet: (b) a paraboliy 
section of the tvpe ib b \ I] HT, where 3 and I] ary 
width and depth respectively of reservoir at the out- 
let, Band 7/7 any width and depth, with the botton 
slope //,/L, and Bb, and /7, the surface width an 


water depth at a cross section VY distance from outlet 


continuity equation Is general} 
Slorag; 
Of th 
the hew infloy 
can be made by a numerics 


2.6. Rectangular Cross Sections 


Figure 9 gives the simple shape of the reservoir an 
the resulting functions along it for the in flow P () 
The steep negative wave has the front height A wit! 


p h/H,. anda | N/L), so that BB H/H,=2 
The total area of any cross section Is y Bo 

and the wavefront area A thB,: f=h/2. so that 

{ 

| , 

_ “Se 

| 

4 t 2 

} 

® 4 
PiGuREY (' eee 

ci de j ng ; roll 

nagiite ) i] 7 the 
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, 9 
from eq. \- 


alloy j 
est iy i a A 16 
sine \ 2 
egion 
erally | with Gy=v glo, the maximum celerity of the upper- 
ora, most part of steep negative Wave hear the outlet 
f the In the case of constant inflow Py and V=P)/A, 
infloy 
Pr cal ‘ » V - 
to ob. ( | t \ Ir- E ; ‘ 
| 2 
MA vith Vy=Po/ Bolly the velocity at the outlet before 
(B the sudden water release. The new velocity along 
AS Its the body of water after the wave passes Is, according 
ith B - 
aboli¢ 0 OF 9s V 
Ty ar V+ 2+ Ka 2 5 Is 
¢ Out- : r\ = 
ra For small values of 4 eq 16), | 17), and (18 have no 
utlet practical meaning, because ( and l ure very high, 
and the resistance plays the major role, so that they 
should be used for values slightly higher than p 
The accelerated discharge plus the inflow P, is now 
: P.=P)+Po, with 
Ir and 
P=( ; 
with P? WBA Cys ! 4 ) 19 
T,=1 \ 2 3 
BH 


that The travel time from the outlet to a cross section 
Y distance, in the case V0, Is 


x #4 “dz 2L 
t L| l—a 20 
; ( Jx Cc ( 

Ifv—0, the travel time along the reservoir for V,—0 1s 

») 

26 
( 
, / . ane , 

and r=1 1—T, where T=t/t With the in- 
ease of ?) the time fy) Increases and the discharge 


though r r 


the ¢—funetion for 


decreases, Increases, 


The Integration of V,+0 m 


lorm 


s best made numerically. 

The functions C0’, and P’,, are given in logarithmic 
scales in figure 9 for the following values: H)—100m; 
By=300 m;h=10 m; Cand’ in m/sec and Q in m‘/see 


2.7. Parabolic Cross Section 


For given h and Hy, and two cases of the inflow: 
a 0 and P=/, figure 10 gives the following fune- 
th ree “ons depending on #=1—CYX/L) for the parabolic 


dy HH, where d= Bo/ Fh; 


ross sections of the type B 


or C—V);U (or U+V); V=P2/A; 
da 
ji (’ ps a (or P,), accelerated 


discharge P,, or P, Po+<A, (C—V); and finally P, 
or ?,) as a function of time ¢, (or f.). 


A: A;; a=A/A,; ( 


f ( “(ort 


2.8. Shape of the Hydrograph in the First Moments 
After Release 


A few moments after the release, at the time in- 
terval Af, a full steep negative wave of height A is 
developed (fig. 1). At that moment the outflow is 
maximum. ‘The approximate solution for the shape 
of the rising part of the hydrograph after the sudden 
water release can be obtained as a straight line for 
an ideal fluid by using the two-dimensional hydro- 
dynamic equations of motion expressed in Lagran- 
cian form [6]. 

If the straight line rise is assumed for the total 
rise, the procedure for the determination of the slope 
of the hydrograph is simple. The total outflow of 
water from the time of opening to the time Af is ap- 
proximately the volume of water confined between 
the outflow section and the section at AX—(C— V)At, 
where (C—V) is the wave celerity at the beginning, 
and between the starting depth of water //7 at the 
outlet and the new depth I1—h. This volume is AS. 
If it is taken that the maximum discharge (J) takes 
place at time Af,< Af, instead of at Aft, the continuity 
equation gives '(At+-At,)Qo=AS, and Af, can be ob- 
tained (fig. 1 The mean slope of the straight line 
for the rising part of the hydrograph is then Qo/At. 
For small A, the values AS and AT are very small, 
and the rise of hydrograph can be taken as vertical. 

Figure 11 gives an example of the difference be- 
tween the outflow hydrograph of the first approxi- 
mation, curve 1 (emptying of the reservoir with the 
negative wave movement neglected) and the outflow 
hydrograph of the second approximation, curve 3 
(with the steep negative wave taken into considera- 
tion, but with the flow resistance neglected) for the 
Savage River Reservoir near Cumberland, Md. A 
rectangular opening (250 m by 18.5 m) with sudden 
water release in the case of a full reservoir is 
assumed Curve 2 gives the fictitious inflow 
hvdrograph computed by the method given above. 
Though the total volume of the two hydrographs, 1 
and 3, is the same, the peak discharge in the second 
approximation is smaller for about one-third 


3. Effect of Flow Resistance on the Free 
Outflow Hydrograph 


3.1. Two Aspects of the Flow Resistance 


The effect of the flow resistance on the free out- 
flow hvdrograph depends on the new velocities along 
the reservoir. The greater the wavefront height h, 
the less the remaining flow area A—A, after the 
steep negative wave has passed, and the greater the 
inflow and the roughness coefficient in the reservoir, 
the greater will be the flow resistance effect 
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Two aspects of flow resistance can be taken into 3.2. Effect of the Flow Resistance on the Steep Nega 
consideration 1) Head losses along the reservoir, tive Wave 
which can be taken into account by introducing the 
~~ S . > cne . ‘ . » ws 
“2 tion lope Into the — - ——— cle oad If a cross section at distance Y from the outlet has 
ys > P . eteas ors P Arie . . { 
a on for un ao =a r* Mamy Varied; 4. following values determined: A. A. fC V.1 
ari <) UEC oo ol ee cia it of the steep nega- and if the resistance effect is introduced. the height 
tive wave, because of the downstream increase of h will decrease along the reservoir, with the slope 
level due to the resistance head as the wave pro 
eTeSses upstream Hence. the effective potential 
. h vn2t(V+i 
energy of wave energy to be transformed into kineth . Z 
energy of accelerated water will be smaller than the j hi 
total energy corresponding to the initial height 4 
The first aspect of resistance can be taken into where the \lanning formula for the flow resistance 
account when the routing procedure is based on the used in metric svstem, where vn is the roughness 
two de Saint Venant equations It is, howeve coefli@ient and FR the hvdraulic radius” of are 
possible to take into consideration the resistances A—A At the distances AX/L) the new 
effect on steep negative wave movement when either wave height 4, is now 
the continuity equation ol the type P—Q=dS/dt 
or both the continuity and dynamic equations are | / j 
used for the routing The effect ol flow resistance / h [ ih “4 m 
on the fictitious inflow hvdrograph, caused by thre . : 
steep negative wave movement, will be dealt with . 
? 
in this Investigation where fois the height near the outlet, Ay the hemhnt oO - 
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m by IS > m ] first approximation / 
neglected: 2 fictitious inflow } jdrogra pI 
seep negative wave 78 considered as the in 
fhe movement of steep negative wave is talk 


and A, the total resistance head 
from the outlet to the section a Since (’ and R are 
he functions of A,, the solution has to be made by 
approximations. As fh, changes slowly, 
soften sufficient to start with 4 constant for the 
omputation of f° and R, and then values 
0 obtained by successive approximations 
In fact, the potential energy of the steep negative 
rave Is not all transformed into kinetic energy of 
ecelerated water, one part of it being left as potential 
ergy to be used as the flow resistance head 
The integration of the last term of eq 24 
from the outlet to the place where 
0 would give the total resistance head. The 
olume S, which is not accelerated by the wave be- 
part of the total volume 


vave at section wv, 


SuUCCeSSIVEG 


eorrect 


along 


reservoll 


ause of the resistance os 


) fig 10) corresponding to wavefront height A. and 
S B h as, 25 
ere Zo Is the place where ¢ l°—0 and BP. 1s the 


ean Width of the area element for total resistance 


ead h at the section g 

lhe analytical Integration of eqs (24) and (25) is 
t tractable even in simple cases like the pyramidal 
of rectangular figure 9. be 


values and R for a section a depend on 


SeTVOolr eross section, 


tUSst the 


esum ol previous dh ; ana because of the expression 





500 600 700 800 ; 300 1000 

, Sec 
River Reservoir for the rectangular opening 250 
jdrograph, when the movement of all waves 1s 
when the water accelerated by the movement oj 
flow second app ori~mation h jdrograph, when 


en into account 


The practical procedure would be to compute U 
and R based on constant value / along the reservoir, 
and to compute the slope of eq (23). By integrating 


"2 dh 
dx 


resistance head 
h along the reservoir is obtained Then h h h, 
gives the effective head of wavefront With 
values of 4, new values U’ and R are again computed, 
and eqs (23) and (24) can be obtained as a second 
approximation. Succeeding approximations can be 
made by continuing the same procedure. 

Due to the fact that the total velocity (V / 
increases in the upstream direction while the acceler- 
ated discharge decreases, it is usually sufficient to 
compute fh, to the second or third approximation. 

With the new values A h h f(r) it Is how 
possible to determine a more exact fictitious hydro- 
graph P 

For the computed function h,=f(x) the width B, 
determined since B,=f(x). The integra- 
tion of eq The product of B.A, along 
the reservoir makes possible the approximate compu- 
tation of the mean distance f, of total volume S, 


from the level I] h 


h ly the of total 


function 


these 


has to be 


25) gives S, 


} = - B h? dr. 26) 


where 
volume S 
weight of water 


is the potential energy corresponding to the 


and the resistance head and y is the unit 
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The energy of the steep negative wave (fig. 10) is 


Lu =vthSr, 27 


where 5, is total volume of reservoir between depths 
H1,—h and F1,, and f, is the distance from the center of 


gravity of S, to the //,—A (fig. 10 
The ratio of energies 
kf 
( : 28 
Be. i 
can be taken as a measure of the flow resistance 
effect on the transformation of the front height of 


steep negative wave. 

As h, increases somewhat faster than P 
in the upstream direction, B_A, 
outlet upstream. The effect flow 
decreases the accelerated flow P at its 
values near the outlet and faster at its lowest values 
at the upstream end of the reservoir. The fictitious 
inflow hydrograph will be now a straight 
line. Asa practical approximation for rapid compu 
tations this hvdrograph can be taken as triangulas 
shaped for a valley reservoir. 


decreases 
from the 
resistance 


highest 


Mcreases 
of 


slow Cl 


closer to 


3.3. Simplest Procedure for the Computation of the 
Steep Negative Wave Effect With Resistance 
Influence 


A simple procedure for the computation of the 
approximate free outflow hydrograph in this case 
would be: 

1) Compute & and ( ; (2) compute A, A), f, ¢ 


V. l’ along the reservoir for constant / te determine 


h, along the reservoir; (4) recompute all values under 
2); (5) recompute /, and all values under (2) until a 
sufficient approximation of the true values is ob 
tained; (6) compute B, as f(r); (7) compute S, for 
given h, and B, curves; (8) determine S, from the 


determine the difference S Pry’ 
triangular shape of fictitious 
with the maximum discharge 


DI 


storage function; (9 

10) determine the 
inflow hvdrograph P 
(), and the total time of fictitious hydrogra 


The sum of the original and of the fictitious inflow 
hydrograph F1IVeS the total inflow hydrograp! The 
new inflow hydrograph, the maximum discharge Q 
ut the new starting level I: fie ), the outflow rating 
curve and the function sufficient 
elements for any suitable routing procedure, based on 
the continuity equation only, of obtaining the free 
outflow hydrograph 

If the routing procedure involves the use of bot! 
the continuity and dynamic equations for unsteady 
flow which is gradually varied, the computed fictitious 
inflow discharges P,,, with their positions in time and 
along the reservoir, and the corresponding reservoit 
depths I, h, for the same times and places are the 
boundary be introduced in those 
equations. 


storage are how 


conditions to 


4. Discussion 


The computed outflow hydrographs from Valley 
reservoirs in the case of sudden wate release Can by 
only approximations of the true hydrographs, T! 
shape and the resistance to flow of valley reservoir 
as well as the water movement itself along them gy, 
usually very complex approximate, 
hvdrograph is to the true hydrograph, the more wo 
involved in the computation The digital COM. 
puters will, in some cases, the 
The use of computers, howevet 


The closer an 


Is 
TIVE Most accurat 
results 
only in cases where very reliable and sufficient record 
and data about the reservoir, outlet, and inflow 
available and where the great precision of results a) 
speed. of computation make their 
Therefore the use ol simpl humerical oy 
graphical procedures for the computation of approy 


IS JUST 


use economicalh 
feasible 
mate outflow hydrographs is often very conveni 
sensible, ana economically desirable 

The movement e steep negative 
the reservon in the 
Many acceleration 
approximations 
depend on the factors eit 


T? 


ol Wave along 


watel reservo 
that 


hvadrograp 


accelerates 


factors affect Lfiis sO 


ol the True outflow 


her taken into account « 


The accel rate cl discharges are conside re 


take 


nealee ted 
here as a fictitious inflow hydrograph to 


1 
De 


riven either when only the storage continuity e 
tion P-—Q=dS dT is used for routing procedure 
when both the conti lity and dynamic equations 
are applied The degree of unsteadiness of wa 
movement is the criterion for the use of only the eo 
tinuity equation or of both it and the moment 
equation for the routing The greatet the height « 
the steep negative wave ll relation to the reservo 
depths, the higher and more unsteady the inflo 
and thre oTreatel the flow resistance then the | 0 
will be the decree of unsteadiness of flow motio 
and the use of the two de Saint Venant equati 
with the computation involved will be more neces 


sary The developme nt of expressions [ol the qui 


titative determination o 
trentedtl tre 


t this eriteriol sa 


\ criterion for the ¢ 
movement of steep neat 


not 


problem 


of flow resistance on thre 


Wave Is given by eq (28 When the ratio ¢ is hig! 
than a selected limit. the flow resistance effeet 
to be taken into account regardless of the routing 
procedure 

The facet that the Lteep wave Is a discontinuity 
not suitable to be considered as an ul steady flo 
vradually varied, and the fact that 


moveme 


and transformation of steep negative wave Is 
affected by smal! downstream disturbances in 1 
reservoir, make the treatment of sudden wat 


release different from an ordinary routing procedur 
The determination of the outflow hydrograph cal 
be divided The effect 


stecp negative wave ae ance 


into two procedures ad 


i of flow resist on 
water acceleration, combined with the inflow hydro- 
movement: 


the 


graph during the period ol wave 


based on positon 


the routing procedures 
place and time ol accelerated dischare Ss and inflow 


discharges along the reservoir, as boundary conditions 
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4s the maximum outflow discharge Qo in the case 
of sudden water release of relatively larger openings 
is significantly lower than the discharge Visas which 
corresponds to the given reservor level and opening 
with already accelerated water In the reservoir, the 
maximum tailwater level for GY is much lower than 
for Qmax: SO that the free outflow hydrograph (out- 


fows without tailwater effect) is much more common 


than one would expect 
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Uniform Asymptotic Expansions for Weber 
Parabolic Cylinder Functions of Large Orders 


F. W. J. Olver* 


April 6, 1959 
of the differential equatior 


dew 
dt: 


u(t L)w, 


ure ought for larze values of |v, which are uniformly valid with respect to arg uw and un 
ricted values of the complex variable f Two types of expansion are found Those of 
first tvpe are in terms of elementary funetions and are valid outside the neighborhoods 

of the points ¢ 1, the turning points of the differential equation. The second are in 
erms of Airy functions and hold in unbounded regions containing one of the turning points 
The special forms of the expansions when the variables are real are considered in detail 


ind isvimptotic expansions for the zeros of solutions of the differs ntial equation are found 
reversiol Numerical ex imples are included 
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Part 1. Introduction 
1. Introduction and General Summary 


The purpose of this paper is to determine asymptotic expansions of solutions of the differ- 


ential equation 


lol laree values ol the complex parametet! uu whi hh are uniformly valid with respect to are’ pu 
and the complex variable ¢ 


The immediate reason for this investigation is that the solutions of eq (1.1) are basi 


functions in the theory of the asymptotic solution of the differential equation 


d we ‘ 9) 
dt i 


lor large values ol uv iha region in which p ‘) has two simple ZeTOs, just is Airy functions are 


< vsic Lal ‘ reddingt Middlese | 
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basic functions in the theory of e 1.2) in a region containing one simple zero of p(t 
; | / i 


knowledge of the uniform asvmptotic behavior of the solutions of (1.1) is required in) developing 
the asymptotic theory of (1.2 

Equation (1.1) is of interest, however, in other connections, for example in the study of 
wave motion inside or outside parabolic cylinders, and wave propagation in an inhomogeneous 
atmosphere. Hermite polynomials of order » are expressible in terms of the solutions of the 
equation when rv 2n - Accordingly, for the sake of completeness Thar \ results are derived 
in this paper which vO bevond those needed for the immediate purpose stated in the preceding 
paragraph. 

The desired asymptotic expansions are obtained by application of the general theory of 
the asymptotic solution of linear differential equations of the second order as developed by the 
present writer in [9] From the standpoint of the general theory, eq (1.1) 1s characterized 
by having turning points at / |. Asymptotic expansions in terms of elementary functions 
can be found in any region bounded away from these points and expansions of this kind are 
obtained in part 2 by application of theorem A of [9] 

In order to determine the behavior near the turning points, expansions in terms of Airy 
functions need to be used. In part 3 these expansions are obtained by application of theorem 
B of [9]; they are valid in unbounded ¢-regions which include one (but not both) of the turning 
points. All combinations of uw and ¢ are covered 

When the variables are real, eq (1.1) has important special solutions. Expansions of 
these solutions are derived in part 4 from the results of parts 2 and 3. Also included in part 4 
are numerical examples illustrating the powerful nature of the expansions 

In part 5 uniform asymptotic expansions for the zeros and associated values are obtained 
by reversion and illustrated by numerical examples Much of the analysis of this part and 
some in earlier parts is similar to that used in [8] in the theory of Bessel functions of large 
order, and advantage is taken of this analogy whenever possible 

Investigations of the asymptotic solution of eq (1.1) for large uw) have been made by 
several writers Their work on eXpansions in terms of elementary functions is deseribed in 
section 6 of part 2, and on expansions in terms of Airy functions in section 10 of part 3. At 
this stage it suffices to remark that the expansions obtained in this paper include practically 


all of the earlier results and extend them in various ways. 


2. Relevant Properties of Parabolic Cylinder Functions and the Gamma Function 


A comprehensive account of the properties of parabolic cylinder functions is given by 
Miller in the Introduction to [4]. For the purpose of reference we collect here the properties 
which will be required in this paper. We use Miller’s notation rather than that of Whittaker 
[13] and adopt Miller’s choice of solutions when the variables are real 

The standard form of differential equation for the parabolic cylinder functions is 


The principal solution U’(a,z) is determined by the condition 
U (a,z)~27* 4e74 as = ie y 
It is an integral (entire) function of 2 and an integral function of a. In terms of Whittaker’s 


notation 1),(z) for parabolic cylinder functions and VW >) for the confluent hypergeometric 


function, we have 


Other solutions of (2.1) are U(a,—z), U(—aa and [ a,—iz). The connection 
formulas are 


For fixed a and large 


l(a, e 4° {140(\2|7° ( jarg 2 < 47 e), (2.6 


where? the 0 is uniform with respect to arg Here and elsewhere « denotes an arbitrary 


positive number which is independent of all other variable Ss. Krom (2.4) and (2.6), we obtain 


Via e~ is ( 1+O(\2|-?)}+—— e114 O(\2|-? 2.7 
r(}+a 
valid when ba t-e<arg 2<ia—e. The same result with the sign of 7 changed holds when 
let+e<arg 2<fmr—e; the apparent discrepancy between the formulas in their common 
region of validity is merely an example of the Stokes phenomenon 
Let the >-plane be divided Into four sectors M : M . M and M . defined by arg zi * ior, 
le<arg 7< 47, arg >) <4aq and im Sarg 2: tm, respectively. Then eq (2.6) shows 
that for large 2, (7(a,z) is exponentially small in My, and exponentially large in M, and Msg. 
Equation (2.7 shows that (a, is also exponentially large in Mb, unless @ has one of the 
values —} s, —, , In which case it is exponentially small. In fact in this event 
n—} n—} (n=0,1,2 2.8 
The solutions U’(a, , '(—a,+iz) clearly form a numerically satisfactory set in the 
complex plane, because in each sector M, one of the solutions is exponentially small and at 
least two others exponentially large. Accordingly, we may expect that a knowledge of the 
asVMpLlolic behavior of U'(a, for large a) in the seetor arg zis ‘Tr, or even the half-plane 
arg <i, will enable us to determine the behavior of any solution of (2.1) over the whole 


-plane by use of connection formulas. 
Solutions when the variables are real. There are two forms of eq (2.1) of importance when 


the variables are real, namely 


dear te) 29) 
tie 

Pw a i‘. 
in which a and z are real. Equation (2.10) is obtained from (2.1) by setting z=ze-ji and 


replacing a by ita 
Real solutions of (2.9) are U(a,z) and U'(a,—r) but they do not form a satisfactory pair 


for all values of a; in particular they are linearly dependent when a ; ; 3 
The fundamental solutions are taken to be [’ aw and V(a.a . where 
¥ 
Viava l r)I(5+< ){sin ra-l (ava l(a ! 2.11 
\J 4 pt x I W 1! I ( 
I 4 
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For all values of a, U(awr) and Via are linearly independent When @ is negative it is often 


eonvenient to use in place ol V aa its multiple lila ] riven by 


When a-+-} is a positive integer ("(a,x) ts infinite 


Wronskian relations for these solutions are 


// [ adi | a LAL // / a / aoa 2 | ( : a ) 4 
The Hermite polynomials 
if 
Ha CO a5! 2.14 
ure related to the funetion / a by 
| 
H, (2) =24"e8"U ( —n—2,2y2 ) 215 
For eq (2.10), the fundamental pair of real solutions is taken to be W(a.r) and Wa, 
defined by 
Wawa 24 | ( Re ( [ awe ? 16 
Wia j 2/ke) 2¢ Im ¢ U (ava YAY 


where 


k= (1+ ) L/k=(1+e ) garg I (5+ia) 2.18 


The value of are T , ia) here is not the principal one but is determined by the conditions that 
it is a continuous function of a and equals a In a—a+ O(a as a 
The Wronskian relation ts 
MW Wia.r) Wea, I | 2 19 
Gamma funetion EL PAnstons We shall need the asVinptotic expanstol of the funetion 
{| ' ; 
r(54 ) and its reciprocal lor large Krom Stirling’s series » chaptet oS we have 
| l , 
il r( ) ln 7 ! Ps 2 .2U 
‘) ‘> 
where 
) ' -_ ' 
7 7s B i | i é | 5 | | = ” 29 
om, 28(28 | 4 IRN) 10320 
uniformly with respect to arg RB denoting the Bernoulli polvnom al Ol degres 2s Hence 
' 
| ( ) ) = -7 ‘ > ire r € 2b 


where the constants are determined by the asVmptotic tdentits 


M 


~exp Z 2 23 


or equly alent ly : 





Numert al calculation vields 


| | L003 - 1027 


. F : : a 
y —y , command Ys ‘ 995 

24 1152 $ 14720 398 13120 

Since Z is an odd function of 2, it follows from (2.23) that 
{= : } {= } ' " oe 
identically and hence that 
| r ~~ ) = 
~ ary T é ier 
r(3 Im)? s=0 


Part 2. Expansions in Terms of Elementary Functions 
3. Asymptotic Solutions of the Differential Equation 


In order to bring the differential eq (1.1) into the standard form required for the applica- 


tion ol theore 1h) A ol i9%l we take new variables é f ana a} hot the | adi of sec, 2 _ defined 


hy 
di 
( ' ) i ——_ ¥ if (?—1)* w 3.1 
dt 
ef. |7] see, 2 Equation 1.1) then becomes 
Pw , 
u j W, 5.2 
dk 
where 
d° >t 2 
Ais / 4 ca 
dk } / | 
Kor convenience we preseribe that “ 0 when f | Krom the first of eqs 3.1) we then 
obtain 
é 1?—1)? dt=3 t(t?—1 > Ini t+(t#?—1 3.4 
and on eNpansion 
| ..5 28 | 
1 1 1 1 \ - _ 
< af * In2? i ed aia Toy, > | ro) 


Temporarily we suppose that the branches of the many-valued functions occurring In (3.4) and 
5) are positive when t->1 and determined by continuity elsewhere. The relation 
be fe < f Lim 3.0 
follows immediately from (3.5) when f¢ >1, and thence by analytic continuation throughout 
the ¢-plane cut along the join of ¢- l and ¢ 
The mapping of the ¢-plane, cut along the real axis from —@ to +1, is illustrated in figure 1. 


The center diagram corresponds Lo the upper half of the t-plane and the right-band diagram 


to the lower half. They are of course conjugate to each other. Features of the mapping are 
that the upper side of the real f-axis Is mapped on three straight lines BA (Im &=0, Re é >0), 
AE (Re &=0, 0>Imé>—3r), ED (Ime br, Re &>0), and that the positive imaginary 
t-axis is mapped on the straight line OC (Im ¢ lor, Re &<0). 


Next, we make the substitutions (ef. [9], see. 1 




















0 ry 
D E’ 
c! 
(-plane E-plar 
FIGURE | Vapping of the netior 
(the wi here is not to be confused with the real variable of sec. 2? TI Cli o.2 ICCOTLES 
PW ‘so ; 
u é ’  , ( "2 W. Pte) 
da 
which has the form of eq 2.1) of [9 The relationship between a and ¢ is vivell by 
/ rif ( "t / ¢ / | df: ».4 


the z-map of the ¢-plane is obtained from the &map by rotation through an angle 26. The 


) , 


function .7 is expressed in terms of ¢ by odd For large ir) we deduce from (3.3) and (3.5) that 


} l6u 


uniformly with respect to bounded arg’ J and bounded a 


In consequence of the last result we could, if we wished, take the domain D(@) of [9], 
section 2 to be the map of the t-plane cut along the real axis from to | with the cireles 
t+] 6° removed, 6’ being a fixed small positive number The preliminary condition (2.2 
of 19] would then be satisfied with oa | Crreater regions of validity in the desired CX PANslons 
will result however, if we arrange that the boundaries ol D 7) are parall | to the positive nag- 
inary g-axis wher possible \ccordingly, we shall now mvestigate what we shall call the 
principal curves in the t-plane, namely the ¢-maps of the straight lines in the &plane passing 
through A (€=0) and E (& Sim) and inclined at an angle —24@ to the imaginary €-axis 
In other words, the principal curves are the level curves of the function « which pass through 
f | 


Near ft=1 we find from (3.9) that 


Hence three principal curves emanate trom f¢ :. and have as tangents the ravs 


Equation (3.6) shows that three principal curves also emanate from / 


For large ¢) we have from (3.5 


Hence the six principal curves are asymptotic to the four rays arg f in —8, tir—?. 
Figure 2 illustrates the principal curves when 0<é@<7. Those through A are the continuous 

curves AP, AQ, AS when 0<@< }a, and AP, AQ, AR when }r<6<r. When @=0, AQ 

degenerates Into the curve KQ plus the join AE, and when @ hr, AS, AR degenerate Into 

AE-+EsS, AE—ER, respectively The principal curves through E are the images in the origin 

of those through A; this can be verified by using the continuation formula (3.6). 

We can now define the branch of the many-valued function £(f) which will be the most 


convenient tor our purposes. When 0<@< 7 we cut the i-plane along the principal curves 
AQ and ER and select the branch of &(f) which is continuous in the cut plane and for which 


i. l 
¢ (0) itr A—()). alee 





/ | VARTET) ait C 
E 0 A B R, ZZ ON A 











When T 8@<0 the cuts are taken to be the conjugates of those for 4. and the branch is 
fixed by 


£(0 hs T i 0), £0 ()) T i () 2 
1 { 14 


‘ “wnder- 


It is evident from (3.6) that (3.13 and (3.14 are consistent when @ () Henceforth aw 
stand t and E(t to de note this branch. Kor a viven value of o T iu T £(/ 


E(?) is a single- 
valued function of ¢ defined everywhere in the f-plane except when @—0 and ¢ lies in the shaded 
region to the left of ER-+ EQ in figure 2(a). When ¢>1, &(f) is positive provided that 4 Sr: 
this result is false if 1 6 r for the cut AQ passes across the real axis as @ passes through hor 

Let a r(t) now denote the branch ¢?'’&(f), We then take the s-domain D(@) of Q). 


section 2 to be the hap of the fun tion «(7 from which circles Surrounaing the points \ and K 


and of radius $6 have been removed, 6 being a fixed small positive number 


The map of 4 has two sheets when @ O and three sheets when 0 0 1 It is illustrated 
In figure 3 for @ 0 and O a ; 7 The passage from one sheet to another takes place “cross 
the positive real axis and, when 0<@< 47, the join AE. The heavy lines correspond to thi 


cuts in the t-plane and are boundaries of D 0 
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The map of «(f) ts clearly a single unbounded domain In order to ensure that D(@) is a 
single domain we have to impose the restriction @ <2r—e, omit sheet (il) when @)< «and make 
§—d(e) small enough. Otherwise, for example, the circle center E radius 36 would intersect 
the boundary AQ 

Next, we define the s-domain G(@) required by [9], section 2 to be that part of the map of 
r(t) Which remains after removing all points at a distance less than 6 from its boundaries 
Sheet (iil) is again to be omitted when @<e. The boundaries of G(@) are the dotted lines and 
circular ares of figure 3; their /-maps are the dotted curves of figure 2. Clearly when |@| <r—e, 
G(@) is a single domain if 6=46(e) is small enough. We readily show that conditions (i) and 


(ii) of [9], section 2 are fulfilled (ef. [9], see. 30 


Following |9], section 2 further, we take a.(@) to be the pomt at Infinit, on the positive 


real axis. We then find that H.(@é)=G(é From theorem A it now follows that there exists a 


solution HW, of eq (3.8) such that for large positive values of wu 


W,~. >> : re Gio 0! <r—e). 3.15) 


— 


uniformly, with respect to v and 6 The coeflicients are given hy the recurrence relation 


[ie -. Ry 


2 da 2 


the path of integration Iving in D(@ 


If we now restore the variables w and & by means of eqs (3.7), and let 


the relations 3.15) and >. 16 become 


, A, (E ; 
Wine »: te S(argy), arg u| <w—e 3.18 
as Lu Where 
| Ty 
/ é 5 A (EV— 5 |-F (E)- Ae (EAE 2.19 
Here S(are p S(#) denotes the ‘-map of G(@): its boundaries are the dotted curves of figure 2 
It is of some interest to determine the asvmptotic behavior of the boundaries of S(é@ 
Consider, for example, the dotted curve AQ lts equation may be written 
be yi ‘Oo¢ 3.20 
wher thie real parameter y ranges irom 0 to Reversion of (3.5) wields 
| | 
/ 2&4 ln Sé+ POs nh: as |§ 5 PP 
eT= eT é 
When y is large, we mav substitute (3.20) in (3.21) and thence obtain 
> cos 24 In s) cos 264 2A r) sin 26 16 In ) 
are 7 7 H+ ae - +_() (  ) , 22 
} s / SY y 


The corresponding equation of the principal curve AQ is of course obtaimed on setting 6==0 


in this result. Clearly the principal curve AQ and the dotted curve AQ, ultimately always lie 


; 


on the same side of their common asymptote, except when @= $2, in which event AQ its its own 


asViInptote 
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4. Identification of Solutions 


The substitutions 


l 
ad > b ive | 
transform eq ) ie Into | This result and the second ol eqs how that 
w=(e—1)tu (—* y2pty2 
5 MMI ) $2 
is a solution of (3.2 Equation (2.2) shows that when yu is fixed and ¢ ( ** this solution 


becomes exponentially small The same is true of the solution H ntroduced in Secllon 5, 
provided that w is sufficiently large and arg uw) <a—e; this can be seen from (3.18) and (3.5 
Since all solutions which are linearly independent ot (4.2) are exponent ally large in the same 
e solution (4.2) is a multiple of 


circumstances, we see immediately that tl 


Hence we have 


as jul->o, uniformly with respect to f and argu. The branch of (7 here and elsewhere is 


understood to be continuous in S arg’ u and to be asvmptotie to the pl neipal value Ol ,/ as 


1 

t{—> wo ¢~ ‘are Similarly for the branch of (f 1)2 used later. The function u) is independent of ¢ 
Its asymptotic expansion for large |u) can be found from the limiting form of the ratio of th 
two sides of the relation (4.3) asf ( é Thus with the aid of (2.2) and (3.5) we find that 

jiu P P —s 2 
~2Z ev” y > arg py "i t.4 

as |uj->©, uniformly with respect to arg uw Here 
lim .A(é as ¢ f 1.5 


That this limit exists and is independent of arg uw can be seen from the lemma of [9], section 5 
The J, are not determined however until the arbitrary constants in the recurrence relation 


3.19) have been fixed, and we now consider how to do this in the most vwdvantageous Manne! 


From (3.19 3.3) and the first of (3.1), we obtain 
1 j y, | ° 4 » 
” ») r fs" +.6 
-(f | SJ (f 


Taking .A 1. we find 


/ + constant 1.7 
24 (1 


A convenient value for this constant Is zero (nen y has the form of a po vnomial in? divided 
3 


by {- 1)2, 


The ec 1.7) suggests the substitution 
{ 


y, w(t / tS 
in (4.6 This leads to the recurrence relation 
f L)u.(t stu (t } f $.9 
where 
Nr (f if 2)u(t 12(s L)f f 1 (7 | $10 


We can satist\ this relation by taking u.(t) to bea poly nomial in ¢ of degree Os which is an odd 
or even Tunction of ft according as s 1S odd or even All the coeflicients are determined aHutlo- 


matically ifs is odd, but there is a degree of freedom if s is even, since the left of eq (4.9 


vanishes when w, is a multiple of (#?—1 2". We remove this by making the coefficient of the 
highest powel of t vanish 
Thus we prescribe that u,(f) is a polynomial in ¢ of degree 3s (s odd), 3s—2 (s even, s>2). 

With w(t | this condition determines the coefficients completely Numerical calculation 
based on (4.9) and (4.10) yu lds 

Molt l, u(t f 67) /24, u(t Q/*+-2497 145)/1152, 

(4.11 
Ma (t 10427" 1SIS9t‘— 282871t?— 1 51995t? —2 592907) /4 14720. 


Hence 
2021 

; Sonn 4.12) 
24 °* 2 07360 

An independent method of calculating the asymptotic expansion of g(u) for large u will be 
given in section 5 

Krpar S10 Tol the derivative Theorem A ol 4y shows that the asVmptotic series 4.55) 
may be differentiated term by term with respect to ¢. With the aid of the first of eqs (3.1) 


we deduce that 


i 1, ANS ; 
( > uty2 )~ > Ul ph f l)te * » > fe S(arg pw), arg ui wwe 4.155) 
- \e~- * la 
as , uniformly with respect to ¢ and arg uw. Here, in the notation of Miller, 
/ a Ol a O 


The coethe hts are Given by 


BE p(ty/(t | 4.14) 
where 
l 
f i(t)+5 t / ) / 1.15 
The first four f) are 
f | / / 67)/24 / L5t iT 145 1152 
1.16 
f 10497 IS1S9¢ 363871 2 384251 2 592901) /4 14720. 


5. Use of Connection Formulas 


In this section we seek expansions which hold in the complement of the domain S(arg ug 
Since © (a,2) is an integral function of a it will clearly suffice to obtain a set of expansions which, 
together with $3). cover the whole t-plane except of course the neighborhoods of the points 

1) when arg wi <2 

Kor convenience we introduce the notation U U.(6 j 1.7. «4 for the closed 
domains into which the t-plan is divided by the principal curves defined in section 3. Figure j 


shows the enumeration of these regions when 0<@<7 When r<60<0 we define 


Where the star denotes the conjugate region Then the boundaries of each region U,(@) vary 


continuously with @ in the interval 1.7 


14] 





- 
f 

\ 5 

U, e\, A q ; Lu, “TA E <A t U, 

3 ° J Us 
i=) ~ 
FIGURE } (-plane domain l a) 
The following relations, suggested bv figure 4 


ire easily verifies 


‘ 
' 
‘ 
‘ 
‘ 
\ olin 
Vo ‘ ae 
‘ -— 
= ' i? 
a "2 
md iter af 
aa ‘ 
o— ‘ V 
‘ ‘oO 
‘ 
C(VoUV>) . 
A 
‘ 
‘ 


The region formed by adding to U.(@ 


points whose «-maps 
less than 6 from the «-map of the boundary of I 
illustrates for example V,(@), V.(@ 


A 

and '. V 

and V.(@), when 0<@<14 The boundaries 
of figure 2 


cl. fig. 5) are at 
Wwe shall det ote 


by V.(@) or V 
the compleme 


a distance 

Kicure 5 
fi V ti 
are of 


t of the union of V 
the same charact 


# 


? 
( 


as the dotted 
The first connection formula we uss 


eurves 


U(- = a yt? ) é il U (=u, if 2 )4 
obtained from (2.5 Replacing u by 


Uiya« 
win (4.3), we have 
' u 
U(5u iuty2)~g 
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, l ‘ ° f ‘ 4 </ - | mf = = 
U(5" iuty2 )~g io > ay (t¢S( +57), =e + eS 65 T e), ».6) 
‘ Ml é 2 


2 t?—1)4 2 
Q und from (4.4 and (4.12 
| } 
ai > * ¢ IS 5* € 
| acti | ” “eo ( I—2>) ) (5.7) 
s { dA 3 :' | 
su € f 57 € 
, , ace _ ; 
When @ : 3 1 e, We can substitute (5.5) and (5.6) in (5 } (‘are is required, however. hecause 
at anv given point in the t-plane the branches of &—£(t) and (f 1)4 may not he the same in 


the two formulas 
When 0<@< 32r—€, S(@— 3x) is the complement of the region enclosed by ‘‘dotted curves” 


surrounding the principal curves AS and ER of figures 2(a) and 2(b 


S(é ir is the comple- 
ment of the region enclosed by dotted curves surrounding AP and KQ To the right of 


AP+AS, ie., in U,(@), the branches of &(#) and (#—1)* are the same in (5.5) and (5.6). Sub- 
stitution of these formulas 1 (5.4) leads to a hugatory result, however, as we would expect, 
' 
: since @* is exponentially large here whereas / pm, wly2) is exponentially small. 
In the central region between AP+ AS and KQ+ER, ie., in C(U,UU,), the two branches 
of c(t have opposite signs Let 
in in (5.5), Ein in (5.6). 5.8) 
J Then » is real and positive when 1<f< 1; from (3.4) we find 
, | l - 
n | dt = cos” *t = #( | / 2.4) 
The appropriate branches of (@—1)4 and .«<4(£) im (5.5) are given by 
| elt (1 7 A(t) —=« u.(t)/(1—t 10) 
ef. (4.8)), in which (1—? and (1—?)** have their principal values when (<7 3 For 
5.6) we change the sign of 
Carrying out the substitution, we obtain 
1 » -9 aa 1 u f ) 1 = Uy f 
( 2M”, hye) ~ 5 | COS (wn 7) 2y Sin (u*n Lie = 
|] — #2) 4 =o (1—t*)™ py =f (1 —f" uw 
5.11) 
as yu uniformly with respect to t¢ C(V,UV.) and @ in the interval »w +e, 5r—e). Here 
9 1 1 ) hy 
e du ~ (27 ] > Ml 24 ip4 | > 
5 
7 . 7 1 
~? . ( LM S i | t] 7 € 
Ss § 2K 
5 12) 
on substituting for ['(3+- 342) by means of (2.22 
The expansion (5.11) is one of the required formulas. Before leaving it we shall show that 
J(u d\p), and in the process obtain a new derivation of the coefficients in the asymptotic 


expansion (4.4) for 1/q(p 
When ¢e<@<39—e the principal curve AQ is included in the region of validity of (5.11) 


this is easily proved by considering the z-map. On this curve y’n is positive. Let us set 


n ” ( m ti) (5.13) 








where mis a positive integer Then 


where _ the value ol f correspond 


Now keep u fixed and let » 


From (3.5) and the first of (5.8 


Hy hee 


5.1] rt 


ne to 


Then 


“luces to 


and substituting these results in (5.15), we obtain 


as the limiting form of the left-hand side of (5 


On the right of (5.14) we have 


ef. (5.10). (4.5) and (4.12 (Compa 


where we use the symbol ~ to meat 


totic expansion for large 


We now substitute (2.27). wit! 


Although this result has been proved 


5.12) shows immediately that it holds in the r: 


~~ 
—“— l 


identically. Thence using > 2b) and 


Multiplying this identity with (5.19 


Q) () if 3 > ] we see that we can take 


ring the 


on thre 


and com 


) 
f 
=p ~~ : 
| } pee | 
( and so from (2.7 
ae 
- i f > 
(2 1 : \<- 
| 
. + ()( ff 
! 
) ‘ 4 
f _- 
4 ) 
=" } 
| 7 3 S 
(1 
I r 
14 
} | f ()} 
two sides we see that 
) 
1 1 » we 
two sides of this rel ive e sSnme 


in (5.18) and obtan 


) +) 


— 
— l 


we 7 €, il 
tive —yor] ol L Wie 
| - 


paring the result w 
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é Compal 


t.4 remembert 


».14 


» r 


wit! 


that 


1p- 


it! 


in the expansion (5.11), as stated earlier. 


We can express the coefficients g,,; in terms of y,, by equating coefficients in (5.21); thus 
. 2021 
d . ’ qd 2) . — —s 
24 vis 2 07360 
on using (2.25 These values agree with those given in (4.12). 
The expansion (5.11) is the first of our expansions for L’(—4 yu’, uty2) in regions comple- 
mentary to S(é The corresponding expansion for the derivative is 
, 1 ») ) i : l \ t 
| rt uly 2 ~My Hp | f . Sin (un sm P 
; 1—ft uu 


+ COS( un ul + 2. P “= = } edad 
® ai os Te = 


valid under the same conditions as (5.11); the coefficients v,(¢) are given oy (4.15) and (4.16). 

A further expansion can be found by substituting (5.5) and (5.6) in (5.4) when ¢ lies in the 
region (U,) to the left of EQ+ER. We shall not carry out this substitution however, because 
the eXpansion so obtained is contained in a more reneral result 5.28). w hich we shall now deri\ e. 
(sa preliminary we record the follow ing relations in which the symbol ~ has the same meaning 
asin (5.18 


() | pe "\wer\? 3 \ pu 7 e~ arg us €), (5.24 
1 
Qe) tein Caw+ 4 os 
i “aE (]\ pe —~ 4 COB muse \ arg pu 5 € 
2 2é ~ & e OF 
(eae 
— ¢ oat Le esargd uw T € 


We use the connection formula 


1 1 1 
| | ee ( "Se "s | . = . 
l (—5 2’, ut 2 ) . his (—; u*, uty2 )4 ‘1 i U(5«, iuty?2 )s (5.26 
2 2 r(a—tn 2 
obtaied from (2.4 We may substitute the expansions (4.3) and (5.5) when 
lt € <0<.7 € 


The branches ot €(f) are the same in both expansions when 





» ) 
se r y, 
Using (5.25) we obtain immediately 
| q , A(=) ) 
2 ) i T os T ‘ “s » > = 
U( till uty?) > ; | (sin 5 = _ 3 ™" ) = 
. INE 
2? cos Tu  & ] 
> 28 
| 
te C'(V \ ( -T+e A 0) 
te C(V,UV;UV,) (O0<0< 5 x—e) 
> 4 











and 





- 4 
The last result That \ be expressed in a form more pe rtinent to thre | ‘ 4 > 1 ly re placing 
uw by we’ and using (5.24 This gives 
7 
U( n2)~or nile 2 2 || 
» * v¢ _ — 
es } })4 
5 () 
, | 
e f \ \ ( , a e ) | 
a J 
Three further formulas, (5.28* Dae and (5.30*), can be d ved immediately from 
5.28 . D 29 and 5.30 respectively, mereh by chat (ring thre sigh ol The ire valid with 
the conditions 
| 7 
eC(V \ Uf ( > e+ es9 0) 


lor (5.2 


te C(VoUV (<<<. 7), for (5.30 


. Summary of Expansions in Terms of Elementary Functions 
6. S a fE T f E] t k t 


The principal expansio1 of thés kind ts (4.3 It is valid for large uw uniform] with respect 
to arg pu Aand ¢t when Tv € i T e and ¢ lies in the unbounded domain S(@), the bound- 
aries of which are the dotted curves tllustrated in fieure Z The hunetio £(7) is Given D\ 3.4 
that branch being chosen which is continuous in S(@) and takes the value 9.1) On » 14) at 
the origin The branch of (7 In (4.3) Is continuous in S(# and is asvmptotr to the prin- 
cipal value of yf as? ( Both and (f are positive when f>1l and @ T The 
coefficients .A(£) are given by tS) and (4.11 the branch of (7 aT defined similarly 


The function g(u) is cali ulable from either of the usVmMptots eXpanslo 


both of which are uniformly valid with respect to arg uw in the inte we. r—e Here 


Qi, gz are given by (4.12); ¥,, Yo, ¥2. ¥4 bv (2.25 


Expansions valid in regions complementary to S(@) are (5.11), (5.28), (5.29), (5.30 


(5.28" 5.29*), and (5.307 When 0<4: aT the eXpansions 5.28) and (5.29) between them 
cover the principal curve ER of figure 2; expansions (5.28*) and (5.30*) cover the principal 
curve AQ, except when 0<é@<e and Re #>—1—é’. The region of validity of (5.11) however 
includes that part of \Q to the mght of Re t> |+-6’ when 0<@<e. Hence the eXpansions 
we have obtained cover the whole ¢ plane except of course the neighborhoods of the points 
f 1) when 0<@<4$a, and similarly also when Lr<a<0 

The branches of the functions &. (f 1)* and ./.(£) in (5.28). (5.29). (5.30). and the con- 


jugate res its are the same asin (4.3), and the function g(a) again has the asvmMpLlolle CXpanslons 
(6.1) and (6.2 

In (5.11 uw) is identical with g(a), and n= (ft) is given by (5.9 The branch of the latter 
function Is contimuous in the region of validity of (5.11) and positive when l / 1: the same 
is true of the branches of (1—f)* and (1—f 

The expansions for the derivative 0” >, wly2) corresponding to (4.3) and (5.11) are 


1.13) and (5.23) respectively, the coefficients -4,(~) and 7,(t) being given by (4.14) to (4.16 


The derivative eXpanslons corresponding to (5.28) to (5.30*) have not been recorded, but they 


can be deduced immediate ly by use of the formal identities 


and are valid under the same conditions as the parent formulas 
UCniform asymptotic expansions of solutions of eq (1.1) for large uw) and unbounded ¢ 
have previously been given by Watson [12], Schwid [10], and Darwin [1] 


Watson applies the method of steepest lescents to integral representations ol 


}) U7 ( —n—5) z ) 


and obtains a set of asvmptotic expansions for the function D 2oyn when |n! is large and 


arg ) lx e. Which cover the whole (-plane other than the neighborhoods of the points 
c For each expansion the ¢-region of validity us fixed, unlike the regions obtained 
in the present paper which are perm'tted to depend on arg 7 For example, the region of 


validity of Watson’s expansion analogous to our principal expansion (4.3) Is 


This is, as We would expect the region Common to s ary uw) as arg’ pu varies over the interval 
I 


4 / t , / t 
Schwid appli s the methods olf Lane r for the usVmMptlotl solution of differential equations 
and his results resemble those we have obtained The main differences are that Schwid confines 
attention to the regions Re#>0O and are u iar, gives only the leading terms in each expansion 
and identifies the asymptotic solutions by considering their behavior at #0 rather than t= o 
Schwid’s results apply to the even and odd solutions of the differential equation (n 
effect of carrying out the identification at ¢—0 is that the error term in the asymptotic repre- 


sentation of each of these solutions Is obtained essentially in the form 


é Olu OO ( 6.4 
in th present notation if we try to deduce from these results the asymptotic form of the 
function / aT if ty2 by means of connection formulas, the error term will necessarily also 
be of the form (6.4 \ccordingly, a nugatory result will be obtained in the region (Up(arg u 
In which / bu, uly 2) is exponentially small for large Indeed, the proof given on page 
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358 of Schwid’s paper of the fact that the error term in the asymptotic expression for the Her- 
mite polynomials is of the form e~**O( u!~*) instead of (6.4) is false. The mistake is apparently 
due to overlooking the fact that the O’s are functions of the indep ndent variable as well as of 
the parameter yu. 

Darwin develops series solutions of the eq (2.10), and in [4] Miller applies the same methods 


to eq (2.8). The results are typified DY the series 


} L a | - 3 j 
l(a. = eXp { ae In .V +4 : > \ } 6.5 


for negative a and real positive So, in which 


: l . 
A V¥ (2 hia). d= Na 


and d;, is a polynomial in z of degree 3s if s is odd, or s if s is even. Explicit expressions for 
ds, ds, . . . , dog are given on page 68 of [4 

The series (6.5) is derived by purely formal methods from the differential equation with- 
out any investigation of its asymptotic nature. Applving the substitutions a yu and 

ut, 2 we find that Y= u(#2—1)2,2 and 6=,7£(t), in the notation of sections 3 and 4 Next, 
replacing {I'(}—a)}? r(4+4y2)2 by its asymptotic expansion for large positive yw (ef 
(2.22)) and expanding exp : d 4 iN) descendmeg powers Ol ¥ we find that the series 
(6.5) reduces to the form (4.3 Henee (6.5) must be a uniform asymptotic expansion for 


large a) of the same character as (4.3 

Darwin remarks that the explicit expressions for the coefficients ¢,, are simpler in form 
and easier to calculate than the expressions for the coefficients in series of the form (4.3). 
This is indeed true, but we cannot however regard the form (6.5) as being entirely superiot 
to (4.3); it does not, for example, lead to convenient forms of expansion when we differentiate 
with respect to 


Part 3. Expansions in Terms of Airy Functions 


7. Asymptotic Solutions of the Differential Equation 


In this section we seek to apply theorem B of [9] to the differential eq (1.1 The prelimi- 


nary step is to take new variables ¢=¢(f) and W, defined by 


de | di — aa 
(ae) é, , HN (“") v ( ) t 7.1 


(ef. [7], see. 2). Equation (1.1) then becomes 

PW 

[ We 7.2 

de 
where 

it d dt ) ) 2 == 
f ; | I« ) iC {(‘.) } It = ib } 

on reduction (ef, Ded 


The first of eqs (7.1) may be mte 


on specifving that C () w] en f | In terms of the va inble © ot sectiol 
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An essential difference between the functions &(¢) and ¢(¢) is that the former has a branch 


point at ¢=1, whereas the latter is analytic at this point. Both functions have branch points 


at? 


g=ar>re wu 
principal curve ER of figure 2 when 0<(6<r, along the curve conjugate to ER when—21<é< 0, 


The branch of ¢(f) which will be the most convenient for our purposes depends on 


We define it to be positive when ¢>1 and continuous in the ¢-plane cut along the 


or along EQ and ER when @=0. This determines ¢(f) uniquely everywhere except in the 


region (U.(0)) to the left of EQ+ER when 6=0 


Following [9], section 1, we make the further substitutions 


Equation (7.2 then becomes 


TT) e¢ mapping of the i-plane on the >-plane is conveniently carried out by using rf (defined 


in sec. 3) as an intermediate variable; in terms of 


In our first application of theorem B to eq (7.7) we make a cut in the 2-plane along the 








part of the level curve of the function exp ($27) which joins the point E of affix—({m)%e*" and 
corresponding to f 1. to the point at infinity on the negative real 7-2 X1S Such a join can 
be made when @ <4 (ef. [7], fig. 2 
FIGURE 6 >-plane: domains 
, l 
D(@) and G6 ( 6 <5" e) 
Figure 6 illustrates the mapping, the lettering corresponding to figures 2 and 3. The cut 
is the heavy curve KQ The principal curves through A are mapped on the ravs arg < us 


and Tr; t} c principal curves through Ke are mapped on the two sides of the cul EQ and the 


= « 


remaining part of the level curve of exp (322) passing through E. Cut in this way, the 2-plane 
Is mapped on the complement of the domain U. e] see fig $) in the t-plane 
When @ <12—e. we take the z-domain D(@) of [9], section 9, to be the cut plane with the 


circle center E and radius 36 removed, 6 being the fixed small positive number introduced in 
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section 3. Then Fi >) is regular in D(@) (Ginelnding 0), and the condition (9.2) of [9) 


is satisfied with o S because 
( I ( ~ - ( 
64 io 
as 2 ©, uniformly with respect to bounded arg 2 and bounded @; this result follows from (7.3 
7.5), and (3.5 


We take the ~lomain G(@ required by 9], section 9, to be thr omplement Ol thre map of 
the t-domain V ), defined in section 5 The boundaries ol G 7) are the dotted Curves Sur- 
rounding the cut EQ in figure 6. All the preliminary conditions of [9], seetion 9, are then 
satisfied if 6=-6(e) is taken to be small enoug! 


Taking the Olt aA of 4! section 9 to he at mfinitv on thre OSITIVE real uXIs, Wwe ting 
| | in) 


H 0 G6 Hence from theorem B, a solution W, of (7 7) exists, such that for lares positive 


values of u 


1(0 Ai’ (u B&O 
a ~ Al ul “ S ( e Gia T é ) 
— Mu ” a / Z 
uniformly with respect to and A, thre coefficients AO and 2 .(¢ bei yroiven by recurrence 
relations derived from eqs (9.12) and (9.13) of [9 When the variables u and ¢ are restored. 
this result becomes 
+ ‘ Ai’ (yu . rb 
Wi ~Ai(pss +> off >> 7.10 
as us ©, when 
. l oe 
te ('V.(0 and , 1 é ON57 é sin 
where Ap(¢ constant and 
7 
B 4 ») \ | 4 I ( | 1c 
Ko 7 () 2 7.72 
_ im i ; 
er Bots | FOBASd 


the paths of integration lving in the Map of CU.(@ 


A second application of theorem B can be made to the differential eq (7.7) by taking the 


cut in the 2-plane for the function F'(¢ to be the part of the level curve of exp (4 joming 
EK to S, the point at infinity on arg 'x. This join is possible if 0<é@<7: the corresponding 
t-lomain is then CU,;(@ 

With similar analysis, again taking a,(@) to be the point , We may show that a 


solution of eq (7.2) exists having the right-hand side of (7.10) as its asymptotic expansion for 


large uw. uniformly valid with respect to ¢ and @ when 
| 
te C'V f and e: on 7 é 7 M 


ntegration now lving 


The coefficients A,(¢) and B,(¢) are again given by (7.12), the paths of 
in the map of C'U,(@ 

By considering the limiting behavior as ¢ D ¢ , we see that when «<@<}2—e, the ratio 
of the new solution to W) is independent of ¢ and has an asvmptotic expansion fo! large yp 
identically equal to unity Hence the expansion 7.10) for a} holds vith the conditions (7.13 
Similar analysis, or an appeal to symmetry, shows that this expansion also holds with the further 


conditions 
te l'V 0 and T e \~@ é 7.14) 
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nd 


We can combine the three f-regions of validity into a single region T(@), given by 





CV,, CV,UCV., CV,, CV.UCV,, OV, ) 
according as @ lies in the respective intervals 
> 1.15 
l l | l 
TT € . T € 5 T € é é€, € 6.5 T € UU €.7 €). 
- ~ - - 4 
T(@) is illustrated in figure 7 for 0<@<a2—e; T #) is the conjugate of T(@).. The dotted curves 
are ol the same character as the dotted curves of figure 2; their r-Inaps parallel the J-Inaps 
of thre prin ipal eurves through Ie ata distance ouway T 0 is the unshaded domain bounded 


bv the dotted curves. 
Collecting together the results of this section, we have shown that there exists a solution 


W, of the differential eq (7.2) whose asymptotic expansion for large |u| is given by (7.10) and 


‘s uniformly valid with respect tot «T(@) and @ in the interval (—x+e«2—e). The functions 
c(t), FiO), ACO and &,(¢) are given by (7.4), (7.3), and (7.12), and are regular when t ¢ T(@ 
R 
Q. \ 
» »\ 
ys ° ff Q Pp »\ Q 
, ; * 
4 4 \ 
> \ 
+! E \ E é. \ 
E » a BA SS \ A 
0 -) je 0 A - 0 \ 
> 7. %2 ‘ \ > 
9 42 4 
- ‘a 
- \ 
y R \ 5 
\. \ 
Me f ‘ 
FIGURE 7. ¢t-plane: domains T(@ 


8. Identification of Solutions 


From } | and the second of eqs r ee , we see that 


. t-—1] | 
W=(—— ) U( —5 #uty2 ) 8.1) 
satisfies the differential eq (7.2 When wu ts fixed and ¢ oe ~** this solution becomes ex- 
ponentially small. If yw is sufficiently large and jarg «4 <a—e, the solution W)of section 7 is 


also exponentially small in the same circumstances, for from (7.4) we have ¢~(})*?t*, and 
1 


SO 3 Hence the two solutions are linearly dependent. Accordingly, 
c \ AE) , Ai’ (uD) SS BYE 
l ( 5M’ uty2 )~h P € ) Al (ui ¢ eu ’ T , hand i 8.2) 
Z ; iv bu s=) 
as yu ©, uniformly with respect to ¢ « T(arg w) and argu in the interval (—2+., r—e 


This is the fundamental expansion in terms of Airy functions. In order to determine the 
function A(yu), we first fix the arbitrary constants of integration in the recurrence relations 


7.12) by specifying the conditions 


Ay(¢ B lim A.(¢)—0 LL a - (8.3) 
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We now examine the limiting form of the ratio of the two sides of the relation (8.2) asf ei arg 
With the aid of (2.2), (3.5), (7.5), (8.3), and the asymptotic forms 


Al M C)~ i wT A(t 7 wih! Ai’ uso ~—3 T rc) te 


(cf. [8], appendix), we find that 1/A(u) has the asymptotic expansior 


P Ly bi — 
“TT y a *u e* (1 S ) ary uw! nT € S.4 
h u , vi 
as iy , uniformly valid with respect to arg wu, where 


h lim {¢? Bie 8.9 


An alternative form of expansion, not involving the quantities A,, may be found as follows. 


Suppose that «e<@< 3r—e The principal curve AQ is then included in the domain T(@ 
On this curve uw’ is negative (ef. fie. 6). Let us set 

¢ ¢ Lu a, S.6 
where dl, 18 the mth negative zero ol the funetion Ai’. Then S.2 becomes 


U( “9 Hy ae 2 )~h(u (" )' Li (a) >> - , 8.7 


] 


In whieh / denotes the value ot f corresponding to ¢ ( Now keep u fixed and let » 
jz—6 —_ ; 
so that to. Then from [8], appendix, we have 
é..==- ri4m—3 1 Olm Ai (a. r a ‘f11Q0im RR 
4 
and from (7.5) and (3.5 
} | 
t= —In 27. c2+—-41(0(|¢, £9 
3 2 
Substituting in (2.7 by means of these equations, we find that the limiting form of the left- 


hand side of 8.7) 1s minus the expression D317). with 7 defined ol course as in the present Sec- 


tion. On the other hand, the limiting form of the right of (8.7) is 


ua, 
h Uo ( ) v a 7™ T Cs h(y)t 
fai 
ef. (8.3) and (8.8)). Comparing the two sides we find that 3 >? uw A(u) has the same asymp- 
rote expansion for large ui as the right-hand side of the relation y IS Hence Wwe obtain 
— I > 1 
h(u)™2r Med \ pp)? wt gp), S.10 
where J Mm), Gi) are defined in seetions 4 and 5: glu has the asvmptotic expansions (6.1 and 
6.2) for large w. Although (8.10) has been proved on the assumption that e<arg u<42r—e, 


8.4) shows that 1/A lu Sand hence also / MM), possesses a single asvmptotic expansion ove! the 
range arg uw <wr—e. Therefore (8.10) must also hold over this wider range of arg u 


Substituting (8.10) in (8.2), we obtain the final form of the fundamental expansion, given by 


Lh 


as uo, uniformly with respect to f € Tiarg uw), arg uw in the interva T-+€, 7-6 
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The coefficients A,(¢) and &,(¢) appearmg in (8.11) are determined by the recurrence 


relations 7.12 and the conditions (8.3 Explicit expressions can be found by the method 
ysed in [8], section 6 for the determination of the coefficients in the uniform expansions of Bessel 
functions of large order. We take ¢ to be fixed and positive, and let rv pom €D. The Airy func- 
tions appearing In (8.11) can then be expanded in their asymptotic series for large positive 
arguments and the result compared with (4.3). Thus we derive the expressions 
| < > > h ( 4 Pg Cc B ( be dd, ¢ Y £ (8.12 
ef. [8 6.6)), where .<4,(£) 1s defined in section 4, a l and 
2m | 2m Brace (OR | 6m-—-1 . 
a, P ~ , 0, =—-— a (8.13) 
m!(144 om ] 
Analvtic contmuation shows that eqs (s 12) hold for all values of ¢. 
We note from (8.10), (8.4), and (6.1) that 
h lim c2h ¢ 4): ‘ O.i. ods S.14 
where g, is given by (4.12) and (5.21 
kerpar S701 Tol the derivative The expansion S11) may be differentiated term by term with 
respect to 7; this is a consequence of theorem B With the aid of the first of eqs (7.1), we find 
that 
De) ty 3h Fe ( : Dc) 
, I . 7 Med \ pb Al ur’ ¢ ‘ ( , Ras ‘ 4 ~ 
( - = um, wty2 )~ : » > " + Aj Me 4 > ; 8.15) 
“ P\ 4 iv M s=) pM 
as , uniformly with respect to ¢ « T(arg uw) and arg yu in the interval w+-€,7—e 


Here we have used the notation 
. 1 di } 
Vee 8.16 
PAS ( = ) (-.) 16 
The coefficients C,(¢) and 1),(¢) are expressed in terms of .A,(¢) and B,(¢) by the equations 


( ¢ X(C)A,(¢ A, (¢ CBC), De A,(¢ X(OB 4 B: i(¢), 8.17) 


¢ ( ¢ . >> h ¢ Db» ” <=), ) ¢ > & Ay ‘Bs : g . S.19 
ef. [8], (6.12)), where .4,,(&) is given by (4.14) to (4.16), and a,,, 6,, by (8.13 


9. Use of Connection Formulas 


In this section we seek « Xpanstons W hich hold in ¢ regions complementary to T(arg uM) Ww hen 
arg wins 7 el. see 5 


In (8.11) replace w by we'*. Then using (5.24) we obtain 


U ( —5 w?,—uty2 )~2etlt ep a(u) (2°) 


153 








when 


TT €  « €, te Ti0+7n Q 9 
In particular, we see from (7.15) ane 5.3 that | holds whel 
l : 
rs6 € te C'V,(0+2 CV.(0 9.3 
) 


Similarly, 


when 
, 
> < 7, € ( \ ‘4 Q 5 
The expansions (9.1 and (9.4) are of course the expansions 101 »M”, wfy2) in terms 
of Airy functions valid at the turning point ¢ | Combined wit S.11) thev cover the 
whole of the t-plane when «< @\<1ig 
When GO6\€ we Cannol expect bo be able to represent the inction { u-, ply2) in 
V.(0 ('T(@) by means of a single Airy function, because as ¢ in this region [ u-, pty 2 


becomes exponentially small or exponentially large according as yu? is or is not an odd positive 
integer (cf, (2.7 
To deal with this remaining Case we use the connection formula 


+( 5 ) ~ — a (- 4’, uty2 ); 06 


obtained in an obvious manner from (5.26). Replacing uw by we in (8.11) and using the first 


of eqs (5.25), we find 


T f - | 4 
l ( L iut 2 )~2r*s ‘OS = Wun: ( ) 
| Ai ‘ = A 4 \ \) ( ur 4 ~ B ( ] 
We substitute this result and its conjugate form on the richi ot (9.6 and on the left we sub- 
stitute (8.11 Using the eq [3] 
Bi Lh? 4 f Alle Lh’ 4 f Ai (¢ a2er 


and its differentiated form, we obtain 


, I ¢ | a ' 1’ 
( ot; uty2 ) ~ rr Mh q Lu ( j ) Siti 7 Tu Al uc : & , | : > . 


| ‘ 1.4 Bi i B ( 
t COS ~ Ty Bi(uss > t » & ved 
on Li Li 
This is valid when /@!< s7—e 1N the ¢-region common to the domains T(@), T(@—! 7), and 
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T a] T Using 7 Ld and 5 ~~), we see thai ihis means (9 7 Is valid when 


A « r—<é te C'V.(8). QR) 


This is the required result We note that when yz? is an odd integer the terms in Bi, Bi’ 


in (9.7) vanish, as we require from a consideration of the behavior of [ Su’, wly2) as t\—>a 
in the sector V,(@). We may also verify without difficuliv that (9.7) agrees with (9.1) and 
9 4) in their common f-regions of validity when lo t+¢e<@< e and «<é<)r—e, respectively. 


10. Summary of Expansions in Terms of Airy Functions 


The fundamental expansion of this kind is (8.11 lt is valid for large (v4) uniformly 
with respect to @=arg yu and ¢, when r+e<0<97—e and ¢ lies in the domain T(@), defined 
by (7.15) and illustrated in figure 7. The function g(x) is calculable from either of the asymp- 
totic expansions (6.1) and (6.2 The function ¢(f) is given by (7.4), that branch being chosen 
which is continuous in T(@) and positive when ¢>1; similarly the branch of {¢/(@—1)}4 is con- 
tinuous in T(@) and positive when t>1. The coefficients A,(¢) and B,(¢) are given by eqs 
8.12); they satisfy the recurrence relations (7.12) and the eqs (8.3). They are regular fune- 
tions of € in the ¢-Ihap of T(@) and real when ¢ lies on the real axis to the right of the point 
of affix ‘*3 

Other asymptotie expansions for / bm, wly2) are (9.1), (9.4), and (9.7); they are 
valid with the conditions (9.3), (9.5), and (9.8), respectively. The combination of these 
expansions and (8.11) covers the whole f-plane when 6@ <$7 

The expansion for the derivative (°’(—3y*, ufy2) corresponding to (8.11) is (8.15), and 
holds under the same conditions. The coefficients C,(¢) and D,(¢) are given by (8.17) and 
S.19 The derivative expansions corresponding to (9.1), (9.4), and (9.7) are not recorded, 
but they can be obtained immediately by analogy with (8.15). They are valid under the 


same conditions as the parent formulas 

Approximations for parabolic evlinder functions of large orders which hold in regions 
containing one of the turning points have been given by Watson ({12], see. 17), Sehwid ({10)}, 
sec. 6), Erdélyi, Kennedy, and MeGregor [2], and Kazarinoff [15]. 

Watson derives asymptotic expansions in terms of elementary functions for the functions 
D 2yn ((—n—}3, 2yn) as no, with the conditions jarg n br and z=o(\n|~*). 
In the present notation these conditions correspond to arg pu trandt l o(\u!~*). Watson’s 
expansions are analogous to Meissel’s expansions for Bessel functions of large order and may 
he obtained from the present uniform expansions by suitable re-expansion (ef. [8], see. 6). 

Schwid gives approximations in terms of Bessel functions of order one-third (Airy fune- 
tions) which are valid in the neighborhood t—1=—O( u\~*) as ju), using the present notation. 
This neighborhood is complementary to the region of validity of Sehwid’s approximations 
in terms of elementary functions, described in section 6. No investigation is made of the 
validity of the Bessel function approximations outside the neighborhood. 

The approach of Erdélyi, Kennedy, and MeGregor is similar to the one we have used. 
These writers obtain a set of asymptotic approximations in terms of Airy functions, to solu- 
tions of eq (1.1) when ww and ¢ are both complex (using the present notation). The results 
are contained in those given in the present paper; the principal ways in which they have been 
extended and improved upon are as follows 

First, the approximations apply to fixed values of arg yu, no investigation is made of their 
uniform validity with respect to arg xu. 

Secondly, the approximations are established with the condition |arg u| <4}; outside this 
interval connection formulas are to be used. This leads to complexity in the number of results; 
eleven are given compared with the four we have found to be sufficient. 

Finally, only the leading terms of the asymptotic expansions are found; no terms involving 
the first derivative of the Airy functions appear in the approximations. An effect of this is 
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No special treat- 


to exclude regions in the ¢-plane containing the zeros of the Airy functions 


ment of such regions is required in the present theory 
which has appeared since the present paper was written contams 


Kazarinoff’s paper 
and NMieGregor to 


a brief investigation of the extension of the results of Erdélvi, Kennedy. 
include uniformity with respect to arg uw and higher terms in the asymptotic expansions. The 


results are rather more complicated and less complete than those viven here 


Part 4. Results for Real Variables 


In this part we suppose yu, f, a, and s to denote real variables except where otherwise indi- 
cated; uw is always large and positive. The variable « is not to be confused with the complex 


function 2(?#) introduced in parts 2 and 3 
ll. The Equation w’’=(1 s*+-a) u 


The principal solutions of this equation are (’(a.r), ((a,—.7), Vilar), and U(a.r), defined 
In section 2. 

(1) a positive When a > © the differential equation has no turning points on the 
those in terms of elementary 


real w-axis. The eXPansions of greatest interest are therefore 


functions, and we shall not consider the forms of the Airy expansions 


In (4.3) replace ut by au, —7t, respectively. Then using (4.8), we obtair 
, | , - ees —" i 
U (5, ut y2 )~g (ue > > 11.1 
2 j27~] | U 


as uw In the complex plane the region of validity of this eCXPansion Is ¢ Ss » 7), and 
from figure 2(c) it is clear that this region includes the whole of the real f-axis 


We have now to interpret the branches of the various functions 
From (6.1) and (6.2) we derive 
g 11.2 


(| pe ¥ ( oe - - dip), 2 


where the function g uw) is real and has the expansions 


as u- x Comparison of (11.4) with (2.27) shows incidentally that 
a ) ‘ ty 
d()o~2 tea” yy we : 5 
where the symbol ~ is defined in section 5 
Next, when ¢ lies on the negative imaginary axis of figure 2(« son the straight line 
O’C”’ shown in the third diagram of figure 1] Hence 
/ itr t f). 11.5 
} 
where &(¢) is positive when ?>0. From (3.4) we find 
£ / f ] dt ; fit | . In / 1] ‘ 


the branches here baving their principal values 


156 


> 


Lastly. we have 


i 1) 4 ge" (F" 1)4 4 l ¢ "(f l 11.8) 
where the functions on the right have their principal values 
Substituting in (11.1) by means of (11.2), (11.6), (11.8), and the relation 
u(t iu it 11.9 
we obtain the required result in real form, given by 
- e—welt . U.lt 
U(5 wuty2 )~G(u tos r<t{<_ @), (11.10) 

Z f? 1)¢s=o (f 1)2% we 

uniformly with respect to / Explicit expressions for the coefficients a,(t), obtained from (4.11 
are 
Uy(t & 71, (t t?— 6t) /24, u(t Of* + 2497 — 145) /1152, 
L1.1] 
Ua(t 10427" +- 181 89f? + 28987F 1 51995 4-2 592908) /4 14720. 
The corresponding result for the derivative is 
~ s a - ~ Olt | 
UO" (5 u?uty2 )~—-S glu) (2+1) te FO SY 11.12) 
“ 4 s=0 (f-+-] ’ 
valid with the same conditions Here 
7 .(f j if 11.13 
Explicitly, 
f | f tf? +-6t)/24 Do(t Lot 3270? +- 1435/1152 
11.14) 
f $0427" + IS1LSO¢ 36387 2 384251 2 592901)/4 14720, 


Since Ula and l(a. J 
equation when @ is positive, it 1 
solutions V(a.r) and l(a 


comprise a satisfactory pair of solutions of the differential 
S unnecessary to give the corresponding expansions for the 


ll ad negative. EL Panstons in lerms ot ele mel fary functions. For {> | é. the eXDansion of 
U(— 3", wty2) is given in real form by (4.3)°. In this expansion g(u) is calculable from (6.1 
6.2), or the relation 
> . 
y (ye 2 ie Qi igi wd 4+ ut) (11.15 
and & .°A(i , bv use of (3.4), (4.8), and (4.11 


, all fractional powers now taking their principal 
values 


When 1+6<t<1—6, the expansion of [ bu’, wly2) is given by (5.11). Here 9(u) is 
identical with g(u) and » is given by 5.9 

For the remaining part of the real axis we have the expansion (5.28). In the present 
CIrPcuMmstances the function £(f) 1s positive and bounded away from Zero, Hence the contribu- 
tion of the series ¢ cos dry?-e7** >) f4(E)w>* is always asymptotically negligible compared 
with that of 2 cos hip” + e » > Af) uy Thus (5.28) reduces to 
: ‘ (ys » a AA§ pwn AE 
(- aK; uty2 ) - —— | 2 cos 5 ~ee ES TSIM 5 we : ¢ i> > ; | 11.16 

2 p2—] 2 ; . 0 M 








valid when t> 1 0 The fun tiohs dip), & and 4.\¢t are the Same as ee . dise usseqd above 


The corresponding ( Xpansions fol the second solution [ ad ] nre Ot interest party ularly 


whenz>0O. It isnot possible to obiain them trivially by substituting in the connection formula 
(2.12) the expansions we have just obtained for (’(a, 2) and U'(a, because although can- 
cellation of terms in the asymptotic series takes place, we cannot assert immediately that the 
remainder terms implied in these series also cancel To put he resul oh a sound basis we 


substitute the CXPansions 5.5) and (5.6 in the formula 


obtained from (2.5) and (2.12 


When ¢>1-+6 the branches of &(f) and «A(£) in both expansions are the principal ones 
| | | 


With the aid of the relation 


proved by use of (2.22), (5.7 5 90). and (6.2). we deduce thet 
a | , : , 
/ ( > I uty2)~2 i \ 19 
When | 5<t<] 6, we interpret the branches of the man a 
of (5.8) and (5.10 With the aid of (11.18) we obtain on reduetior 
= | Pa(y * aT | | 
U(—5 0, uty2)~ > cos (uta+-ie)% on (ste 1 te) 8 
Z , , — (| : + > : 
11.20 
valid when 1+-6<t<] 
=e mane for the derivatives correspondi Yr to the expansions oDbtamed in this sub- 
section can be written down immediately by use of (6.3), or by analogy with (5.23 
“3 : negative expansions im terms of A / funet is Kor / the expansions of 
Ty yu? uty 2 and | “a. ul, ? mre TiIVel by XN 1] ana 9.7 respectively I) seen 
expansions Cf is given by P | and is real when l-n more conve ent form ol his relation 
when | f 1 is 
) . 
~ | ; at COR ; } | r 4 1] 2] 
: > 
(el. 5g : The coefficier ts cl ana BR re ren sna CLVel yy 819 Q | and 1s 
An expansion for the solution (’(a.r) mav be derived by use of 17). as in the ease of 
the exponential-tvpe expansions 7 ‘ result Ve 
> ; 
/ ( I ; »)~ = “ ( ) Bi \ | Bi < i} , 5 
vy \< é i r 7 ‘ S : 
valid when ¢>—1--6 
The corresponding CXDPDUANsLoO! IO! / 4 i ul, 2 we \ bye writtel dow! hy\ analoe with 
(S.15 
' Hermit polynon vals of large degree Expansions for the Hermite polynomial /7, (2 
MI | 
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when 7 is large can be obtained from the results in subsections (ii) and (iti) above by use of 


the equal 


Tae 


H pt) = 20 te eT | Su? yuty2 ) 11.23 


hus (4.3) and (5.11) vield expansions which hold uniformly in the ranges z>(1+6) (2n+1 
nT 2n-+-1 r<(l 2n+-1)?, respectively. An expansion uniformly valid through- 
out the t Ca 1—6)(2) is furnished by the Airy expansion (8.11 
12. The Equation w’’=(a—j.s*)w 
( wwilire CF PAnsions n terms of elementary functions Krom 2.16), (2.17 and (2.18 
we have 


Fo hare 


bd replace ubv uw Then we obtain 
| ( : u(t | 
U' (5 tae uty? )~g (ue ye 12.4 
- f | | Ma 
as i The t-r FIO) of validity of this expansion, S rs includes the segments f | 0 
ena | | > ot the real axis (cl fie 2 I) 
When g we have from (12.2) and (2.22 
| 5" 7 
ewe (5) [Ss s : 12.5 
y # hat (3 howe §, 7 
- ‘ MA 2's 
the chores ot smbicuous sigh on taking’ Live square rool beme resolved by use of 12.3 : Krom 
h.] 
(ue ~2e" ei es ielaalT | LL ‘ 12.6 


where the svmbol / denotes Lhe asvmptoti series 


L=>> 2.7 
Further from (5.21), with « replaced by wet”, we have 
- ) | yi 
We ie ene -° 12.8 


Introduce Ing a Tunction [(u). defined by 


yh dau } 12.9 


and usine (12.5 12.6) and (12.8), we find that 





Thus /(z) is a real function, having an asymptotic expansion of the form 


24 i] 
(g)~=2,—; a8n 12.11 
uss { es 
From (12.10), (12.7), and (4.12), we find 
/ l, / - Uk " ’ / qd ,qd + qe . |? » 
2° 1152 oe 8 ° 98 13120 sai 
We are now in a position to combine (12.1) with (12.4 When ¢>1-+-6 each of the fune- 
tions &(f), (®—1)* and (fF—1 is positive. Substituting by means of (12.9) and the first of 
(12.3), and equating real and imaginary parts, we find 
, . 2 ¢ idl ul : . u f ‘ uf f 
f- | ‘ - l bl / | u* 
12.13 
, Qres™ |] i Us (Tt : u f 
WW bu? uty 2 ~ ~ Sin ue ial pe : = + COS (ut 7 > > - ’ 
* f? ] 4 g = 4 l ue : f= l ue 
12.14 
uniformly with respect to ¢ when f 1+-6. The function /(u) is calculable from (12.11 
When —1+6< ¢ <1—6 the appropriate branches of — and (@—1)* are given by ¢ im and 
o 1 i, a > - . 
(f—1)*=e'"(1—f)*, where 7 is the real function (5.9) and (1—f)* is positive. Hence we derive 
° l . / Lu enn" . u | i 
W( Sut .uty2 )~ = >* 7 t<1—6 12.15 
“ Z7es™™ (1—f*)* « l—f we 
The expansions for the derivatives corresponding to (12.13), (12.14), and (12.15) may be 
obtained by term-by-term differentiation with respect to ¢ (ef. (5.23) and (4.13 
It may be noted that (12.13), (12.14), and (12.15) are respectively equivalent to the formal 
series (317), (318) and (324) given on pages S4 and 85 of [4] (ef. see. 6 
il a positive 4 rpansions in terms of 
we obtain 


Airy functions In 


$.1] replace ubv ue 4 Then 
| »). . 
lx 9 on L ; 7 ( ‘ 
U( 5 in uty2 ) an F 3 
«| Aiv Ar) > AACS), An ( ait TY B a 12.16 
The f-region of validity ol this CXPansilon, T a includes the segment f | el. fig 
7(b)). Now from [3 
Ze Al f MA \ Bi MA 4 » Al MA 2 7 
Hence substituting (12.16) in (12.1) and equating real parts, we find, with the ud of (12.9 
and the first of (12.3), 
ait 3 riusl (uy z ‘ " A, (C B) U . Boi 
VW (5 u?,uty2 )~ is ( )" | Bi(—uty) > 1. S 12.18 
- 2 es*h / l ae Lu u sand ll j 
valid when ¢t> ; 


The corresponding expansion for the derivative 
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) 


of 


Bi ' C'.(¢ , , ‘ De) 
| — + Bi’ (—ptt) >> - 12.19) 
‘ é u 


ef. (S Ld 
The expansion for the second solution may be obtained by equating imaginary parts. 


This gives 


w(; u?,—uty2 )~ ae (- ya wie) > s u pAi(—a'd) -- oa 


lad rv 3 ( M 


12.20) 


However, this result is only proved in this way for ¢>1, i.e., non-negative values of ¢. This 


is because the error term on truncating the series inside the square brackets at the mth terms is 
Bi ure / Al Mee O(n” Bi’ ure rAd’ wee 1+ | ¢]? ‘~ Mh = 


When ¢ is negative the terms in Bi and Bi’ cannot be absorbed in the terms in Ai and Aji’ be- 
cause they are exponentially larger 
That (12.20) as valid in the wider range ¢> 6 can be proved as follows. In (8.11 


replace u by ue” Then we obtain 


<{ Ai ute y > A 4 Al Mb ¢ > B . ] 12.91) 


bL iT MM 


valid when ¢ € T(} m)and so including t>—1-+-6 (ef. fig. 7 Equation (12.1), with ¢ replaced 
by —?f, shows that the left-hand side equals 


Q-ig— bet -i(de ac W (5 u?,—mty2 )+ik W (5, uty2) | 12.22 


Also from (6.1) and (12.7) we derive 


(] | jae ~24" t¢ ( L if 
ee (we 
y 4 es™*'« / 12.23 
ef. (12.6) and (12.9 The expansion (12.20) may now be obtained by substituting (12.22) 
and (12.23) in (12.21), equating real parts and using the first of (12.3 
lili) a negative. When a is negative the differential equation has no turning points on the 
real axis and we ma confine our attention to the exponential-ty pe expansions 
Krom (2.16), (2.17) and (2.18), we have 
l 9 ) tru ; ri 9 99 
I "| we, uty 2 ) / W( , uly = “*e€ * é 5 U( 5 tv, a pty 12.24 
where 
k= (1+e7™")3—¢7 2° —14O(e- 2), (12.25) 
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and 
e's r (5. : iu® ) Ir (3 = ip ) et (Su? oe (14+ O(n 12.26) 


(ef. (12.2)and(12.3)) 


In (4.3) replace yu, t by wet™, te” 2", respectively. Then we obtai 
, Bist . ( " ul if 
[ ( — = bp yf : pty2)~g ue * >: ( : 2.27 
os f | f | lu 


The region of validity of this expansion, ¢?"S(4m), includes the whole of the real f-axis; this is 


clear from figures 2(b), 34) and 1 
Substituting (12.25), (12.26), and (12.27) in (12.24) in the manner of (i) and (il) and 


determining the branches of the manv-valued functions, Wi obtain finally 


/ ) - u f 
y 9 ‘ ia ) \ 
W (—4u2,uty2) ~ cos (w+ 4x) SG et 
t* 4 | < » l lu 
. UW f 
sin (yré r) >> ’ 12.28 
| rv 
uniformly when o<t<@m. Here f=é¢) and wu,(t) are given by (11.7), (11.9), and (11.11 


and /(u) is calculable from (12.11 
The corresponding expansion for the derivative is 


- . vi : = . h f 
HW } u~,uty2 ~ ia / u)(t- l | sin Ti | rs : [ 


v(t) being given by 11.13) and (11.14 


13. Numerical Examples 


Here we illustrate the use of the Airy-tvpe expansions by two examples. It may be 


remarked that applications would be greatly facilitated by the preparation of suitable numerical 


1 


tables of the functions g(u), [(u), (ft), /& (P—1))* and AY(o), BYO), CO), Die) (s=0,1, 
In the present examples these functions are calculated directly from formulas given in this 
paper 

It will appear from the examples that the early coefficients in the expansions decrease 
in magnitude as s increases. This property, shared with uniform asymptotic expansions for 
Bessel functions of large order ({8], sec.6), makes the expansions well suited to numerical work. 


E-ram ple .. To evaluate W a, r), W’ a, I for a So, 2 7 


We use the expansions (12.18), (12.19), and (12.20 Here 
u—4 t=7/(4V2) = 1.23743 6867 
Hence 
2 ¢* | 05527 596 10 | ) 1) 
and from (12.11), (12.12 
Hu ~ (24/9 1—O.00000 3391—-0_00000 OO06 0.59460 1538 
so that 
wel (py 1.67296 798 129 
Next, from (7.4) we have 
C - f(t = In ¢ f QO.16932 Q5S] 
} | 


6) 


ea 


~é 


nd 


be 
al 


LIS 


Hence 


Krom (4.8 


and from (S.15 


Substitution in 


0.30607 0793. ¢€2—0.55323 6652. c/(t 1)}*#—0.87122 5817 (13.3) 


and $.1] we compute 


s J 0.59504 451 y, 2? 92479 70, y, 26.48264 7, 
a 385 85085 
a l, a ’ a : ’ a - pores’ 
18 - 4608 6 63552 
7 155 Q5095 
Kf § a , db i 
48 ; 1608 6 63552 
8.12 vields 
Ali 0.00649 0. Ble 0.03637 97, B,(c¢)=0.01202. 


Next, we have pi f= 1.94342 840, and by interpolation in table VII of [3] 


0.39524 090, 


0.26162 602. Ai’ ure 0.59097 966 


The content of the square brackets on the right of (12.18 is accordingly 


0.39524 09061 


0.03637 97—0.00004 70+ .. .) 


0.39554 353. 


0.00002 535 4 0.32386 134 


Substituting this value and (13.1), (13.2), (13.3) in (12.18), we obtain 


W(8.7 


Similarly, 


W(s,—7 1.67 


1.67296 798 *& 0.87122 5817 * 0.39554 353 ‘ re 
al 0.14216 485 «10 


$.05527 596 ™ 10 


© ~ 


296 798 0.87122 5817 * 0.26109 871 * 4.05527 596 10 1.54327 57 « 10°. 


For the derivatives, we have first 


From (4.14) and 


ZA 


and thenee from 


Hence 


2” 24 pil rv 1. S7784 306. 
t.16), we compute 
2991942 3, 


a 4, =1.00283 725, Z. 3.52326 59, 4 


0.07833 70. Ce 0.00971, DA(e | tae 0.00760 4 


1.87784 306 * 0.31897 315 


8.7 —————— —— 0.16953 613 « 10 
t.05527 596 10° *« 0.87122 5S14 
1.87784 306 * 0.59419 143 * 4.05527 596 X10" gan eo « 
Tf BO 24 ' 5.19368 20 10°. 


0.87122 5817 


The results of this example agree adequately with the entries given in [4], table I. 
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Kram ple Zs. To evaluate the Hermite polynomial H] I for 7 UJ : by means of the 


expansion (8.11 


Here (ef. (11.23 


ei V2] } 58257 Db6Y, / } \ 21 ) S72S7 156] 
From (6.1) and (4.12 
l/a(u) ~ 2°-%¢° ?] 1+ 0.00198 4127—0.00000 1052 2 11580 14910 
and hence 
or usd (m 2.78291 471 10°. 
Next, from (11.21 
e > eee 't...— t(h—f 0.06286 6477 
} 
Hence 
c QO.15S10 5174, c 0.39762 4413, c/(t l ‘—().90271 123 
When ¢ is negative, real forms of equations (8.12) and (4.8) ar 
A(t > bm (— 5) Ph, ms OBE ‘F< , / 13.4 
where 
J Ti f f | 3D 
In the present example we find 
o l. y, 1.63978 319. J, 91.19031 25. J 516.52363. 
and thence 
Ay lt 1, A(t 0.01024 3, Bye 0.04317 57, B 0.0166 
Next, wtf 1.20344 324 Interpolation in table I of [3], vields 
Ai(urc 0.52582 209, Al’(uic 0.10920 SLO 
Substituting these results in (8.11), we obtain 
bi >, 4V¥2 2.78291 471 10 0.90271 125 * 0.52572 856 =1.32071 865 0) 


and thence using (11.23 
Hi) (4) =2AU( —S> 42 )= 1.25984 22 10° 


The exact value of this polynomial, computed from explicit expressions gi 


of [4], is 1259 84224 
Part 5. Zeros and Associated Values 


14. Zeros of U(a, z), U’ (a, z), Ula, z), and U’ a, z 


l a negative 7—2(] 6 \ (a)* . 24 a). Let the sth real zeros, counted in de- 
scending order away from the point 2,(—a), of the functions 0 (a ('(a U(a,z), and 
L’ ,z be denoted by u F u., u : and ii, ; respectively 

‘ The 


The expansion (8.11) is of exactly the same form as the expansion (4.24) of 


analvsis of section 7 of that paper is therefore immediately applicable, and we find 


' 2a p a) Pp PE | 14.1 
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as u \ aa : ~~. & fixed. Here 


a i a 14.2 


a. denoting the sth negative zero of the function Ai and not to be confused with the a,, of (8.13 


The cot flicients Pp ¢) are given by 
p C U(¢), P C a (Ct ee Pr ¢ Qt. ¢ ‘ C TT ay 4 "; OF oe eg 14.3 


where ¢(¢) is the function inverse to ¢(¢), defined by (11.21), and the a,(¢) are given by rela- 


tions of the form of (7.13) of |S] Explicitly, 
/ of 7) 14 { 
Wt ¢ + . 4) 
f 24(—1)? ' 4g/(#—1e 
The corresponding expansion for the associated value L’(a,u,.,) can be found by sub- 
stitutine (14.1) im (8.15). More convenient expressions for the coefficients are obtained how- 


ever, if we substitute instead in the equation 


UL’ (au +- (Dar {r(; a ) ( - = ) ! 14.5 


obtained from eq (2.15) of [6] and the second of (2.13) in the present paper. This yields the 
required result in the form 


; i I , « Xe a 
l au “rat d r(5-+5 0) } Al’ (a Via I + >) a | 14.6 


where 


and the coefficients 7?,(¢) satisfy the asymptotic identity 


}i+35 : ~1+ {w(g) }? D3 A. 14.8) 
1 r=1 bM 


rs ¢ pent y ¢ “Dp chs Ps ¢ C a y ¢) * Do C). 14.9 


_ 
Thus 


Explicitly, 


14.10) 


¢ tT - 


18(t-—1 96 (¢ 1)?¢ 64¢ 


For the zeros of the derivatives of the parabolic cylinder functions, we may show by 


ul aul g s) +" P +2 a . — 14.11) 
Mh bh 


B=p *aj, 14.12) 


similar analy sis that 


where 


, ; ] | ” 
Jol ¢ t(¢), Ti¢ Biot (oO), Jo\¢ B(O)C(O +4 Teer @ (lls «es (14.13) 


and the 8,(¢) are given by eq (7.22) of [8]. 


The equation corresponding to (14.5) is 


l A,Ua, an) { | Sia )} ‘( —<a- Ue ) ‘( a ) J (14.14) 
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and from it we derive 


tar) 4 ; wa Y,(B 
l(au ~ r( “ent . ) Ai(a,)o(B | Pa ’ 14.15 
u” - — 
in which the coefficients satisfy the identity 
ov ¢ Ole ~ f 4 , G ( 
i+ 2g 7 a >- Lj} > 14.16 
| Lh 4 M Li 
For the function 7 a. the expansions corresponding to (14.1 14.6 14.1] and 
(14.15) are obtained on replacing the symbols U, vw, Ai, and a, by (7, a, Bi, and 6,, respectively, 
(ul a negative : other €1OS Since the zeros of the function Ai are all real and hegative, it 
follows from (8.11) that for all sufficiently large positive values of the parametes the zeros of 
a at < uly 2 in the f-domain T () the unshaded region ol fie 7 il all lie 1) t| ce interval 
1<t<1. This is also true of the zeros of the derived funetion / Fu pty 2 The fune- 
tion [ a uty 2 and its derivative however, have complex Zeros i) this recron corresponding 
to the complex zeros of Bi and Bi’ (ef. [8], appendix 
In the complementary region C'T(0) the distribution of the zeros of / 1? uty2) can 
be investigated by reversion of the expansion (9.7 The situation is more complicated; the 


pattern of the zeros depends largely upon the non-integer part ol 
expansions of zeros in the interval |< ¢< 1 of fixed enumeration 


come within this category 


It should be noted that 


counted away from ¢ . 


(il a complea The extension of the results of this section to general values of a presents 
no difficulty. Wereadily see, for example, that when e<arg p<a2—ethe zeros of / bu, pty 2) in 
the ¢-domain T(arg u) are asvmptotically close to the principal curve \Q illustrated in figure 4. 

15. Zeros of W(a,x) and W’ (a,x 

We denote the sth positive zeros of the funetions Wa wr), W’ (aca Wha r)and W’(a j 
by w Wy», We, and W,.., respectively 

1) a positive Writing wu \ (2a ana 

a lu bh e) h 15.1 
where b il id h are the stl negative zeros ol Bi and Bi’ respectivess ve find b reversion. of 
(12.18) and (12.20 

rT] ws |» a : : | i= 9 
) ) / 1 yy }? . ~ 
2) é i aqiau “ és T u® Bi h y a I ’ | Deed 
, q 3 | D | - 
wr ~ ify q ) | i.) } 
, . : (/ D 
22627" Wawa») ~24e Bi (b/)¢(—B) | 1+ : ) 15.5 
as ( uniformly with respect to Here ¥ O(C), p g q.(€) and Y(€) are the 
functions defined in section 14(i 
With a. B given by 14.2) and (14.12) instead of (15.1) and Bi, 6 replaced by Ai. a. the 
right-hand sides of (15.2 15.5), (15.4), and (15.5) are the corresponding expansions ol Ww 
I . : — — ! , 
- ( , i ad. Ww i and 2 ( 1 ad Ww respectively 


| 
I 


il a negative Let vi \ 2a . and 


t= r(€) (15.6) 
be the function inverse to &(t), defined by (11.7 Then by reversion of the expansion (12.28 
we find 
e T Col T ~— 
lh ~ 2 Mu T + © © «© Ld.¢ 
we LU B 
as a uniformly with respect to s. Here 


T { * S } 15.8) 


and the leading coefficients are given by the formulas 


T ir 567°+- 247 7 2917 2520r 37807 


; 7 — 15.9 
24(7*+-1)° 5/60( 77+ 1 


The expansion of the associated value of the derivative can be obtained by substitution 


of (1loa.4 in the equation 
i ‘(aw T ( 15.10 
zs) 


obtained from eq (2.15) of [6] and (2.19 Thence we derive 


: Q i (7 
W' (aw . 27 4y2(7?+1 l S |: 15.11 
7 
where 
r*+15r?—6 1127 1{897°— 98r°— 349957'* + 804607? — 7020 


18(72+-1)3’ "aX" 22040( 724-1) 


In a similar manner we obtain from (12.29 


where how 


aii 
t) “O67 QATZ OH5l7 2H407 34?0r7 " " 
’ / T ——— 15.15 
24(7 l 5760 (77+ 1 
Finally. from 
| b (dw oe 
Wau T (a a) ( ds ) ’ 15.16 
Wi fina 
. I T 
1 aq. ~ 2 T I > > ] 15.17 
vi 
where 7 1s again given by (15.14). and 
r§—97°+6 1127 1797°— 1427 257257'— 6750077 + 7380 
J , n |: r)\< : : 15.18 
18(77-+-1 23040(77°+-1 
If in the meht-hand sides of (15.7) and (15.11) we take the value (15.14) for r instead 
of (15.8), we obtain the asymptotic expansions of %,., and —W’(a,—@,.,), respectively. Sim- 


larly the right-hand sides of (15.13) and (15.17) with + given by (15.8) are the expansions 


of Ww and 4 ad. we 
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lt may be noticed that some of the numerical coefficients oc urring in the expressions 


(14.4) and 14.10) for the coefficients in the reverted Lirv-type expansions are the same as 
those in the expressions (15.9) and (15.12) pertaining to the reverted exponential-type expan- 
sions. This is no mere coincidence; we may expect that by suitable re-expansion of the reverted 
Airy forms, and subsequent comparison with exponential-t vpe expansions obtained by reversion 
of (5.11), we shall obtain explicit expressions for the coefficients p, and /’, in terms of the 


e, and F, respectively, analogous to the formulas (8.12 Similarly for g, and 4 


16. Numerical Examples 


Evample 1. To compute the smallest positive zero of W (10 and the corresponding 
W’(10,a 

We LUS¢ LS 2 and Ld K with Lu y 20 } 1721 OOH. | o) termmg table \ ot 3], 
we find 

h L-37371i 322. Bi’ (h 0.60195 789 
Hence 
‘ a 1 *b 0.15929 Q 
The value of ¢ is vivel by 
tif nif f | c 0.08477 S64 
Solving this equation by successive approximation, using Newton's rule, we fin 
{= 1.12489 6229= p a 
and thence from (14.4 
a 0.02960 O04 
Hence 
Wio.1™ 2? $.47215 595(1.12489 6229—-0.00007 4001 7.11400 040 
For the derivative, we compute from (14.7) and (14.10 
ly a o(—a)—0.88019 8053. P a 0.00942 2 

Substitution in (15.3) then vields 

= 2” 472 (20 0.60195 789 é x : 

W’ 10.0 i one - - 1—0.00002 356 1.39416 865. 10°°. 
f O.SSOL9 SO53 
By nverse and direct interpolation 1th table II] of $]. we fina 
uw 7.11400 000. W100 1.39416 8610 

The discrepancy in the values of w s undoubted], due to neglect of the third term of the 
asvmptotic series 

[ram ple 2 To compute the fifth positive Zero ot I’ 10. SLT ra the corresponding 


W(—10,2 
We use (15.13) and (15.17 Krom (15.14), (15.6), and (11 
ate value of r satisfies 


7) we see that the appropri- 


r(r?+-1)2?+-In{ r+ (7?+-1 2¢=2(5—4)w/20= 1.49225 6510 


Solving by successive approximation, we find 


r=(0.69390 0720. r?+-1)3—] 


1S and thence from (15.15 


* f(r 0.07269 508, f.(7 0.0332 
d Thus 
nn 1 
yy ~ 92 $.47213 595(0.69390 0720—0.00018 1738—0.00000 0208 ey See 4 38746 276. 

ic 

From (15.18) we compute 

Fy(7 0.00919 198. I(r 0.07882. 
and thence 

, , I 1 i 1 

W(— 10. wey -)~ —24*(20)-*(1.21716 8110)-2(1 + 0.00002 2980—0.00000 0493 


, 0.50972 3297. 


The results of this example agree with values obtained by interpolation in table IIT of [4]. 


References 


} c.. ta Darwin, On Weber’s function, Quart J Mech (ppl Math ie $11 1949 
\ | rdélyi, M Kenne ly, and J L, MeGregor, Parabolic evlinde r functions of large orde! J. Rat. Mecl 
Anal. 3, 459 (1954 
3] J. C. P. Miller, The Airy integral, British Assoc. Math. Tables, Part-Volume B (Cambridge University 
Press, Cambridge, Eng. 1946 
4] J. ¢ P. Miller. Tables of Weber parabolic cylinder functions, London: H. M Stationery Office (1955). 
5] L. M. Milne-Thomson, The calculus of finite differences, London: Macmillan (1933 
F. W. J. Olver, A new method for the evaluation of zeros of Bessel functions and of other solutions of 
second-order differential equations, Proc. Cambridge Philos. Soc. 46, 570 (1950 
7] F. W. J. Olver, The asymptotic solution of linear differential equations of the second order for large values 
of a parameter, Phil Trans Rov soc London, Per \ 247. 307 1954 
8] F. W. J. Olver, The asymptotic expansion of Bessel functions of large order, Phil. Trans. Roy 
Ser. A 247, 328 (1954 
9) F. W. J. Olver, Uniform asymptotic expansions of solutions of linear second-order differential equations 
for large values of a parameter, Phil Trans Roy Soc. Lon lon, Ser. A 250, $79 (1958). 
Math. Sor 37. 339 


Soc. London, 


{10} N. Schwid, The asymptotic forms of the Hermite and Weber functions, Trans. Amer 
1935 

ll] G N \V atson, \ theory ol asvimptotic series, Phil Trans Rov soc London, Ser. A 211, 279 1911) 

12] G. N. Watson, The harmonic functions associated with the parabolic cylinder, Proc. London Math. Soc. 17, 


116 (1918 

13! E. T. Whittaker, On the functions associated with the parabolic cylinder in harmonic analysis, Proce. 
London Math. Soc. 35, 417 (1903 

l4]) bk. T. Whittaker and G. N. Watson, A course of modern analysis (Cambridge University Press, Cam- 


bridge, Eng., 1952 
15] N. D. Kazarinoff, Asymptotic theory of second order differential equations with two simple turning points, 


Arch. Rat. Mech. Anal. 2, 120 (1958). 


Wasuineton, D.C Paper 63B2-14) 


169 





Publications of the National Bureau of Standards 


including papers in outside journals) 


Selected Abstracts 


Diffraction of electromagnetic waves by smooth 
obstacles for grazing angles, J. R. Wait and A. M. 
Conda, J. Re search NBS 63D, No. 2, 18 (1959) 


rhe diffraction of electromagnetic waves by a convex evlin 
irical surface is consid red. Attention is confined primarily 
to the region near the light-shadow boundary The complex- 
ntegral representation for the field is utilized to obtain a 
correction to the Kirchoff theory Numerical results are 
presented which illustrate the influence of surface curvature 
and polarization on the diffraction pattern. Good agreement 
with the experimental results of Bachynski and Neugebauer 
is obtained. The effect of finite conductivity is also con- 
sidered 


Path antenna gain in an exponential atmosphere, 
W. J. Hartman and R. E. Wilkerson, J. Research 
VBS 63D. No. 8, 273 (1959). 


[he problem of determining path antenna gain is treated here 
n greater detail than previously. The method used here takes 
+ into account for the first time the « xpone ntial decrease of the 
gradient of refractive index with height, and a seattering cross 
section inversely proportional to the fifth powel! of the seatter- 
ing angle. fesults are given for all combinations of beam- 
vidths and path geometry, assuming that svmmetrical beams 
re used on both ends of the path and that atmospheric 
turbulence is isotropic, The result as a function of 
both of the beamwidths, in addition to other parameters, and 
thus the loss in gain cannot be determined independently for 
The values of the 


appears 


transmitting and receiving antennas 
8 in gain are generally lower than the previous estimates for 


hich a comparison is possible 


Pattern synthesis for slotted-cylinder antennas, 
J.R. Wait and James Householder, J. Research NBS 
63D, No. 3, 303 (1959 


rhe radiation from a evlinder excited bv an array of axial slots 
\ procedure for synthesizing a given radiation 
pattern is developed with particular attention being paid to a 
lchebyscheff type pattern. Specifving the side lobe level and 
the width of the main beam, the required source distributions 
ire computed for a number of cases The effect of using a 
finite number of slot elements to approximate the continuous 


= discusse d 


source distribution is also considered 


Convexity of the field of a linear transformation, 
A. J. Goldman and M. Mareus, Can. Math. Bul. 2, 
Vo. 1. 15 (1959 

The field F(A 


U’,, is the set of complex numbers ir 
ver the unit sphere in U, The convexity of F(A) (as a 
subset of the z-plane) is proved by a simple inductive argument 


vhich reduces the essential computations to the case n=2 


1 of unitary n-space 
where x ranges 


of a linear transformation 


Graphical diagnosis of interlaboratory test results, 
W. J. Youden, Jnd (dual. Control XV, No. 11, 1 


1959). 


Chis analysis of int rlaboratory test results depends upon the 
ivailability from a number of laboratories of a 
res} ilt for each of two materials The two results from each 
‘aboratory are used to plot a point using the z-axis for one 
material and the y-axis forthe other. The resulting configura- 
tion of the points from the several laboratories permits deduc- 
tions on the prevalence and extent of laboratory bias, the 
presence of sampling variation and the occurrence of blunders 
In addition the graph provides an estimate of the precision of 
the test proced ire results, 


single test 
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The calculation of the field in a homogeneous 
conductor with a wavy interface, J. R. Wait, Proc. 
IRE 47, 1155 (1959). 

The field at any depth in a homogeneous conductor with a 
wavy interface is calculated. For purposes of illustration, the 
field above the interface is taken to be a uniform plane wave 
traveling in the horizontal direction, 


Lower bounds for eigenvalues with application to 
the helium atom, N. W. Bazley, Proc. Nat. Acad. 
Sea. U.S. 15. No. 6. 850 (1959). 


Let A be a self-adjoint operator with domain D in a Hilbert 
A=A’'+A where A is self-adjoint and A’ is 

The eigenvalue problem for A, whose solu- 
bounds. If u 


space Suppose 


positive definite 


tion we assume known, gives rough lower 
=I, , k) are k discrete eigenvectors of A and if p;=(a’)" 
uitti= 1 , k) exist then one can substantially improve the 


lower bounds, The theory is applied to helium atom operator. 
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